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Measurements on X-Ray Production and Absorption in the Range 0.7 to 2.5 Mv 
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In the voltage range from 0.7 to 2.5 Mv, the intensity 
and the total mass absorption coefficient of the heterogene- 
ous radiation produced in a thick gold target have been 
measured as a function of tube voltage, angle with respect 
to the electron beam, and thickness and material of the 
filter. Curves giving the dependence of absorption coeffi- 
cient on voltage are presented for lead, tin, copper, alumi- 
num, carbon, and water as absorbers. A method is de- 
scribed for calculating the absorption coefficient of the 
heterogeneous radiation produced at a given tube voltage 
and subjected to a given amount of filtration, utilizing 
accepted values for the absorption coefficient of homogene- 


ous radiation. Measured and calculated values of the 
absorption coefficients for heterogeneous radiation are in 
good agreement. Curves giving the dependence of intensity 
on angle are in qualitative agreement with approximate 
theory. Curves giving the dependence of absorption coeffi- 
cient on angle and the dependence of intensity in the 
forward direction on voltage are presented. A method for 
calculating the efficiency of x-ray production from these 
data is described and used to obtain results in satisfactory 
agreement with theory in the voltage range from 0.90 
to 2.35 Mv. 





I. INTRODUCTION 


HIS research is concerned with the intensity 
and quality of the continuous radiation 
produced at a thick target of high atomic number 
by bombardment with high energy electrons, and 
with the absorption of this radiation by various 
materials.! Such information has application in 
high voltage x-ray research and therapy, and 
makes possible a limited comparison with existing 
theory. 

Most previous precision measurements of ab- 
sorption coefficients have been made with homo- 

* Now in the Physics Department, University of Notre 
Dame, Notre Dame, Indiana. 

** Now in the Radiation Laboratory, Massachusetts 
Institute of Technology, Cambridge, Massachusetts. 

*** On leave from the Physics Department, New York 
University, and now at the Frankford Arsenal, Frankford, 
Pennsylvania. 

1 Rough preliminary results were presented earlier: L. 
C. Van Atta and D. L. Northrup, Am. J. Roent. and 


Rad. Therapy 51, 633 (1939); present authors, Phys. Rev. 
59, 688A (1941). 
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geneous radiation either at low voltage where 
homogeneous x-rays are obtainable by crystal 
reflection, or at high voltage as the homogeneous 
gamma-rays from radioactive materials. Such 
measurements have been made in various ma- 
terials? up to 0.5 Mv with a probable error of one 
percent, up to 2.5 Mv with a probable error® of 
four percent, and in a few elements at isolated 
higher voltages‘ with less accuracy. Experiment 
and theory concerning the absorption of homo- 
geneous radiation are in satisfactory agreement 
throughout the voltage range involved. 
Measurements of the intensity and quality of 
the heterogeneous radiation produced at a thick 
target have been made up to 0.5 Mv, and to a 
limited extent up to 1.0 Mv. In the low voltage 


‘- Read and C. C. Lauritsen, Phys. Rev. 45, “ os 
Y. Chao, Phys. Rev. 36, 1519 (1930) ; . P. 
Pin Proc. Roy. Soc. 128A, 345 (1930). 
4L. A. Delsasso, W. A. Fowler, and C. C. Lauritsen, 
Phys. Rev. 51, 391 (1937). 
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Fic. 1. General view of the x-ray target, collimator, and thimble chamber, 
and the structural frame supporting the collimator. 


range the efficiency of production of x-rays, the 
angular distribution of the radiation from thin 
targets, and the absorption coefficients in various 
materials have been measured with considerable 
precision. In the region between 0.5 and 1.0 Mv, 
where characteristic radiation no longer plays an 
important part, much remains to be done. The 
variation of the intensity of heavily-filtered 
radiation with angle has been determined at 0.5 
Mv.*® The dependence of intensity on thickness in 
various absorbers has been measured for the 
heterogeneous radiation produced at voltages up 
to 1.0 Mv.® 
Il. APPARATUS 


The general scheme is to make intensity meas- 
urements in a severely collimated x-ray beam 
with an air-equivalent thimble chamber. The 
collimation serves to define the angle with respect 
to the direction of the original electron beam, to 
exclude x-rays scattered in the absorber, and to 
exclude the background of x-rays scattered from 
various parts of the room. By varying the direc- 


5 K. E. Corrigan and B. Cassen, Am. J. Roent. and Rad. 
Therapy 37, 811 (1937). 

®R. Phillips and G. S. Innes, Brit. J. Radiol. 11, 498 
(1938) ; J. C. Hudson, Radiol. 32, 530 (1939). Also Folsom, 
Stone, Rose, Bouwers, and Hermann and Jaeger. 


tion of the axis of the collimator with the x-ray 
target as a center, the intensity or the absorption 
coefficient of the radiation can be measured as a 
function of angle. Figure 1 shows the x-ray 
target, the collimator, the thimble chamber, and 
the aluminum structural frame which supports 
the collimator. 

The collimator was subjected to two perform- 
ance tests. Absorption measurements in lead for 
thicknesses up to 10 cm were made with 2.5- and 
1.3-cm diameter defining diaphragms. The re- 
sulting curves (J/Jo versus absorber thickness) 
were identical, showing that the 2.5-cm diameter 
diaphragm provided sufficient collimation. As the 
thickness of lead absorber was further increased 
to 15 cm the x-ray intensity at the thimble 
chamber continued to decrease exponentially, 
which indicated that the background radiation 
inside the collimator was negligible. 

The thimble chamber was adopted for the 
measurement of x-ray intensity on the basis of 
considerable previous experimental’ and _ theo- 


70. Glasser and L. Rovner, Am. J. Roent. and Rad. 
Therapy 36, 94 (1936); L. S. Taylor, G. Singer, and A. L. 
Charlton, Am. J. Roent. and Rad. Therapy 41, 256 (1939) ; 
and others. 
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retical*® work. This work has shown that the 
thimble chamber is the only practicable method 
for measuring radiation above 0.4 Mv, that 
measurements made by different observers with 
chambers of different geometries and different 
“air-equivalent’’ wall materials are accurately 
comparable, and that the corrections for wall 
effect necessary to reduce reading to roentgens, 
according to the accepted definition, are small 
and calculable.* An ‘‘air-equivalent” wall ma- 
terial is one in which the x-ray beam has the same 
equilibrium complement of electrons as it has in 
air. Theoretically the wall thickness should be at 
least as great as the range of secondary electrons 
in the material® if the chamber is to be used in a 
collimated beam, but experimentally the wall 
thickness can be somewhat less than this. In the 
case of carbon which is the element most nearly 
air-equivalent, the correction factor to obtain 
equivalent roentgens from the actual readings is 
1.023 in the case of gamma-rays from radium in 
equilibrium with its decay products, and is 1.027 
in the case of the total radiation at 0.6 Mv. Thus 
the x-ray intensity measured by a carbon-wall 
thimble chamber in the energy region from 0.6 to 
2.5 Mv may be written r= 1.025(J/V) roentgens 
per min. where J=ionization current (stat- 
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Fic. 2. Detailed view of the carbon-walled thimble 
chamber with its mounting strap. 


8G. C. Laurence, Can. J. Research 15A, 67 (1937). 
*L. H. Gray, Proc. Roy. Soc. 156A, 578 (1936). 
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Fic. 3. Absorption curves in lead for the radiation 
produced at various voltages. 
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FILTER THICK 


coulombs per min.) and V=volume of cham- 
ber (cc). 

The thimble chamber used in these experiments 
is shown in cross section in Fig. 2, with a mount- 
ing strap for fastening the chamber to the lower 
end of the collimator. The volume of the chamber 
is 1.94 cc and the wall thickness is 2.4 mm on the 
sides. High voltage is applied to the carbon 
thimble by means of a wire leading to the side 
binding post. Ionization current is collected on 
the central aluminum wire and conducted to an 
amplifier. An aluminum tube connecting to the 
cable shield acts as a guard ring to reduce the 
effective leakage current. Saturation voltage for 
high x-ray intensities proved to be 270 v. This 
thimble chamber was subjected to several per- 
formance tests which proved it to be entirely 
satisfactory for the measurements involved. 

The ionization produced in the carbon thimble 
chamber was measured with a current amplifier 
rather than by a condenser discharge method, in 
order to have a continuous reading of x-ray in- 
tensity and to avoid the necessity of maintaining 
constant conditions. The sensitivity ratios for the 
various ranges were measured with a constant 
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x-ray beam under actual conditions of use and 
were found to agree within the accuracy of 
reading with the ratios calculated from the input 
resistor values. 

The construction of the target is indicated in 
Fig. 1. It is evident that much greater absorption 
must be expected in the neighborhood of 90° and 
that this effect must be considered in connection 
with angular distribution curves. Gold was 
chosen as the target material because it is the 


element of largest atomic number practicable for . 


such use. The thickness used was required both 
for strength and for stopping the high speed 
electrons. An earlier target! of lead electroplated 
onto copper was not satisfactory because of a 
tendency of the lead to evaporate away from the 
target spot. 

For the absorption coefficient measurements 
the absorber consisted of disks of the pure ma- 
terial 10 cm in diameter and each of such thick- 
ness as to reduce the intensity by 30 to 50 per- 
cent. In the case of water, distilled water of 
known temperature was used in a 15-cm diameter 
glass cylinder. The total thickness of absorber, in 
addition to a lead pre-filter, varied from 10 cm in 
the case of lead to 30 cm in the case of water. 


Ill. MEASUREMENT OF ABSORPTION 
COEFFICIENTS 


For the measurement of absorption coefficients, 
thimble chamber intensity readings in the colli- 





"a | 


k\ 


7 sion 
. 


VERSUS FILTER THICRNESS 














MASS ABSORPTION COEFFICIENT 
































aos Pam aA 
ave ss id 
se al fs — 16 a 
eee Pan tended 
ase J 
im 4 4 J a s 





FILTER THICKNESS § (CM OF LEAD) 


Fic. 4. Mass absorption coefficient in lead .versus filter 
thickness for the radiation produced at various voltages. 





ATTA, AND MYERS 


mated beam were taken as a function of thickness 
of absorber placed above the defining diaphragm. 
Comparable relative intensity readings were ob- 
tained by dividing each intensity reading by the 
corresponding electron current, averaging the 
results obtained with increasing and decreasing 
filter thickness, and dividing successive intensity 
values by the initial value. When such relative 
intensity (I/Jo) values are plotted against filter 
thickness on semi-logarithmic paper, the slope of 
the resulting curve provides a measure of the total 
mass absorption coefficient, um=In (1o/Z)/px. 
Curvature in the initial part of these curves is an 
indication of a decreasing absorption coefficient 
because of the selective absorption of the softer 
components in the heterogeneous radiation. After 
sufficient filter has been used the wave-length 
distribution of the radiation in the beam becomes 
so narrow that the corresponding absorption 
coefficient is approximately constant. The re- 
sulting straight line portion of the logarithmic 
curve permits an accurate determination of the 
final constant value of absorption coefficient. 

Typical absorption curves are shown in Fig. 3, 
which gives the absorption in lead of the radia- 
tion produced at various voltages. It appears that 
the absorber thickness necessary for establishing 
the equilibrium condition increases with voltage, 
but that the final slope is constant over a wide 
range of absorption. This is more evident in 
Fig. 4, in which the mass absorption coefficient 
obtained from the slope of the logarithmic curves 
is plotted against filter thickness for the radiation 
produced at several voltages. From Fig. 4 we see 
that the initial absorption coefficient at each 
voltage is somewhat more than twice the final 
absorption coefficient and that the filter thickness 
necessary for obtaining the equilibrium condition 
is roughly proportional to the voltage. 

In Fig. 5 is given the variation with voltage of 
the final value of absorption coefficient for 
several absorbing materials. The curves are 
drawn through the experimental points without 
regard to theory. In the case of absorbers other 
than lead a sufficient thickness of lead pre-filter 
was used to establish approximately the equilib- 
rium condition. 

The data presented in Fig. 5 represent the 
summary of approximately 1600 intensity read- 
ings on 83 absorption curves. For each absorption 
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Fic. 5. Experimental values for the mass absorption coefficient in several materials of 
heavily filtered heterogeneous radiation plotted as a function of the tube voltage. 


curve an estimated probable error in the slope can 
be obtained by considering the spread in the 
curves which might reasonably be used to repre- 
sent the points. Such estimated probable errors 
for individual curves (i.e., for individual determi- 
nations of absorption coefficient) ranged from 0.3 
percent to 4 percent with an average value of 1.2 
percent. 

It is evident from Fig. 5 that in four of the 
curves (H,O, C, Al, and Cu) the same type of 


absorption is present, whereas in two of the 
curves (Pb and Sn) there is an additional absorp- 
tion which becomes important at low voltage. In 
the higher part of the voltage range the order of 
increasing mass absorption coefficients (excepting 
Pb) is Sn, Cu, Al, C, and H,O. These results are 
in agreement with theory in which the total ab- 
sorption of energy from the x-ray beam is con- 
sidered to be due to true absorption and scattering 
by free electrons (Compton absorption and 
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scattering),!° absorption by bound electrons 
(photoelectric absorption), and absorption due 
to the production of positron-electron pairs (pair 
production). 

Photoelectric absorption and pair production 
are appreciable in this voltage range only for lead 
and tin. Photoelectric absorption explains the 
steeper rise in the lead and tin curves of Fig. 5 at 
the low voltage end, and is appreciable through- 
out the entire voltage range in the case of lead. 
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Fic. 6. Distribution in voltage of unfiltered (1) and 
filtered (2, 3, 4) radiation assuming Kramers’ law and 
accepted values for lead absorption coefficients. 


With the exception of lead, the absorption coeffi- 
cient in the middle part of the voltage range is 
almost entirely due to Compton effect and there- 
fore should be proportional to the number of 
electrons per gram of absorber, as is found 
experimentally to be the case. At the high voltage 
end, pair production" is responsible for about 12 
percent of the absorption in lead and 7 percent in 
tin. 

It appears therefore that the measured ab- 
sorption coefficients for heterogeneous radiation 


QO. Klein and Y. Nishina, Zeits. f. Physik 52, 853 
(1929). 

1H. R. Hulme, J. McDougall, R. A. Buckingham, and 
R. H. Fowler, Proc. Roy. Soc. 146A, 83 (1934). M. T. 
Jones, Phys. Rev. 50, 110 (1936). 

2H. Bethe.and W. Heitler, Proc. Roy. Soc. 146A, 83 
(1934). 

18 For the calculation of absorption by pair production, 
we are indebted to Dr. R. D. Evans, ‘Introduction to 
the Atomic Nucleus,”” M.I.T. course notes. 
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are in general qualitative agreement with ac- 
cepted theory for homogeneous radiation. In 
order to make a quantitative comparison, it is 
necessary to know the distribution in wave- 
length of the heterogeneous radiation. Kramers’ 
theory" yields the relation 


Ivy =constant- V?(V»— V) 


where V)=tube voltage. This relation was used 
to obtain curve (1) in Fig. 6 for the distribution in 
voltage of the unfiltered radiation produced 
either at 1.0 or 2.0 Mv. The change in the shape 
of this curve as the result of lead filtration can be 
calculated by using accepted values for the ab- 
sorption coefficient of homogeneous radiation in 
lead at the different voltages. Curves (2) and (3) 
were obtained by assuming the thicknesses of 
lead pre-filter actually used at 1.0 and 2.0 Mv, 
respectively, in measuring the absorption coeffi- 
cients of other materials. From these curves it 
appears that the distribution of the radiation 
after pre-filtering was very similar in the two 
cases. Curve (4) gives the distribution of the 
radiation produced at 2.0 Mv and filtered 
through 11.5 cm of lead. From curves (1), (3), 
and (4), together with accepted values for the 
absorption coefficient of homogeneous radiation 
in a given material, it is possible to calculate the 
absorption coefficient in that material of the 
heterogeneous radiation produced at 2.0 Mv and 
filtered through 0, 3.8, and 11.5 cm of lead. The 
results of these calculations for lead and alumi- 
num absorbers are compared with the experi- 
mental values in Table I. ; 

These calculations show the same effect of lead 
filter on heterogeneous radiation as was observed 
experimentally, i.e., that the initial high absorp- 
tion coefficient in lead is reduced to roughly half- 
value and then becomes approximately constant 


TABLE I. Calculated and measured values for the mass 
absorption coefficient in lead and aluminum of the hetero- 
geneous radiation produced at 2.0 Mv after various 
amounts of lead filter. 











Thickness of Mass absorption coefficient (cm*/g) 
lead filter Lead Aluminum 
(cm) Calculated Measured Calculated Measured 
0 0.090 0.059 
3.8 0.054 0.056 0.052 0.052 
11.5 0.050 0.048 








4H. A. Kramers, Phil. Mag. 46, 863 (1923). 
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Fic. 7. Variation of intensity with angle for the unfiltered and filtered radiation at three values of tube voltage. 


as the filter thickness is further increased. The 
absorption coefficient at small filter thicknesses 
could not be determined accurately because it 
changed so rapidly with added absorber. There 
was not sufficient intensity left after 11.5 cm of 
lead filter to permit an accurate measurement of 
the absorption coefficient of the resulting radia- 
tion. However, for lead thicknesses between 3.8 
cm and 11.5 cm no appreciable change in absorp- 
tion coefficient was observed. The calculated and 
measured values for both lead and aluminum 
absorbers after 3.8 cm of lead filter are in very 
satisfactory agreement. This agreement tends to 
confirm the values for the absorption coefficient 
of homogeneous radiation and to indicate that 
the distribution function given by Kramers’ 
theory represents a reasonable approximation in 
spite of its theoretical limitations. Furthermore 


it tends to justify the method outlined above for 
calculating the absorption coefficient of the 
heterogeneous radiation in a severely collimated 
beam of x-rays produced at a given tube voltage 
and subjected to a given amount of filtration. 


IV. MEASUREMENT OF ANGULAR DISTRIBUTION 


To measure the x-ray intensity as a function of 
angle with respect to the forward direction of the 
electron beam, the collimator was pivoted in its 
support to swing about the x-ray target as a 
center. The thimble chamber in the collimated 
beam then provides an approximately true meas- 
ure of the intensity as a function of angle. Having 
the collimated beam also permits insertion of 
absorbing material into the beam to obtain a 
measure of the quality of the radiation as a func- 
tion of angle. 











With the tube voltage held constant the 
collimator was set at various angles, and the 
intensity and electron current were recorded. The 
intensity readings were divided by the corre- 
sponding electron currents and by the specific 
intensity in the forward direction to obtain com- 
parable relative intensity readings. This proce- 
dure was followed with no filter and with a heavy 
lead filter for three values of tube voltage. 

The results of the angular intensity measure- 
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Fic. 8. Relative softness of the radiation produced at 
several voltages as a function of angle, where relative 
softness is defined as the ratio of unfiltered to filtered 
intensity. 


ments are presented in Fig. 7. These curves show 
the strong predominance of the radiation in the 
forward direction. The intensity ratio between 0° 
and 90° varies from 4.2 for the unfiltered radia- 
tion at 0.90 Mv to 7.3 for the filtered radiation at 
2.35 Mv. No attempt was made to pass these 
curves through the 90° experimental points be- 
cause of the excessive absorption introduced by 
the target in that direction. 

In Fig. 7, the unfiltered and filtered relative 
intensity values are compared for three values of 
tube voltage. This comparison shows that the 
radiation is considerably softer at larger angles. 
The softening of the unfiltered radiation with 
increasing angle is an effect which increases with 
voltage, as is shown by Fig. 8 in which the ratio 
of unfiltered to filtered relative intensities is 
plotted against angle for three values of tube 
voltage. The ratios have been adjusted to unity 
at the angle of hardest radiation, rather than at 
0°. These curves of an arbitrarily defined quan- 
tity are presented primarily because of the way in 
which they emphasize the angle of maximum 
hardness of the radiation. 
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A detailed examination of the curves of Fig. 7 
at small angles reveals an interesting fine struc- 
ture. There is a maximum in the curves which 
becomes more pronounced and which moves in to 
smaller angles as the voltage is increased. This 
effect is so small at 0.90 and 1.63 Mv that the 
small angle points on these curves were taken 
with great care. The additional maximum at 0° 
on the 0.90- and 1.63-Mv curves can be obtained 
by combining two slightly overlapping angular 
distribution patterns of conventional shape, and 
is therefore assumed to be of no especial theo- 
retical interest. 

The angular positions of the intensity maxi- 
mum at the several voltages are given in Table II 
with an estimated probable error of 0.5°. Exami- 
nation of the curves shows that the maximum 
moves in to smaller angles as the voltage is in- 
creased, but for a given voltage the harder com- 
ponents of the radiation have their maxima at 
larger angles. If these curves are thought of as 
the combination of two slightly overlapping 
angular distribution patterns on opposite sides of 
the axis, resolving the lobes would increase 
appreciably the angle of the maximum. 

The dependence of x-ray quality on angle for 
the heavily filtered radiation from a thick target 
was obtained by measuring the final value of the 
mass absorption coefficient in lead at a number 
of angles with a fixed tube voltage. The depend- 


TaBLe II. Angular position of the intensity maximum at 
several tube voltages. 








Positions of 
Positions of the intensity maximum hardness 





Unfiltered Filtered Sommerfeld maximum 
Voltage radiation radiation theory (Fig. 8) 
0.90 9.0° | yg * ig 
1.63 3.8° 4.5° 8° ha 
2.35 2.8° 3.3° 6° 3.7° 








ence on angle of this absorption coefficient is 
given in Fig. 9 for the radiation produced at 1.09 
and at 1.63 Mv. A pronounced dip in the curves 
at about 7° is followed by a regular rise as the 
angle is further increased. These curves illustrate 
the advantage of working in the forward direction 
rather than at 90° when penetrating radiation is 
required, and the additional advantage of work- 
ing at a small angle from the forward direction if 
a small portal is used. 
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There is no available theory strictly applicable 
to the problem of calculating the angular distri- 
bution of the radiation produced in a thick target 
of high atomic number at high voltage. This is 
not surprising in view of the extreme difficulties 
imposed by these conditions. Sommerfeld’s clas- 
sical theory'® is relativistic but assumes a thin 
target and complete stopping of the electron 
without change in direction. His wave mechanical 
theory'® is non-relativistic and therefore not 
applicable. Later theories!’ are too complicated 
to justify calculation in that they are approxi- 
mations for elements as heavy as gold. In spite of 
its limitations, therefore, Sommerfeld’s classical 
theory has been used to calculate the positions of 
the intensity maximum in the angular distribu- 
tion curves (last column, Table IT). Actually the 
angular position of the hardest radiation (mini- 
mum in Fig. 8) is in better agreement with 
Sommerfeld’s predicted angle than is the position 
of the intensity maximum. In general, theory 
predicts that the maximum will move toward 
smaller angles and will become more pronounced 
as the voltage is increased, but that for a given 
voltage the harder components of the radiation 
will have their maxima at larger angles. In these 
respects theory and experiment are in qualitative 
agreement. 


V. EFFICIENCY OF X-RAY PRODUCTION 


In all the measurements described previously 
the results have not depended on the absolute 
calibration of the thimble chamber, and not 
sensitively upon the variation of this calibration 
with x-ray quality. The results presented in 
Fig. 10, however, depend on the carbon thimble 
chamber for the absolute measurement of in- 
tensity in r-units. The x-ray intensity in the 
forward direction is plotted on double logarithmic 
paper as a function of voltage for the unfiltered 
radiation (except for the 2.3-mm lead equivalent 
thickness introduced by the target) and for 
radiation filtered through 2.5 cm of lead. The 
unfiltered radiation increases approximately as 
the third power of the tube voltage. The filtered 
radiation shows a more uniform increase as the 


18 A, Sommerfeld, Physik. Zeits. 10, 969 (1919). 

16 A. Sommerfeld, Ann. d. Physik 11, 257 (1931). 

”W. Heitler and F. Sauter, Nature 132, 892 (1933); 
F. Sauter, Ann. d. Physik 20, 404 (1934). 


4.7 power of the voltage. This emphasizes the 
importance of having high voltage for producing 
radiation which is to be used through consider- 
able absorbing material. 

From a consideration of the x-ray intensity 
available at high voltage and the low value of the 
absorption coefficient after moderate filtration, it 
is evident that radiographical examination of 
very thick specimens, such as heavy steel 
castings, becomes feasible. If the target-to-speci- 
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Fic. 9. Dependence on angle of the final value of the 
mass absorption coefficient in lead for the radiation 
produced at 1.09 and 1.63 Mv. 


men distance is made large, sufficient pre-filter is 
used, and stray radiation is carefully eliminated, 
good resolution should be obtained for thickness 
as great as 12 in. of steel. 

The efficiency of x-ray production may be 
defined as the ratio of the power radiated from 
the target as x-rays to the power received at the 
target as the result of electron bombardment. 
The power radiated as x-rays at a given tube 
voltage and beam current can be determined from 
a knowledge of the x-ray intensity (in roentgens) 
in all directions from the target, provided that 
the proper correction be made for true absorption 
in the target. 

The rate at which energy passes through unit 
area at a distance r from the target and at an 
angle @ with respect to the forward direction of 
the beam is given in terms of the reading of an 
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Fic. 10. Dependence of absolute intensity in the forward direction on voltage for 
unfiltered and filtered radiation. 


air-wall thimble chamber at that point by 
Pia= fi Idr=TI,,o/u roentgen cm/min. 


where u is the linear coefficient for true absorption 
in air of the heterogeneous radiation incident 
upon the thimble chamber. If we accept a value 
of 45.5 ev per ion pair as the energy required for 
x-ray ionization in air,!® we may write this ex- 
pression in terms of a more conventional unit for 


1” The generally —- or value of 32 ev per ion pair is 
2 based on Geiger’s early value of 33 ev 
obtained from alpha- 4%. ionization experiments (Proc. 
Roy. Soc. 82A, 486 (1909)). The value of 45.5+1.0 ev 
obtained from the excellent work of J. F. Lehmann and 
T. H. Osgood on ionization in air by electrons in the energy 
range 200 to 1000 ev (Proc. Roy. Soc. 115A, 609 (1927)) 
more properly applies to the case of x-ray ionization. We 
are indebted to Dr. M. S. Livingston, who brought these 
references to our attention as well as his own confirming 
calculation based on ionization by radium gamma-rays. 


radiation density 
P,,g=2.53 X10-"°T,, o/u watt/cm?. 
The total rate at which energy passes through a 


sphere of radius 7 is 


P,=1.59X10-%r? | (J,,6/u) sin 6d@ watt. 
0 


If we introduce J;, the x-ray intensity in the 


TABLE III. Efficiency of x-ray production: a comparison 
of experiment with theory. 








Experiment Theory 
Target 
Voltage Correction 
(Mv) | Is Ig sin 6d0 K Factor ” > 





2.35 | 750 0.52 3.08X10-§ 1.22 10. 
1.63 | 240 0.61 3.30K10~§ 1.35 5S. 
0.90 44 0.68 3.48105 2.00 3 
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forward direction expressed in roentgens per min. 
and milliampere at one meter from the target and 
Ip, the intensity at angle @ relative to that in the 
forward direction, and if we recognize the fact 
that the coefficient for true absorption in air does 
not depend appreciably on angle, we may write 


P= 1.59X10-*(I,/u) Ig sin 6d@ watt/ma. 
0 
Then the efficiency of x-ray production is 
n=1.59(1[;/p vf Ig sin 6d@ percent 
0 


where V is the tube drop in volts. 

Experimental data already presented make it 
possible to calculate the efficiency at three 
voltages: 2.35, 1.63, and 0.90 Mv. Values of I; 
are obtained from Fig. 10 and values of the 
integral are calculated from the dashed curves of 
Fig. 7. The values of u, the linear coefficient for 
true absorption in air of the heterogeneous radia- 
tion incident upon the thimble chamber, were 
calculated as weighted averages from the corre- 
sponding coefficients for homogeneous radiation. 
From curves similar to those shown in Fig. 6 it is 
possible to obtain directly factors to correct for 
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true absorption of x-ray power in the target. The 
ratio of the area under the Kramers’ distribution 
to that under the reduced distribution gives the 
factor by which the power has been reduced in 
passage through the target. The correction is not 
serious except at the lowest voltage. Values of J,, 
the integral, uw, the correction factor, and 9 are 
presented in Table ITI. 

The theoretical values for the efficiency given 
in the last column were computed from curves 
calculated from Bethe and Heitler’s theory by 
Arcimovic and Chramov.”® In view of limita- 
tions in the theory when applied to targets of high 
atomic number, the agreement between theory 
and experiment throughout this voltage range 
may be considered quite satisfactory. 

The writers are indebted to Dr. G. G. Harvey 
for helpful discussions and to Mr. Zigmond 
Wilchinsky, Mr. W. B. Nowak, and Mr. Howard 
Rowland for much capable assistance in making 
measurements and calculations. The research 
was supported in part by grants from Research 
Corporation and the Carnegie Corporation of 
New York. 


*L. A. Arcimovic and V. A. Chramov, Comptes 
Rendus Acad. Sci. U.S.S.R., 18, 415 (1938). 





Errata: Quantization of Molecules, Inter- and Intramolecular Forces 
[Phys. Rev. 63, 309 (1943) ] 


KASIMIR FAJANS AND THEODORE BERLIN 
Department of Chemistry, University of Michigan, Ann Arbor, Michigan 


HE sentence on page 310, left column, lines 8—6 from the bottom should 
read: “‘In Fig. 1A, two Li* cores are placed at the internuclear distance 


0.74A of He (n=I).”’ 


On page 311, left column, line 7 from the bottom should have n=3 instead 


of n=IITI. 





Erratum: The Carbon Arc in Oxygen for the Spectrochemical Determination of Potassium 
[Phys. Rev. 63, 322 (1943) ] 


L. T. STEADMAN 
Department of Radiology, The University of Rochester, School of Medicine and Dentistry, Rochester, New York 


HE last paragraph, beginning line 8 from the bottom, should read: ‘‘To the 

cathode is added 0.01 ml of serum, which amount in the human normally 
contains 33.0 gamma (micrograms) Na, 2.0 gamma K, 1.0 gamma Ca, 0.27 
gamma Mg, and 1.3 gamma P. Also, 2.5 gamma Rb is added as the internal 
standard for K and 100 gamma Cd as the internal standard for all the other 
elements. The spectrum lines measured are Na 2680.3A, K 4044.1A, Ca 3006.9A, 
Mg 2783.0A, P 2535.6A, Rb 4201.8A, and Cd 2677.6A.” 
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Schwinger’s modification of the Mgller-Rosenfeld theory 
of nuclear forces based on a mixture of pseudoscalar and 
vector meson field is investigated in the strong coupling 
theory. The isobar separation turns out to be } of its value 
in the pseudoscalar theory. The interaction energy between 
two nucleons is calculated, and in the most interesting case 
of the symmetrical theory, it is the same as the corre- 
sponding expression in the weak coupling theory with the 
spin and isotopic spin vectors replaced by the e* vectors of 
Pauli and Dancoff. A classification of the states of the two- 
nucleon system is made, and the lowest state for the 
deuteron is found to be the triplet state. An estimate is 
made for the values of the constants involved, and we find 
that a suitable choice of the coupling constant and the 


masses of the mesons can be made in such a way that the 
binding energy and the quadrupole moment of the deuteron 
agree with the observed values, and such that the condi- 
tions for small source (compared to the Compton wave- 
lengths of the mesons) and for small effect of the higher 
spin states on the ground state are fulfilled. However the 
theory gives for the magnetic moment of the deuteron a 
value only a few percent of the observed value, and ac- 
cording to this theory highly charged nuclei would be 
unstable. We therefore conclude that the strong coupling 
theory, based on the assumption of an extended source, 
should be abandoned in favor of a weak coupling theory, 
based on a point source with the singularities of its field 
eliminated by means of a subtraction formalism. 





































_ I. INTRODUCTION 


HE usual theories of nuclear forces derived 

from one type of meson field give rise, 
both in the weak and strong coupling approxi- 
mations, to tensor forces with inadmissible r~* 
singularity for small r, which consequently has 
to be cut off at a certain arbitrary radius. 
Mgller and Rosenfeld! showed that the singu- 
larity can be removed by taking a mixture of 
pseudoscalar and vector meson fields with the 
same coupling constant. They, however, also 
took the masses of the two types of meson equal, 
and found that the tensor force vanished in the 
first approximation. Schwinger? then pointed 
out that the tensor force can be retained with 
only an admissible r~' dependence at small 
distances if the masses are taken unequal, and 
that the right sign for the quadrupole moment 
is obtained if the vector meson has the larger 
mass. This is in agreement with the hypothesis 
that the vector meson is highly unstable and 
that it is responsible for the 8-decay in the 
nucleus. 

It may be argued that the introduction of the 
vector meson field is just another way of re- 
moving the inadmissible singularities by an 
arbitrary constant, and in a much more compli- 

1C. Mgller and L. Rosenfeld, K. Danske vidensk. Selsk. 


17 (1940). 
2 J. Schwinger, Phys. Rev. 61, 387A (1942). 
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cated way at that. The introduction of a new 
particle is admittedly undesirable, but we feel 
that further investigation in this direction is 
justified for the reason that this may be a first 
step in the direction of a more far-reaching 
theory in which the vector meson enters, not as 
a new independent particle, but as an excited 
state of the pseudoscalar meson. 

We have investigated in this paper the mixed 
meson theory in the strong coupling approxi- 
mation. The calculations of Section II for the 
isobar separation follows closely the development 
for the pseudoscalar meson given by Pauli and 
Dancoff.* The theory of nuclear forces in the 
strong coupling approximation has been de- 
veloped by Serber and Dancoff‘ who treated 
the charged scalar and the neutral pseudoscalar 
meson fields, and the calculation in Section III 
follows their procedure. Hence the computations 
in these sections are not given in detail. However 
we wish to point out that in contrast to Serber 
and Dancoff, we do not consider the nucleons to 
be at rest, and consequently we do not discuss 
here the small oscillations around the minimum 
of the potential energy which occurs for suffi- 
ciently small separation of the nucleons. The 

*W. Pauli and S. M. Dancoff, Phys. Rev. 62, 85 (1942). 
This paper will henceforth be quoted as P-D. 

4R. iter and S. M. Dancoff, Phys. Rev. 63, 143 
(1943). This paper will henceforth be quoted as S-D. 


We wish to thank these authors for letting us see their 
calculations before publication. 
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results of these sections are applied in Sections 
IV and V to the deuteron and the heavy nuclei, 
respectively. 


II. ISOBAR SEPARATION 


We describe the pseudoscalar part by the real 
quantities g*(x) and the vector part by y,*(x), 
where the subscript ¢ denotes the components in 
the ordinary space and takes the values 1, 2, 3. 
The superscript @ in both cases denotes the 
components in the isotopic spin space. In the 
symmetrical theory it takes the three values 
1, 2, 3; in the charged theory the values 1, 2; 
and in the neutral theory only the value 3. In 
the following we do not restrict the values which 
a may assume so that the results for these three 
types of theories may be obtained by making 
suitable restrictions on the values for a. How- 
ever, unless otherwise stated, we restrict our 
attention to the most interesting case of the 
symmetrical theory. The Hamiltonian of the 
field and its interaction with a nucleon is given by 


H=3E f {(x9)?+-(9 o8)2+02(g8)3}dV 


+4 [leet (ovetteveriav 











dg" 
said 4 ammenities 
we » J Iz, 2 
g(a aye 
a —_ Ci; re V, 
(~ ~~) | ee 


where x* and w;* are the momenta conjugate to 
g* and y,*, respectively; «x and yw are the rest 
masses of the pseudoscalar and the vector 
meson, respectively (we are using natural units 
where A=c=1); o; and r* are the spin and 
isotopic spin matrices, respectively, with o;; 
= —oj;=0, if 1, j, k are cyclic permutations of 
1, 2, 3; and f and g are coupling constants for 
the pseudoscalar and vector field, respectively.® 
For the time being we do not make the assump- 
tion f=g. U(x) is the source function of the 


5In this paper we find it more convenient to use these 
constants which have the dimension of «~!. Thus our f is 
equal to g/\2« where g is the coupling constant for the 
pseudoscalar field used in P-D. 
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nucleon which we assume to be spherically 
symmetrical and normalized according to 


fu@av=1. (2) 
It determines a radius a of the nucleon 


-= ff hd aye (3) 


where R= |x—x’|. 
As in P-D, we split the fields into two parts 
and write 











on _ f oo" f S, V 
saa (4) te. gy be, ’ 
Oa Oy; dyp* 
ij = (2)! _ U(x)dV 
v (2m) f Ox; Xj | @) (4) 


U dU 
= (2) f vee ——|d £ 


Ox ; Ox, 


0U 
[ e—av=o, 
Ox; 


0a 


Wi; is antisymmetric in 7 and j, and we shall 


; ‘ oa 0 — ; 
sometimes write ¥;; =v. for 1, 7, k cyclic. The 
source function generates the potentials X (x) and 
Y(x) according to 


(—A+nx)X=4rU, (—A+yp?)¥=4rU; (6) 
so that 
X(2)= f Ue ye" ‘R)dV’, 


Y(x)= f U(x')(e*®/R)dV', 


with R= |x—x’|. We define 


0X aU 
16;;= ——dV, (8) 
Ox; OX; 


and if we assume both xa and wa<1, we have 


0X dU dY au 
—_—dV= | ——¢dV. (9) 
Ox; OX; Ox; OX; 
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This approximation gives some formal, although 
unessential, simplification and we shall make use 
of it. We therefore write 


E(x) =X(x)/I, n(x) = ¥(x)/T; (10) 


then 
d& AU On OU 
fs ze- — -—=§,;. (11) 
Ox; OX; Ox; OX; 


We thus obtain from (4) and (5) 








re) 
"(x)= (44)? 2D oo a + 92), 
(12) 
¥i2(x) = = (8a)- ; ~ vis om ty, ‘a(x D 


The corresponding decompositions of the 
momenta are given by 


; at 
r= (4) [ ea V, 
Ox; 


ta On 0 
aes a (8n)-*f fo ws : lav, 
Ox; Ox; 


- 0a 0a PP on . 
with w;; = —w 5 =w, (4, j, k cyclic); 


dg 
fr dV=0, 
Ox; 


(13) 








(14) 





aU aU 
Nouy=a4n f —~wr: (15) 





a(x) = (42)! > mr P* 


0a dU(x) 
w*(x) = (2m)! ~ Wij 





+w;’*(x). (16) 


j 


The commutation relations are 
iL ri, oj? |= dbap6ii, 


a. 0g 
tL wis, Wer) = bap ( bind j2— 5:25 jx). 
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The Hamiltonian (1) becomes 
N 
— ss { (32)? + (w,°*)?} 


1 
tr - { (g.°*)? + (W.%*)?} 


+2 { feP@+v2gy°*} oi7* 


+(40)) + a “av 


a,i,j 

wij dU 
wav 

v2 Ox; 


F4E fleet (oav 


+42 [ to'(-ate0" 


+y;'*(—A+yp*)y,'*}dV, (18) 


and the potential energy in the zero state E° is 
given by 


1 
E =| M6. )?+ (w.iP*)?} 


+ [fete+vieri}o.r*| (19) 


The total angular momentum of the field and 
the nucleon is 


Oy* dg* 
hige~ f| “( - ) 
. 2 , Ox; "Ox; 

OW, OW,.* 
me {ex( ~ _? Y )Jav 
a,k Ox; Ox; 


+2 [ Wewr—vroaldV+houi (20) 








dV 








and the part L;;° due to the zero state of the 
field is 


Oa Oa 0a 


L;i= 2 (giten;" = 9 ;*r®*) + D. (Winww jx me Vinwin r), 
a,k 

or 

Li?= Le { 9° i = ene 


tp Ptw 9 — P %w,*} . (21) 
Now E° takes its minimum value when 


gi=fle*, ye*=v2gle;*; (22) 
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so that it is convenient to introduce new vari- 
ables defined by the orthogonal transformation 


B88 (fP+2¢)-H{ for tvIgde}, 


(23) 
v= (f?+2g7)- t{v2g9,°" — fir}; 
with the inverse transformation 
oe= +2 2 —if } «+ v2 W ; 
got = (f?+2g7)-4{ f evo} (24) 


Voa= (f+ 2g!) Vgboa—fU.} 


The corresponding transformation equations for 
the momenta are 


11,9" = (f2-+ 2g") frs9*-+ V2ge0*} 
0,9¢ = (f?+ 2g?) 4 { V2gr 9 — foi} ; _ 

with the inverse 

mt = (f?+-2g)-4{ 11°" + v2g0,"%}, 

wi = (f?+ 2g?) { V2gIT,°* — f2,°%}. ei 

Then we have | 
pm Lis= List Lis (27) 
Lis= Xa (PM oe — 411%}, (28) 
Lij= Cal Vor Pn— ¥ P20}. (29) 

As in P-D, we write 

$,* = De*+ dis qveF, (30) 


where g*=q**, e;* are the components of an 
orthogonal matrix introduced in Eqs. (51), (51a) 
of P-D and can be considered as components of 
an orthogonal system of three unit vectors, and 
according to (22) 


ai = (P+ 2¢¢)MT; (31) 
IT,°* = Dia{ pe +(1/2D)L* je, (32) 

where p= p** and 

[8 => e; ae fL.,= —D dL lerllY— Pll}. (33) 


i,j 


In (32) and (33) terms of higher order in g%*/D 
have been neglected [cf. P-D, Eq. (72)]. p® 
and g® satisfy the commutation relation 


ip, g” ]=}(da75¢s+ 5as5sy), 


[cf. P-D, Eq. (57)]. If we insert these new 
variables in (18) and perform the S-transfor- 


(34) 
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mation, we obtain for the Hamiltonian 


N 
H=—s?— }(f?+2g*)I 
4D (f?+ 2g”) 








nif aU , 
+——___ | _ Late f —r'ad | 
D(f? +2g7)) 43: Ox; 
7. - pms { [VUXet ued V 
D(f?-+2g2)) pi a 
HE Jleereoraay 
ab 35 
+E )*, (35) 


where - esis 
=3 DL (L%*)?. 
a, B 

The terms which do not influence the isobar 
separation have been omitted; in particular the 
terms due to W,°* and 2;°* which describe free 
mesons. By shifting the origin of ’* and w,’* 
according to 


nif 











g'*@=_"e~— ee L8e Pad 
D(f?-+2g2)4 5 
Ie 30 /Ox; 
Ox; Vi)'dV 
xi S(V8&) (36) 
mig 
w,;/*=w,'!* Las 
D(f?+2g*)''s 
x| ex (¥ =e )| 36) 
Sivo-av) |; * 
we can get rid of the terms linear in 2’* and 


w;* and of the order 1/D, and still retain the 
orthogonality relations (14) [cf. P-D, Eq. (75) }. 
The final result for the isobar energy is 


3 s? An 
4 (f?+2g*) S(VX)*dV 
and the isobar separation AE is given by 


3 {s(s+1)—#} 4 








““. (Pte Sx ©? 
Now for a small source (xa and wa both<1) 
4 1 
T(wX)dVa 


* We have s instead of L as in P-D since we want to 
use the usual spectroscopic notation in the discussion of 
the two-nucleon system in Section IV. 











404 Ww. 


[cf. P-D, Eq. (80a) ], hence 
3a 

a 

4(f?+ 2g”) 

We obtain the result (80a) in P-D for pseudo- 

scalar meson alone if we put g=0 and f=g/v2k. 

For the vector theory alone, we put f=0 and 

find that in this case AE is 1/2 of the value in 

the pseudoscalar theory (for the same coupling 

constant). For the Mdller-Rosenfeld mixture, 

we put f=g and find AE in this case is 1/3 of 
the value in the pseudoscalar theory. 


{s(s+1)— Z}. (38) 


III. NUCLEAR FORCES 
A. Pseudoscalar Theory 


We make the simplifying assumptions that 
the source size a is small compared to both x“! 
and the distance r= |x;—xn| between the two 
nucleons. Let Ui(x) and U(x) be the source 
functions for particles I and II, respectively, 
and normalize them according to 


facoav=1 for A=I, II. (39) 


The potentials X1(x) and Xn(x) which they 
generate according to 


(—A+x«?)X4=4rU, for A=I, II, (40) 


are given by 
eR 


Xa(e)= f Uae’ R 





dV’, (41) 


where R= |x—x’|. Also let 








0X4 9U 4 
I6;;= —dVforA=I,II; (42) 
Ox; Ox; 
0X4 0UB 
J i= ——dV for A, B=I, Iland A+B 
Ox; Ox; 


(43) 
oe? ek 
=— ff cence ‘ae R )avav’ 


It is evident that J;;=J;;. Now J is of the order 
a~*, and J;; is for small r of the order r~* so that 
in the following we can neglect all quantities of 
the second and higher orders in J;;/J.7 More- 
over we can substitute for J;; its value for point 








7A similar situation holds for the integral /(0X1/dx;) 
X (@Xy1/dx;)d V in comparison to §(0X1/dx;)(0X1/dx:)dV 
which are of the order of magnitude r~ and a™, respec- 
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sources ;° i.e., 


J s (—) (44) 
a Ox1,0x1;\ 7 


where r= |X;—Xr1}. 








To obtain the interaction energy between two 
nucleons, it is sufficient to consider just the 
potential energy of the system E which is 


E=> fe(-atereav 


re) 
~(4e)'f £ 5 eairatUa (x)dV, (45) 


A, a, i Xi 


where the subscript A is summed over I and II. 
As in S-D, Eq. (71), we split the field as follows: 


0a 
g*(x) = (49)! DS paitai(x)+e'(x), (46) 
A, é 
where ¢’? satisfies the condition 


OU 4 

fea V=0 for A=I, II, (47) 
Ox; 

and where the functions £4; span the linear 

subspace of the functions 0X,4/dx; in such a 

way that 





dUz 
f ba av = bani (48) 
Ox; 
From these relations we obtain 
= 0U 
eai= (4) fo “dV. (49) 
Ox; 


The above conditions are fulfilled up to terms 
of the order J;;/J if we put 
1 OX, 1 OX p 
co —— ij 
Tox; I? ; Ox; 
for A, B=I, Il and A#B. (50) 


In terms of these variables, the potential energy 
(45) becomes 





0a 
+f 7 GA ‘TaiTa® 


A, a, i 


1 0a. 2 Jj Oa 0a 
E= —(pai) ——en.¢1 
Za —_— en 


+z fom(-atererar. (51) 


tively. It is shown in S-D that, due to this circumstance, 
the neglect of the dependence of the isobar levels on r is 
justified as long as the latter integral is large compared to 
the former. 

§ The substitution of the point source is in reality not so 
simple and holds only for distances r considerably larger 
than a. This was pointed out by Oppenheimer and Serber. 











a 


W 


I 











The justification for the assumed decomposition 
of ¢* lies in the circumstance that no cross terms 


between ns and yg’ occur in (51). g’* describes 
free mesons and can be omitted for our purpose. 
For the case of infinite separation of the two 
nucleons (J;;=0), the lowest eigenvalue of the 
interaction energy is obtained analogously to 
the treatment of the one-source problem in P-D 
(cf. Eqs. (82), (64), (65), (67) there]. We put 


vac= flea: (52) 


a . . . 
where the e4; satisfy the conditions 


aa a Bg 
dia Calaj=5i;, Doi CaCai=Sap, 
for A=I, II; (53) 


and by selecting the state where 


¥ caitateac=—3forA=I,II; (54) 


a,i 


we obtain for A=I, Il 


o4'*t4*= —1 for each value of a (55) 

with 
oa*=) 5 CATA: (56) 

and 
o4'*748+04%7r4*=0 for a XB. (57) 


To obtain the minimum value E° of E to the 
first order in J;;/J in the absence of free mesons, 
we simply insert (52), (53), (54) in (51). The 
result is® 


P= —3f2l—f?d Jijeneny, (58) 


a,i,j 


or omitting the first term which gives the 
self-energy, we have for the interaction energy 


Rue ~*T Sede (59) 


a,t,j 
Inserting for J;; the expression (44), we obtain 
0 , 
Ei 1n=f? Dia (€1*° V) (enr* V)(e*"/r), +—(60) 
and if the differentiation is carried out, 


, f2{_«? 1 
Ey w»=—} l--e-*"" +A—(3 + 3ar+ x*r?)e—*" (61) 
' 3l or - 


* If we had taken a more general expression 


0a a 
¢Ai = Lal flbap+qa™)eai 


instead of (52), we would have obtained in the expression 
for E® terms in ga which are of the relative order J;;/J, 
and hence give corrections of the order (J;;/J)? to E°. 
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where 
T=) >. (e1*- er"), 


3 
A= =| (er -X) (€11* -X) — (1% - @r1*) (62) 


with x=x;—Xr1. 

It follows from the method of derivation that 
our theory differs from the usual perturbation 
treatment in the weak coupling theory only in 
the fact that e474* is replaced by e4* (A =I, II). 

As stated at the beginning of Section II, the 
symmetrical theory is obtained by letting a take 
the values 1, 2, 3; the charged theory for a taking 
the values 1, 2; and the neutral theory for a@ 
restricted to the single value 3. 


B. Mixed Theory 


As in Part A we consider two source functions 
normalized according to (39). They generate 
the potentials according to 


(—A+)X4=40U 4 


(—A+u) ¥.=4eUs for A=I,II; (63) 


and hence 


e7*k 
X 4(x) = [ Uae a Vy’ 
for A=I, II. (64) 











"on 
Ya(x -fu x’ dV’ 
a(x) A( > 
We have 
0X4 dU dY, dU 
Tou= f *—“ave[ a 
Ox; Ox; Ox; Ox; 


for A=I,II; (65) 

0X, dU oe mer 

Jninn f 2 8eaye_Y (X) 
Ox; Ox; OX} ,0X1; r 

for A, B=I,I1andA#¥B. (66) 

eects eS. 


Ox; OX; 0x1;,0X%1;\ 7 








K,;=K,= 








The potential energy of the system is given by 


E=> f ye(—A+e2) 9*-+y(—A+y?)ye}dV 





dg" 
— (4dr)! CAi 
(4n)'f X {= 


A, a, i, j Xi 








1 ee ayia 


? - Joa} ratUad V, (67) 


Ox; OX; 














406 WwW. 


where we have assumed the equality of the 
coupling constants. 
As before, we split the field as follows: 


(x) = (40) gaitai(x) +¢'"(x), 
A,i 


(68) 
Yit(x) = (8x) D Paisnas(x) +¥i"(x) ; 
A, j 


with the orthogonality conditions 


0U. 
fe ad V=0 


Ox; 





for A=I, II; (69) 
0U4 
for avac 
Ox; 


dUz 
fe dV =6486;;, 


Ox; 


dUz 
fms dV =6,436;;. 


Xj 


and 





(70) 





Then 


a 0U. 
we f o—aV 


xi 





for A=I,II, (71) 
0a U4 
VAij= (8x)! f y—“av 
Ox ; 


where we do not assume yet the antisymmetry 


Oa , ° ° — . 
of ¥4:; in 7 and j. The above conditions are satis- 
fied up to term linear in J;;/ZJ and K,;/I by 


10X,4 1 OXp 





Ome, Fe, 
for A, B=I, IIland A+B. (72) 
10Y,4 1 OY, 
NAi ij 


“? Ox; 7} j OX ; 


The potential energy now becomes 


1 Oa. 2 da .2 
E=— DL {(¢ai) +3(Wais) } 


2I A,a,i,j 
1 Oa 0a 1 0a Oa ’ 
-— D> {Jisgrigu yt 3 Kirundi} 
I a,i, j,k 


a 1 a 
a | casratonrt oastavau} 
+2 f te'(-at ee" 

| +yi'(—A+u%)yi'"}. (73) 
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To obtain the minimum value of E to the 
first order in J;;/J and K;,;/TJ, it is sufficient to put 


eai= flea: 
for A=I, II andi, j,k cyclic, (74) 


Pau™ — Wai =Var= Vi fleas 
4 E hescratacteactatas) = bag 
for A=I, II. (75) 


We find the interaction energy to be 
0 aa aa 
Ein=—-f?Dd {((Jis—Kipenenj;—Kiserjen;}, (76) 
a,i,j 


and inserting the expressions (66) for Ji;, Ki; 
for point source, we obtain 





0 e-*'—e-#r 
Ein=f? E{ (er-¥)(eut-¥)(* ) 


r 


+ (err-en2)a(“) ; (77) 


r 





The evaluation of the differential operators 
yields 


Et n= }{TI(r)+AK(r)} (78) 


with T and A as defined in (62) and with the 
radial functions 
f? 
J (r) =—(x*e-*" + 2y?e"), 
r 
(79) 


2 


K(r) wks (3+3xr+x*r?)e-*" 
r 
—(3+3ur+p*r?)e"}. 


It is to be noted that for small 7, both J and K 
behave as r~', and K(r) =0 for x=y in agreement 
with the results of Mgller and Rosenfeld for the 
weak coupling theory. 

The remarks made at the end of Part A also 
hold in this case; i.e., by letting a take the values 
1, 2,3; 1, 2; or 3, we obtain the symmetrical, the 
charged, and the neutral theory, respectively. 


IV. THE DEUTERON 
A. Classification of the States 


The stationary states of the one-nuclear sys- 
tem in the symmetrical theory can be described 
completely by the numbers s, m, and nm which 








al 


al 


eis 
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are the eigenvalues of the following operators :!° 
s?=t?=s(s+1), ss=m, =n. (80) 


s is a positive half-odd integer, and m and n are 
»half-odd integers such that —sSm,nsSs. It is 
‘to be noted that there is complete symmetry 
between spin and isotopic spin. 

For a system containing two nucleons, the 
Hamiltonian is 


1 L? a : 
' n=—|p++—} +— {s1?+s1?— §} 
M 


4f? 
+43{TJ(r)+AK(r)}. (81) 


The first part gives the kinetic energy of relative 
motion," the second the isobar energy, and the 
third the interaction energy. The motion of the 
center of gravity of the system is not considered, 


“ and of course this is in the non-relativistic 


approximation. 

Let us first consider states with no orbital 
angular momentum (L=0); then there is still 
complete symmetry between spin and isotopic 
spin as in the one-nucleon system. Their opera- 
tors S=s;+sy and T=t;+ty are still constant, 
but now their magnitude need not be the same. 
Hence the stationary states can be described by 
the eigenvalues S, 7, M, and N of the following 
operators: 


S?*=S(S+1), 

T?=7(T+1), 

S;=M, (82) 
T?=N, 


where S and 7 are positive integers or zero such 
that |s;—sn| =S, T=(s1+5n), and M and N are 
integers or zero such that —S=M=S and 
—T=NST. If there is an orbital angular 
momentum L, then the spin angular momentum 
S is no longer constant, and it has to be replaced 
by the total angular momentum J given by 


J=L+S, (83) 
and the eigenvalues S and M by J and Jy, 


1° Cf, P-D. We use s instead of j there to denote the 
eigenvalue of the spin. Also we denote the spin operator 
by s instead of L as‘mentioned already in footnote 6, and 
the isotopic spin operator by t instead of T. 

1! We identify M with the empirical mass of the nucleon 
from which the part of the mass due to the meson field has 
been subtracted. It was pointed out to one of us by Dr. 
Schwinger that in the neutral theories the latter is a tensor 
with respect to the direction of e*. In the symmetrical 
theories, however, the field mass is also a scalar. 
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FIELD 


respectively, where 


PH=J(J+1), Ju=Lut+M. (84) 


In addition there is the quantum number pro- 
vided by the parity, the eigenvalue of the spatial 
reflection operator. The Hamiltonian (81) is also 
invariant under any rotation in the isotopic spin 
space. Such a rotation leaves the s4; unchanged 


and transforms the ¢4* and ea; (for each 7) in the 
same way. In particular, a rotation which 
changes the sign of the 3-component, as for 
instance the rotation about the 1-axis through 
180° which results in the transformation 


tts, ty2—-ty’, th®—t,', 


(85) 


e,’—e,4!, e42—>—e4’, e4*——e,', 


will turn out to be very useful. [Note that this 
transformation conserves the commutation rela- 
tions satisfied by the e,* given in (135) of the 
Appendix. ] The transformation (85) multiplies 
the spin function by a phase factor e*°°™**, and 
changes the sign of m, myx, and N. Thus the 
operator shows the degeneracy in the energy 
levels of the excited spin states which differ only 
by the sign of N. In particular for the case 
N= 0, that is for the deuteron, the phase factor 
is unity, and the transformation does not change 
the state to which the spin function belongs. 
Consequently this operator provides another 
quantum number which is the signature of the 
spin function under the substitution m;j——,, 
nyi— — ny. 

With the aid of these quantum numbers the 
complete wave function of a stationary state 
can be written in the form 

ur, s(7) 





V(r, v, gi J, T, Ju, N)= > cf . 
Ls oan r 


LM+M=Jy 
L) 
x Vz (8, g)z(S,T, M,N), (86) 


where c¢ is the normalized coefficient of the 


J 
LMM 
Clebsch-Gordon series; uz,s(r) is the radial 
function which we normalize, for bound states 
where this is possible, according to 


2 


> tz, s(r)ux, s(r)dr=1; (87) 
L,8 


0 


Y,/“(8, ¢) is the normalized spherical harmonic; 








- 
: 
| 
: 
: 
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and =(S, T, M, N) is the spin-isotopic spin 
function (we shall refer to this hereafter simply 
as the spin function when no possibility of 
ambiguity exists). Of course the numbers L and 
S do not determine the radial functions uniquely 
since there is still the radial quantum number, 
but in the present work we consider only the 
radial functions belonging to the lowest energy 
level. In accordance with the exclusion principle, 
we assume W to be antisymmetric in the inter- 
change of particle I and II. Then since the 
interchange of the space coordinates gives rise 
to a factor (—1)¥”, the spin function = must be 
symmetric or antisymmetric in the interchange 
of particles I and II depending on whether L is 
odd or even, respectively. Further since the 
parity quantum number is (—1)¥, any state is 
restricted to spherical harmonics of even or of 
odd orders. 

In order to evaluate the eigenvalues of the 
operators [ and A, it is necessary to evaluate 
the matrix elements of e;* in the (s, m, nm) 
representation and express the spin function 
=(S, T, M, N) in terms of the one-nucleon spin 
functions &(s1, m1, my) and én(su, mu, mu). The 
calculation of the matrix elements of e,* is 
given in the Appendix. The expression for 
=(S, T, M, N) in terms of &(s1, m1, mr) and 
tn(su, mu, mu) is obtained by means of the 
Clebsch-Gordon series. In fact they can be 
expressed as products of the spin and isotopic 
spin functions 


=(S, T, M, N)=X(S, M)Z(T, N), (88) 
£4(Sa, Ma, Ma) =xaA(SA, Ma) Sa(Sa, Ma) 

for A=I, II; (89) 
and 


x(S, M)= i cmrmrt {x1(St, my) x11(St1, M11) 


ST SII 
mi +my1 =M 


+x11(S1, m1) x1(Si1, M11) }, 
(90) 


a77TN= 2 carers { $1(St, m1) f11(Si1, M11) 
SI, SII 
ni +niI=N 


+f11(S1, m1) £1(Si1, 11) }. 


The symmetrization of these expressions in 
particles I and II is necessary for the reason 
given above. However since I and A have matrix 
elements between states with values of s; and su 


which differ by 1, s: and sy are not constants of 
the motion, and hence the series (90) are actually 
infinite series, and the calculation of the eigen- 
values of I and A are very complicated. 

We therefore make the following approxima- 
tion : we replace I and A by Ip and Ao which are 
the parts of [ and A diagonal in sy and sy. 
This treatment is of course rigorous in the limit 
of infinitely large isobar energy. We shall show 
for the ground state of the deuteron, by treating 
the off-diagonal parts [,; and A; of I and A as 
perturbations, that even when the isobar energy 
is not too large, the effect of the perturbing 
terms is small. This justifies our method of 
approximation at least of the ground state of 
the deuteron in which we are interested. 

In the unperturbed system, s; and sy are also 
quantum numbers of the system, and hence, 
combining (88), (89), and (90), we obtain for the 
spin functions of the unperturbed system 


S T 
=(sr, SII; S, T, M, N)= Om Cmim[Cnyntit 
mi +mi[=M 
ni +ni,[=N 


X (E1(S1, mr, m1) E11 (Str, Mr, M11) 
+ £11 (S1, mr, M1) E1(St1, Mir, M1) }- (91) 


It is readily seen that these spin functions are 
eigenfunctions of the operator given in (85) with 
the eigenvalues (—1)*!+*"~7, It is only necessary 
to use the following property of the Clebsch- 
Gordon coefficients : 


st+st1-T T 


pe oe Suteien, (92) 


and to remember that the operation (85) on 
the spin functions (91) for the deuteron simply 
replaces m; and ny by —m; and — mn, respectively. 

The spin functions (91) are also eigenfunctions 
of the operator Io, and its eigenvalues are found 
by simple calculation with the matrix elements 
of the e,*. They are given by the general formula 

| S(S+1)—s1(s1 +1) —Ssu(su+1)} 
—_ x {7(T+1) —si(s1+1) su(sut hj (93) 
451(s1+1)su1(sir+1) 


Thus if s is the smaller and ¢ the larger of the 
two numbers s; and sy, we have 


s s(t+1) 
—-—_ Sh= ; 
s+1 (s+1)t 





IIA 


(94) 

















— we . 











TABLE I. Classification of the deuteron states. 











Sy Sty Ss - (LID (n-—n) To 
1/2 1/2 0 Oo + ~ +1 

1/2 1/2 0 1 - + —1/3 
1/2 1/2 1 O - = —1/3 
1/2 1/2 1 1 + + +1/9 
3/2 1/2 1 1 . < oes - +5/9 
3/2 1/2 1 2 . i + —1/3 
3/2 1/2 2 1 +,- - —1/3 
3/2 1/2 2 2 +,- + +1/5 
3/2 3/2 O O + al +1 

3/2 3/2 O 1 = + +11/15 
3/2 3/2 0 2, + - +2/5 
3/2 3/2 0 3 - + —3/5 
3/2 3/2 1 0 - = +11/15 
3/2 3/2 1 1 + + + 121/225 
3/2 3/2 1 2 - - +11/75 
3/2 3/2 1 3 + + — 11/25 
3/2 3/2 2 0 + - +2/5 
3/2 3/2 2 1 - + +11/75 
3/2 3/2 2 2 + - +1/25 
3/2 3/2 2 3 - + —3/25 
3/2 3/2 3 0 - - —3/5 
3/2 3/2 3 1 + + — 11/25 
3/2 3/2 3 2 - - —3/25 

3 3 + + 





3/2 3/2 +9/25 





so that for large s; and sy, Ip lies between —1 
and 1. Furthermore, Io is negative when either 
S or T is near the maximum value s;+5s,; and 
the other near the minimum value |s;—s|, 
while Ip is positive when both S and T are either 
near their maximum or near their minimum 
value. This result is very useful since 4I9J(r) is 
the dominant part of the interaction energy, and, 
in general, the sign of Ip determines whether the 
interaction is attractive or repulsive. Since the 
operator Ao contains angular coordinates besides 


the operators eax, the calculation of its eigen- 
values is too complicated in the general case, 
and they are computed only for the special 
states considered in Part B. 

The first few states of the deuteron system 
are classified according to their quantum num- 
bers in Table I. The fifth column under (I, II) 
indicates the symmetry of the spin function 
under the interchange of particles I and II. 
For s:#su, both symmetric and antisymmetric 
functions are possible for each value of S and 7; 
but for s:=sy, only a symmetric or an anti- 
symmetric function is possible for a given value 
of S and 7, depending on whether (S+7) is 
even or odd, respectively. The sixth column 
under (n—>—n) indicates the signature of the 
spin functions under the operation (85). As 
mentioned above, it is positive or negative 
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depending on whether (s1+sn—T7) is even or 
odd. The last column under Ip gives the eigen- 
values of this operator. 


B. The Ground State 


The lowest state of the unperturbed system is 
for s:=sn=4, and for the spatial part of the 
wave function to consist mostly of the S-wave 
(L=0), we must consider antisymmetric spin 
functions (in the interchange of particles). We 
see from Table I that both the singlet (S=0) 
and the triplet (S=1) states have I'p=—}. 
Thus these states will evidently be bound states. 
Now the tensor interaction given by Ao vanishes 
for the singlet state, and we shall find below 
that it contributes additional attractive potential 
to the triplet state. Thus in agreement with 
experiment the ground state is the triplet state. 

To obtain the matrix elements of Ao, we use 
the matrix elements of the e;* and the expansion 
(91) for = in terms of & and &,. The result is” 


4(2x)! 





AoZo(1, 0, M) = — 


x YL cwem¥2 Zo(1,0,M"). (95) 


M'+M''=M 
Now if we use the spectroscopic notation to 
designate the states; i.e. write *Sy' for the 
triplet S state with J=1 and Jy=M, *Dy' for 
the triplet D state etc., and if we let Zo(*Suas'), 
2o(*Dy'), etc. denote the spin and angular parts 
of the wave function, we have 


Lo(®Sa') = (44)—4Z0(1, 0, M), 


‘tue (96) 
Yo(?@Dy') = is CmM’M"* Y2 Zo(1, 0, M ). 
M'+M’''=M 
Thus combining (95) and (96) we obtain 
2v2 
AoXo(?® Sx!) = 2D’), 
(97) 


2v3 
AoXo(@Dy') => #> @Dy') = sy 20 Su!) 


No higher values of the angular momentum 


2 Instead of writing the spin values sy and sy as in 
(91), we now use the subscripts 0, 1, and 2 to denote 
quantities belonging to the states sy = syy = 4; s1= 4, sur =4; 
and sj =s11 = 4, respectively. Also we drop the designation 
0 for N since we do not consider states where N has values 
different from 0. 











> ce RE 
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occur in this state, and we have 


uo(r) 


Vo(r, 3, 9; 1,0, M)= Do(®Sar') 





Vo 





Os.¢Dw). (98) 


r 


oa 


Inserting this expression into the Schrédinger 
equation 


HV=EW, (99) 
and using the expression (81) for H, we obtain 


with the aid of (97) the simultaneous set of 
differential equations for uo(r) and v9(r) 


+z +=J( MK 
a? 0 9 r) ug >= — 9 (r)v0, 
(100) 
dv 6v0 1 2 
| Ext oJ) —K) bn 


2v2 
= ——p MEG ae. 


Thus we see that the tensor force contributes 
additional attractive potential. The above equa- 
tions are to be solved with the normalization 
condition 


f (uo?+v9?)dr=1, (101) 
0 


and the boundary conditions uo, v%=0 at both 
r=0 and r=. The solution of these equations 
will be studied in Part C. 

First we want to see whether (98) is a good 
approximation to the actual ground state wave 
function ; i.e., to see whether the approximation 
we have made in taking only the parts of I and 
A diagonal in s; and sy is justified. We do this 
by calculating the matrix elements of the pertur- 
bations I; and A; between the state described 
by (98) and the other states. The first excited 
spin state occurs when one of the particles has 
spin 3/2 while the other has the unexcited 
value 1/2. We see from Table I, however, that 
none of the states for this case has the same 
eigenvalues as any of the states with sy=sy =}, 
and since neither I’, nor A, change these quantum 
numbers (except s; and sy), both I’; and A; have 
no matrix element between these states. 


For the second excited spin state, s;=5y=3, 
the states which combine with the ground state 
(S=1, T=0) are only the triplet and the septet 
terms with T=0 (cf. Table I). The spin-angular 
functions which occur in this state are 


LY2(*Sar') = (4) $Z2(1, 0, M), 
1 M’ 
Y2(?Dxy') = > CarmeYe 22(1,0, M”), 
M’+M''=M 
(102) 
.M’ 
ZCDe= YS caeme Ys Z2(3, 0, M"), 
M’+M?'=M 
1 M’ 
27Gu')= DL cum Vy 22(3,0, M”). 


M’+M’''’=M 


The operator Ao applied to these functions gives 
the following results: 



































40/7 344/2 
AoZ2(2.Saz') = 35 Y2(7Dy') — 75 Y2(*Dy'), 
12/6 4/2 
Ao®2(?Dy') = Y2(7Gyz') — v Y2(7D yy") 
25/7 25/7 
4, (6D aes 4 (Sy!) 
757 M 75 Le(*Saz'), 
Aot2(?7Dy') = — 12V3 5 Gu) + I SA (D u') 
or 175 175 
(103) 
4/2 4/7 
- D.(8Dy! >(3Sy'), 
25/7 2(*Dy')+ = 2(®Saz') 
AoX2(7Gay') = OS A(°G uw’) _2V3, 0D uw) 
: 7 — 
12./6 
+ >2(*Dy'). 
254/7 aintel 


Writing the wave function in the form 











u2(r) 
V2(r, v, Y; Bi 0, M) = Y2(*S4') 
4) dy) Os, cDy) 
r r 
4225 0G), (104) 





we find the following set of equations for the 








ri 








te 
et 
ar 
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radial functions: 


aus vy ule 3a 4] 
dr? “oe aT" 


2 
= Fg MK 17V2 v2— 6/7 wy}, 





dv. 62 3a 34 
+M{ E.-——s¢ )—-—K(r) lr 














ir 2f? 15 225 
2 6/2 18/6 
= ——_MK(r)1174/2 i 
325 | V2 ust 7 We V7 ya}, 
d*We 6w2 3a 
—_— M|B.- = 
dr? sr? 2f? 
(105) 
+20) (0) |w 
175 : 
4 V2 3/3 
_ }- r nas eh 
ns 1 “4M Es os + K 
oe = Jin) -= ly. 
4 V6 V3 | 
—! a 
25 |- a an 


No higher spin states need be considered since 
neither I; nor A; has matrix elements between 
states whose spin value differ by more than 1. 
Hence the excited state considered above is the 
only state which combines with the ground 
state (98). 

To calculate the matrix elements between 
them, we have 


V5 


T20@Sa') = 3 22(?Saz'), 





V5 
P\2o(@Dyy') = a 2a’), 





v2. 3 1 
35" Du +106) 


3/6 2 
gg out 35 


AiXo(@Sy') = — V1S 0D u')+ 
; V5 : 








AiZo(@Da7') = >2(7Dy') 


>. (3.1 v2 
Wy alae ae 


Hence if we write (2!H|0) for the matrix 








22(8Syz'). 


element of the perturbation 
A, =3{TiJ(r) +Ai1K(r)} (107) 


between these states, we have 


(2| H|0) =~ SS [tant oan 700 








4 | Famts 3/5 -Bo(4/2 Uo — Vo) 
1 , 
_ "ani Vo) 
~=¥ 90] K(0 fir (108) 


From (105) we see that for the dominant S 
part “2 of the wave function the direct interaction 
gives rise to a repulsive potential 11/15/J(r). 
The tensor force gives rise to an attractive 
potential when coupled with the triplet D part 
ve, and a repulsive potential when coupled with 
the septet D part we, but these contributions 
will be small compared to the direct interaction. 
Hence it seems that the state (104) will not be 
a bound state. Assuming this to be the case, 
we can estimate the correction to the wave 
function (98) due to the perturbation (107) by 
approximating the wave function (104) with a 
plane wave. If we denote this correction term 
by Vo, then we find 


fivaitav~o2, (109) 


by using numerical values for the constants a, f, 
«x and yu which will be found in Part C and for 
isobar separation of 20 Mev. Thus the correction 
is quite small, and it will not alter appreciably 
the results obtained for the unperturbed ground 
state. 

On the other hand, if the state (105) is a 
bound state, the effect of the perturbation is 
proportional to the ratio 


(2|H|0)/(E2—Eo)=(2| H|0)/Eisov (110) 


where Eisop = 3a/2f?. We would certainly obtain 
an upper bound for (2|H|0) if we insert in 
(108) 20, vo for #2, v2 and neglect the other small 
quantities. This gives for (2|H|0) a value a few 
times the binding energy of the deuteron, 
Evina= | Eo|. Hence in this case the approxima- 
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tion is justified since the upper limit of the 
ratio (110) is about 0.3, and actually it will be 
much smaller. 

In either case, the validity of the approxima- 
tion depends on the smallness of the quantity 


Epina/Eisov- (111) 


C. Determination of the Constants 


In the present theory the following constants 
occur: the size of the nucleon a, the strength of 
the coupling f, and the masses of the pseudo- 
scalar and the vector mesons, «x and yu. Of these, 
f must be chosen so that the binding energy of 
the deuteron is equal to the observed value of 
2.17 Mev, and «x equal to the observed value of 
the mass of the meson of ~200 electron masses 
(we take for convenience x=177 m). The tensor 
force depends on the ratio of the masses u/x 
and vanishes for 1=«x. Thus we shall choose u 
so that the quadrupole moment of the deuteron 
comes out equal to the observed value of 
2.73 X 10-7 cm*. However we have the following 
conditions to satisfy for the validity of our 
calculation : 


(1) small source, i.e., ax<1; 

(2) strong coupling, i.e., (fx)*>(ax)?; 

(3) small effect of the higher spin states on the 
ground state, i.e., Epina/Eison1. 


It is not possible to say a priori whether these 
conditions can be satisfied with the above 
determination of the constants, but as we shall 
see below, this is possible. 

In order to determine f and yw, we must 
consider the differential equations (100). f is to 
be determined so that the lowest eigenvalue E» 
of these equations is equal to — pina, and u so 
that the quadrupole moment" 


— ( : ‘)a (112) 
e-— the uv wi r 


agrees with the experimental value. Since the 
Eqs. (100) cannot be solved in terms of any 
known function and since the conditions deter- 
mining f and uw are interdependent, we have to 
assume some reasonable value for yu, solve (100) 
approximately, and see whether it gives a good 


18 Since from now on we consider only the ground state 
wave function (98), we drop the subscript 0. 


value for Q. Let us first simplify the equations 
by introducing the dimensionless variable 


x=2(M|E,|)'*r (113) 
which transforms (101) to 


d*u/dx?— {4—j(x)}u= —v2k(x)v, 


d*y 6v : 
—— =~ 14-5 (x) +h(x) jv = — v2k(x)u, 
dx? x? 
with 
; . 
(x) =—{a*e-** + 2B%e~F*}, 
x 
(115) 
2y 
k(x) = aa (a2x?+ 3ax+3)e- 
— (67x? + 38x+3)e-#*}, 
where 
. M 
° 2M Bal)" 20ar| Bal" 
pe y (116) 


2 
1= 5 f°-M(M | Eo})}. 
The normalization condition (101) becomes 
f (u?+v?)dx=1, (117) 
0 


and the quadrupole moment (112) becomes 


Q — ( ‘)a (118 
=> uv — ———V P 
io J, x v i x ) 


and it is now measured in units of (4M|Eo|)-. 

We have made an estimate of the constants 
by carrying out a variational calculation to 
obtain the minimum value for y, and approxi- 
mate expressions for u(x) and v(x) in order to 
calculate Q. Hulthén'‘ has shown that for the 
simple Schrédinger equation with the potential 
e~**/x, a good approximate solution is 


{a:(1 —e~*) +-a2(1 —e~*)*}e-*/8. 
Hence we have taken this expression for u(x), 


and an analogous expression for v(x). That is, 
we take 


u(x) = {a,;(1—e-*) +a2(1 —e-*)?} e-*/?, 


(119) 
v(x) = {b,(1 —e-*)*+-52(1 —e-*) 5} 


6 12 
x (14-45 ows, 
x x* 


4 L. Hulthén, Arkiv fér mat. astr. och fysik A28, No. 5 


(1942). 


(114) 





ns 


3) 


4) 


“4 








They have the required behavior at both small 
and large values of x. We insert these expressions 
in the energy which is 


1 fl , d*v Y* b(n) 
-= u—+v— —-— 4+ u7j(x 
4 0 dx? dx? x? J 


+ 2v2uvk(x)+-v*{ j(x)—k(x)} Jdx, (120) 


and determine the constants a), @2, b;, and bz so 
that y is a minimum and the normalization 
condition (117) is satisfied. The calculation is 
tedious but straightforward, and we give only 
the results. Taking | £o|=2.17 Mev, «=177 m, 
and p/x=2, we obtain 


7 =0.433, ; (121) 
a,= 3.18, a2.=> —2.01, 
b,=0.095, b2.=—0.072. (122) 


(121) gives (fx)?=0.375, and the evaluation of 
the quadrupole moment with (119) and (122) 
gives Q=1.6X10-*’ cm*. The latter seems to 
indicate that the ratio u/x should be larger to 
give the right value for Q. However, it is well 
known that the variational method may give 
quite good results for the quantity which is 
minimized (in this case f), but that the use of 
the approximate wave functions to calculate 
other properties may give very bad results. In 
our case this is quite possible since Q is propor- 
tional to v(x) which contributes very little to 
the energy, and moreover the main contribution 
to Q comes from the parts of the wave functions 
at large distances from the origin which do not 
affect the energy very much. In fact the same 
calculation made with u/x=3 actually gives a 
smaller value for Q. Thus it seems that it is not 
possible to obtain anything but a rough estimate 
u/x~2 without going into a much more careful 
computation. 

We can, however, show that the conditions 
given in the beginning of this part can be 
satisfied with our estimated values for the 
constants. Now since Ejso, = 3a/2f?, 


Evina 2 (A)= 0. 0056 
Eisop 3(ax) (ax) 


Thus there is a region of values for a (and 
consequently for Eis») for which all the in- 





(123) 
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equalities in the conditions given above hold by 
a factor of about 10. 


D. The Magnetic Moment 


The calculation of the magnetic moment of 
the deuteron is exactly analogous to that given 
in P-D for a single nucleon. The current vector 
is!5 
j=¢’V¢'—¢'V¢’? 

+(49)*f S04 oa(g¢'ra?— g? ra) U4 

+¥°X(V Xq')— YX (VX q’) 

+(4r)*f Dal(oaXy')ra? 
—(o4X*)ra'}Us. (124) 


Let us first consider just the part due to the 
pseudoscalar meson. Inserting in the expansion 
(46) the value corresponding to the absence of 
free mesons (g’*=0) and the lowest eigenvalue 
of the interaction energy (52), we obtain with 
the aid of (50) 

f Fe a 0X4 (125) 
~ (4m) ae or 
We have neglected terms in J;;/J in (50) since 
they only give rise to higher order terms. 
Putting (125) in the pseudoscalar part of (124), 
we obtain 


rh BOR) (a) 
=>— CAi ——i Cp; 
” 4a A, B, i, j ? Ox; J OX "ox; 
( 1 ~*) 0 2 OX, 
-(ee 
Ox; OX, Ox; 
OX OX 
~4r| (ci “ein (ci )ers| Us| 
Ox; Ox; 

















or pf 
: es a 1 2 
R= > {@a@pj—Catn;} 
4m 4. Bij 
OX, O°Xe OX 
tata +46 j,—Upz}. (126) 
Ox; OX,OX; Ox; 


The magnetic moment is given by 


é€ 
Mux J (xejr—xj)dV, (127) 





16 We have not written down here the terms which 
would occur if the vector meson had an anomalous mag- 
netic moment, but this will not change our results in any 


way. 
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and inserting (126) for j:, we obtain 


ef? i a 
Mu=— : i eT a bY (128) 


A, B, i,j 


where 


ij 0X4 0?X zp 0X4 
Au={ Xk +426 ;—— Us, V. 
Ox; Ox,0X; Ox; 
(129) 


In the double sum over A, B in (128), the 
terms with A =B give the moments of the indi- 
vidual particles, and they are just the expressions 
calculated in P-D. Since it was found there that 
the moments of the proton and the neutron are 
equal in magnitude and opposite in sign, these 
terms will cancel out, and we need only take the 


terms with A #B. Since A, is not altered by the 
interchange of particles I and II, (128) can be 
written 
ef? ii ai 
Mu=—- 2 Pi(Au—Au) (130) 


Wij 








with 
21 12 21 12 
Pi 5= lere1j—eriernjt+enierj—enierj}. (131) 


We see that P;; is antisymmetric in i, 7 so that 
it has only 3 independent components. In order 
to compute the component of M in a given 
direction, we need the diagonal matrix elements 
of the P;; in the (S, 7, M, N) representation. 
However the P;; change sign under the operation 
(85) so that they cannot have any matrix 
elements between states with the same signature. 
Thus this theory gives zero magnetic moment 
for the deuteron. , 

There is no need to continue the calculation 
for the vector meson since the reason that the 
magnetic moment vanishes in the above case is 
so general that it is clear the same result will be 
obtained in all cases of the strong coupling 
theory. 

Of course the existence of the tensor force 
gives a magnetic moment due to the orbital 
motion of the proton, but this contribution is 
very small (of the order of a few percent of the 
observed value). 


V. HEAVY NUCLEI 


It was shown in S-D that the charged scalar 
and the neutral pseudoscalar meson theories 
give rise to nuclear forces which have no satura- 


tion property. These results are due to the fact 
that in these theories, any number of nucleons 
can be arranged in a symmetrical configuration 
in such a way that each particle interacts with 
an attractive potential with all the others: in 
the charged scalar theory by having the iso- 
topic spin vectors of all the particles in the same 
direction, and in the neutral pseudoscalar theory 
by having all the particles lie in a plane with 
all their spin vectors parallel and perpendicular 
to this plane. In the neutral mixed theory such 
configurations do not lead to attractive inter- 
action between the particles.'° In the sym- 
metrical pseudoscalar and symmetrical mixed 
theories, we have a system of 3 orthonormal 
vectors e* instead of a single spin or isotopic 
spin vector for each particle. The only sym- 
metrical configuration in these cases is obtained 
by having all the vectors parallel, and the 
interaction for this arrangement is repulsive. 
Thus for these theories the potential energy will 
be proportional to the number of particles and 
not to its square as in the cases treated in S-D. 

However, for the forces to have the saturation 
property, it is also necessary that the kinetic 
energy of the particles increase faster than the 
attractive potential energy as the radius of the 
nucleus decreases in order that there exist a 
minimum for the total energy of the system. 
This condition is fulfilled in the usual weak 
coupling theory” since, if we denote by ro the 
radius corresponding to the volume occupied by 
each particle, the kinetic energy increases as 
ro*, and the dependence of the potential energy, 
of course, is the same as the behavior of the 
radial potential function, but this cannot increase 
faster than r~. In the strong coupling theory 
however, the existence of isobar states increases 
the number of allowed states in a given volume 
of space with the same upper limit for the energy, 
and this will change the dependence on fo of 
both the kinetic and potential energies. 

In order to see how this feature changes our 
results, we have made a calculation similar to 
that carried out in Bethe and Bacher. We make 


16 The second configuration gives attraction only when 
K(r)>J(r). For u/x~2, this condition holds only when r 
is greater than about 2/« and for such distances, both J(r) 
and K(r) are very small. 

17H. A. Bethe and R. F. Bacher, Rev. Mod. Phys. 8, 
82 (1936), in particular §25. 
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a Hartree approximation for the wave function 
of the system with free particle solutions for the 
individual wave functions, both for the spatial 
and for the spin parts. The latter are hyper- 
geometric functions and have been investigated 
by Rademacher and Reiche'* in connection with 
the quantum theory of a symmetrical top. The 
integrals which we need are all tabulated there. 
We shall not give the details of the calculation 
but just make a few remarks on the points which 
differ from the calculation given in Bethe and 
Bacher. 

The number of particles N is now the product 
of N; and N, where N; is the number of states 
with kinetic energy less than ko?/2M, and N, is 
the number of states with isobar energy less 
than (a/4f*){so(so+1)—2?}. The total energy of 
the system is 


E=E,+E,+E, (132) 


where E;, E,, and E, denote the kinetic, isobar, 
and the potential energy, respectively. If we 
consider so to be large compared to 1 and keep 
only the highest powers of it, then N;, N, and 
E;, E, depend on ko, So in the same way. Thus 
for (132) to be a minimum, 


ko?/2M = (a/4f?)s0?, (133) 
and 
E.=E,. (134) 


Moreover they are proportional to N/ro. E, on 
the other hand turns out to be proportional to 
N/(ro)*, and there does exist a value of ro at 
which E is a minimum. 

In the above consideration we have neglected 
the Coulomb energy, and though this may not 
affect the saturation property, it does bring in a 
serious difficulty in connection with the stability 
of heavy nuclei with high charge. As pointed 
out by Fierz,'® in a nucleus like that of uranium, 
the transition of a proton or a neutron to a 
particle with negative charge is energetically 
favorable unless the isobar separation is about 
50 Mev. However such a value for the isobar 
separation is inadmissible since we would then 
have to take the size of the source of the same 
magnitude as the Compton wave-lengths of the 

18H. Rademacher and F. Reiche, Zeits. f. Physik 39, 


444 (1926), and 41, 453 (1927). 
19 M. Fierz, Helv. Phys. Acta 14, 105 (1941). 


mesons. That is, the size of the source would 
have to be taken as the same order of magnitude 
as the range of the nuclear forces, and hence the 
shape of the source would completely determine 
the nature of the nuclear forces. Furthermore, 
such a value also removes the main reason for 
the introduction of the strong coupling theory; 
namely, to give a small value for the meson 
scattering cross section. 


VI. CONCLUSIONS 


Our results show that the theory considered 
here suffers from two grave difficulties; it gives 
a magnetic moment for the deuteron a value 
only a few percent of the observed value, and it 
predicts instability of highly charged nuclei. 
These results seem to be fundamental properties 
of all strong coupling theories, and there does 
not seem to be any way of overcoming them. 

These difficulties are not present in the weak 
coupling theory, and it thus seems advisable to 
go back and reconsider the arguments which 
led us to take up the strong coupling theory in 
favor of the weak coupling theory. The main 
difficulties in the weak coupling theories are the 
divergences due to the treatment of the heavy 
particles as a point source, and the large scatter- 
ing cross section of the meson. As already 
pointed out by one of us,” the first difficulty 
can be overcome by using the A-process developed 
by Wentzel and Dirac, and the second by using 
the theory of radiation damping developed by 
Heitler and Wilson. In addition, the weak 
coupling theory developed in this way has the 
advantage of relativistic invariance which the 
strong coupling theory does not have on account 
of the finite size of the source. Thus there is no 
reason now to consider the strong coupling 
theory, and we should go back to the weak 
coupling theory. 

As stated at the end of Part A of Section III, 
the transition from the strong coupling to the 
weak coupling theory can be made by a simple 
replacement of the spin and isotopic spin vectors 
in place of the vectors e*. Thus the radial wave 
functions for the ground state of the deuteron 
in the weak coupling theory satisfy exactly the 

20W. Pauli, Bull. Am. Phys. Soc. New York Meeting, 


Jan. 22-23, 1943, Abstract No. 25: Phys. Rev. 63, 221 
(1943). 
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same differential equations (114) except that as 
shown in the Appendix, 7 is redefined as 


y= 2f"?M(M|Eo|)}, 


where f’ is the coupling constant in this theory, 
and there is no longer any approximation due 
to the higher isobar states since they do not 
exist. Hence the computation in Part C of 
Section IV applies just as well for this case, and 
we obtain (f’x)?=$(fx)?—0.042 for the strength 
of the coupling. The ratio of the masses is 
unchanged. 


APPENDIX 
Calculation of the Matrix Element of e;* 


The e;* satisfy the following commutation relations with 
s; and é* [cf. P-D, Eq. (55) ]: 
[si, e;*]=te.%, i, 7, k, cyclic; 


i, eP] =ie;", a, B, Y cyclic. eames 


It is more convenient to introduce the imaginary com- 
ponents 
$4=51+152, S_=5,—152, S3; 
tt=f+i2, f=f-it, 2B; ee 
and 
e,* =e)!—e2+71(e2'+e;"), €,3=e,5+7e:', 
e,~ =e;'+e2+1(e2!'—e;"), e-?=e;5—1e2', 


e_+ =e;'+e?—i(es!—e:"), e3+=e3'+1e;7, (137) 
e_~ =e;!—e2?—i(e2'+¢;"), e3” =e3'—1e;’, 
e33. 


For the matrix elements of the e;* between the states 
(s, m,n) and (s’, m’,n’), we have the following selection 
rules: 


(1) for any e;*, ’=s—1,s, ors+1; 
(2) for e,%, m'=m—1; 

(3) for e3*, m’=m; 

(4) for e_%, m'’=m+1; 

(5) for e;*, n'=n—1; 

(6) fore?, . n'=n; 

(7) for e-, n’=n+1. 


Matrix elements which do not satisfy these conditions all 


vanish. 
The non-vanishing matrix elements of s; and ¢* are: 


(m|s,|m—1)=[(s+m)(s—m+1)}, 

(m|s3|m) =m, (138) 
(m|s_|m+1)=[(s+m+1)(s—m) }, 

with exactly the same expressions for /*, #, & with m 


replaced by m. Hence in order to satisfy the commutation 
relations (135), we must have the m dependence of the e;* as 


follows: 


(s, m, n\e,*|s—1, m—1, n’) 
= —(s, n|c*|s—1, n’)[(s+m)(s+m—1)}, 
(s, m, n|e,*|s, m—1, n’) 
=(s, m|c*|s, n’)[(s+m)(s—m+1)}, 
(s, m, n|e,*|s+1, m—1, n’) 
=(s, m|c*|s+1, n’)[(s—m+1)(s—m+2) }, 
(s, m, n|es*|s—1, m, n’) =(s, n|c*|s—1, n’)(s?— m)}, 
(s, m, n|e3*|s, m, n’) =(s, n|c*|s, n’)m, 
(s, m, n|es*|s+1, m, n’) =(s, n|c*|s+1, n’)[(s+1)?—m?}, 
(s, m, n|e_*|s—1, m+1, n’) 
=(s, n|c*|s—1, n’)[(s—m)(s—m—1)}, 
(s, m, n|e*!s, m+1, n’) 
=(s, n|c*|s, n’)[(s—m)(s+m+1)}, 
(s, m, n|e_*|s+1, m+1, n’) 
= —(s, m|c*|s+1, n’)E(s+-m+1)(s+m+2) }. 


Similarly to satisfy the commutation relations, the n-de- 
pendence must be as follows: 


(s, m|ct|s—1, m—1) = —(s|b|s—1)[(s+n)(s+n—1)}, 
(s, n|ct|s, n—1) =(s|b|s)[(s+n)(s—n+1) }, 
(s, m|ct|s+1, n—1) =(s|b|s+1)[(s—n+1)(s—n+2) }, 
(s, n|c|s—1, nm) =(s|b|s—1)[s?—n?}, 
(s, m|c?|s, nm) =(s|b|s)n, 
(s, m|c?|s+1, m) =(s|b|s+1)[(s+1)?—n?}, 
(s, n|c~|s—1, n+1) =(s|b|s—1)[(s—n)(s—n—1)}, 
(s, n|c~|s, n+1) =(s|b|s)[(s—n)(s+n+1) }, 
(s, m|c~|s+1, mn +1) = —(s|b|s+1)[(s+n+1)(s+n+2) }. 


To obtain the s-dependence, we use the normalization con- 
dition which the e;* satisfy [cf. P-D, Eq. (S1a)]. Using 
(137) we find 
e,%e_3+(e;3)?=1, 
e,te_~+e_te,-+2e;te;- =4, (139) 
and from these conditions, we obtain the equations 


s3(2s—1)(s|b|s—1)2+(s+1)3(2s+3)(s|b|s+1)*=1, 
$*(s—1)(2s—1)(s|b|s—1)?+5%(s+1)*(s|b|s)? (140) 
+ (s+1)2(s+2)(2s+3)(s|b|s+1)?=2. 


Then from the relation 


(s-e*) = —#8 (141) 
we have 
(s|b|s) = —1/s(s+1). (142) 
Thus (140) can now be solved, and we finally obtain 
(s|b|s—1) =1/s[(2s—1)(2s+1) }, (143) 


(s|b|s+1) =1/(s+1)[(2s+1)(2s+3) }. 
In particular for s=s’= 4, we have the relation 
(4, m, n|e*| 4, m’, n')=4(m, n\or*|m’,n’) (144) 


which gives the connection between the strong and the 
weak coupling theories; namely, the only difference for the 
ground state is the occurrence of a factor § in the interaction 
energy which merely changes the values of the coupling 
constants. 
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Forms of the beta-decay interaction are sought which are totally symmetric and totally 
antisymmetric to the interchange of any two of the particles involved in the emission. It is found 
that there is no symmetric interaction and one antisymmetric. The antisymmetric interaction 
allows transitions without change of parity and with a change of nuclear spin of either 0 or +1, 
including 0-+0. First forbidden transitions entail a change of parity and Ai=0, +1 or +2. 
Second forbidden transitions are characterized by Ai= +2, +3 and no change in parity. The 
contributions of electron-neutrino states of definite angular momenta are presented explicitly 
as an aid in the study of the influence of the nucleus on the light particle waves. The influence 


of the Coulomb field is discussed. 





I, INTRODUCTION 


N the first paper describing a field theory of 

beta-activity Fermi! drew an analogy with 
electromagnetic theory and assumed the inter- 
action to be the scalar product of a four-vector 
constructed from the wave functions of neutron 
and proton with a four-vector constructed from 
electron and neutrino waves. It is well known 
that besides the four-vector one can construct 
from products of two wave functions an axial 
four-vector, a skew-symmetric tensor, and two 
scalars, one of which changes sign under reflection 
and is commonly called a pseudoscalar. Invariant 
products of any of these five covariant quantities 
with the corresponding contravariant quantities 
formed from electron and neutrino wave func- 
tions present a fivefold infinite manifold of 
choices for an interaction between heavy and 
light particles. The experimental material that 
may make possible a particular choice of an 
interaction consists of (1) the selection rules 
obeyed in beta-emission and (2) the shape of the 
beta-ray spectrum, particularly the shape of the 
spectra in forbidden transitions. 

An alternative choice of interaction has been 
proposed by Gamow and Teller? who noted that 
the selection rules imposed by the four-vector 
interaction were inconsistent with the beta- 
activity of ThCC’. The same violation of Fermi’s 
selection rules appears to be made in the beta- 
activity of He® and kindred nuclei. The Gamow- 


* Now at the Geophysical Laboratory, The Carnegie 
Institution of Washington, Washington, D. C. 

1E. Fermi, Zeits. f. Physik 88, 161 (1934). 

2G. Gamow and E. Teller, Phys. Rev. 49, 895 (1936). 


Teller selection rules are less stringent than the 
Fermi rules in that they allow Ai=0, +1 (except 
0—0), whereas the Fermi rules allow only Ai=0 
(including 0-0). Interactions that give G-T 
rules are the tensor and the axial vector inter- 
actions. 

Allowed transitions are defined as those transi- 
tions that are independent of the small com- 
ponents (i.e., the particle velocity) of the Dirac 
waves for the nuclear particles and in which the 
electron and neutrino waves have no node 
through the nucleus. The shape of the allowed 
beta-ray spectrum is the same for any one of the 
five possibilities of interaction listed above. It 
may be slightly different for linear combinations. 
Only in forbidden transitions will the spectra of 
the five possibilities be different. 

Uhlenbeck and Konopinski* have compared 
the shapes to be expected from each of the five 
invariants with the experimental results on Na™, 
P®, and RaE. Their conclusion is that there is 
evidence against all these interactions taken 
singly except the tensor interaction. It is possible 
that actually none of these interactions is suf- 
ficient by itself to explain the experimental 
results but the tensor interaction cannot be 
excluded at present because of the great flexi- 
bility of this interaction which contains several 
unknown matrix elements. It is also possible that 
a linear combination of the five possibilities will 
account for experiment although no one taken 
singly is able to do so. 


3 E. J. Konopinski and G. E. Uhlenbeck, Phys. Rev. 60, 
308 (1941). 
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The supposition that a particular one of the 
five invariants discussed might be a complete 
description of fact is a kind of extension of 
Fermi’s analogy to electromagnetic theory. It is 
what we would reasonably expect to obtain if 
beta-decay occurs in two steps, i.e., the proton 
emits a meson which then creates a positron and 
a neutrino. In this description one would natu- 
rally form covariant quantities with proton and 
neutron waves and contract them with the 
analogous light particle functions. The analogy 
with electromagnetism then leads to the hope 
that one of these covariants is enough. If, on the 
other hand, the beta-process is not analogous to 
radiation theory in this sense, linear combina- 
tions of the five invariants are as acceptable 
possibilities as the invariants taken singly. 

It is the purpose of this paper to present the 
theoretical predictions of a particular linear com- 
bination of the five invariants. The basis on 
which this combination is chosen is as follows: 
The beta-process is considered to be an essen- 
tially four-particle process, and the requirement 
on the interaction is that it shall be either totally 
symmetric or totally antisymmetric to the 
interchange of role of any two of the participating 
particles. 


Il. THE INTERACTION 


The symmetry properties of the interaction 
between fields of the four particles are most 
easily determined if the beta-process is described 
as the simultaneous creation of neutron, proton, 
electron, and neutrino. Since all these particles 
have spin 3h, this description can be permitted 
by assuming the Dirac wave equation and the 
hole theory extension of that equation for all four 
particles. A hole in the theoretical neutron sea 
of negative energy is then the experimentally 
observed neutron. In negative beta-emission 
such a neutron hole is filled at the same instant 
that a proton appears to take its place in the 
nucleus and electron and neutrino are created. 
Positive beta-emission is then a simultaneous 
disappearance of four particles: a proton of 
positive energy and ‘‘neutron,” electron, and 
neutrino of negative energy. In this description 
it is evident that symmetrization of the inter- 
action producing beta-activity is identical with 


symmetrization of an invariant product wave 
function of four particles of spin one-half. 

Four particles of spin $4 cannot form a totally 
symmetrical state that is invariant to Lorentz 
transformations. This may be seen at once with 
the aid of the calculation rules of Van der 
Waerden’s spinor calculus.‘ If two of the four- 
wave functions have dotted indices and two 
undotted, the invariant forms are not sym- 
metrical because of the skew-symmetry of the 
e-operator. If all indices are undotted (or all 
dotted), the symmetric invariant vanishes be- 
cause of the identity: 


eH err + ee + eke = 0. 


The symmetrization in this discussion is re- 
stricted, of course, to space and spin coordinates. 
If isotopic spin were introduced, symmetrization 
would have no significance. 

There is one and only one totally antisym- 
metric invariant wave function constructed of 
the four Dirac y-functions, viz., 


VM, &% wi gi 
V2 D2. Yo 
x(x) = , (1) 
Vs; 3; vs $3 


Ve Dy Ws dalx 








Different Greek letters are used to distinguish 
between particles and, in the following, WV shall 
stand for the neutron wave function, ® for the 
proton, y, electron, and ¢, neutrino. All are taken 
at the same point in space, x. The subscripts in 
Eq. (1) refer to the four components of a Dirac 
wave equation, the particular form of which is 
immaterial as long as it is the same for all par- 
ticles. That x is an invariant is readily seen from 
the fact that the determinants of the matrices 
that represent Lorentz transformations are 
unity. Let S be the unitary matrix that repre- 
sents a certain Lorentz transformation: 


y= Sy. (2) 


This S applies to all four wave functions in the 
determinant, Eq. (1), so if we write 


x=|DI, 


*B. Van der Waerden, Géttingen Nachrichten (1929), p. 
100. 
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then the new antisymmetric wave function, x’, is 
x’=|S-D| =|S|X|D|=|D|=x, Q.E.D. (3) 


For the sake of definiteness the numerical sub- 
scripts in Eq. (1) will be assumed to refer to the 
four rows and columns in their usual order in the 
a and 6 matrices in the familiar form of Dirac’s 
equation: 


ha 
[i —— (a, p) —Ame |W =0. (4) 


c ot 
If we consider the wave functions in Eq. (1) 
to be quantized, i.e., as operators, the form of the 


interaction between fields of the elementary par- 
ticles becomes: 


H' =¢ J x(x) +x(x)* dx, (5) 


where g is the “Fermi constant’’ having dimen- 
sions of energy times volume, and x(x) gives rise 
to the disappearance of four particles at the 
point x; x(x)* gives rise to their appearance. The 
transitions induced by H’ are between states of 
equal energy and therefore the probability of 
beta-emission is given by the well-known formula 


P(E)dE=2rh"| H’'|*p,(w—E)p(E)dE. (6) 


In Eq. (6) p»a(w—Z£) is the density of neutrino 
states of energy w—E and p,(E) the density of 
electron states of energy E. The energy released 
by the nucleus is w. 

An obvious difference between the form of this 
theory and that of preceding ones is that H’ is 
an integral over products of YW and ® or over 
products of ¥* and &*. This is only a formal dif- 
ference from the customary method of calculating 
matrix elements which contain ¥* and ® to- 
gether, etc., but it will be convenient all around 
to introduce ¥* in place of ¥ in Eq. (1). Then the 
physical interpretation of Eq. (5) will be the 
usual disappearance of a neutron when proton, 
electron, and neutrino are created or vice versa. 

The replacement of ¥ by ¥* at the same time 
that positive and negative energy levels are 
exchanged has been studied by Furry® and 
Pauli.* These authors have shown that there is 


5 W. H. Furry, Phys. Rev. 51, 125 (1936); 54, 56 (1938). 
®W. Pauli, Inst. H. Poincaré Ann. 6, 130 (1936). 


a unitary operator @ such that 
y*—0 (7) 


when positive and negative energy levels are 
interchanged. In the spin coordinates pertinent 
to Eq. (4) 

6= —iBar. (8) 


Applying @ to the neutron waves in Eq. (1) and 
substituting in Eq. (5), we get the interaction for 
positron emission: 


Wi(x)* Pi(x) ilk) ¢:(x) 
——— f —W3(x)* Do(x) a(x) 2(x) 

—W2(x)* 3(x) ws(x) $3(x) 
Wi(x)* a(x) Pa(x) ga(x)| (9) 


dx. 








Electron emission is given by the complex con- 
jugate of Eq. (9), but since in the expression for 
the probability of emission only |H’|? appears. 
it is sufficient to consider Eq. (9) alone. 

It is convenient to expand Eq. (9) for H,’ as 
a linear combination of the familiar invariants. 
For this purpose we shall use the Dirac operators, 
p and a, acting on the heavy particle waves only 


a= 10: QtyO2, a: hae Baziy&s, be B, 























(10) 
Or =MtyAz, Ty=1Wz, Tz=WizpMy. 
In this formulation Eq. (9) becomes 
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From this expansion selection rules can be 
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deduced at once. The terms proportional to the 
o’s and to p; constitute an axial vector inter- 
action and give At=0, +1 (0-0 excluded) with 
no change in parity in allowed transitions. The 
ps term is the scalar interaction and gives Ai=0 
(0-0 included) no change in parity. The pz term 
is the pseudoscalar and gives no allowed transi- 
tions. 

A curious form of the rule for determining the 
change in parity is evident from Eq. (11). It will 
be noted that the light particle functions that 
are associated with p3 and @ (which mix large 
terms of proton waves with large terms of 
neutron waves) contain products of one small 
and one large component of the light particle 
waves. Thus, if we follow the customary method 
of giving the electron (and neutrino) the orbital 
quantum number L that applies to the spin 
eigenfunction, i.e., the same as the large com- 
ponent, we obtain the rule: There is no change 
in the parity of the nucleus if the sum of orbital 
quantum numbers of electron and neutrino is 
odd; the parity changes if the sum is even. This 
is just the reverse of the rule that applies to 
particles of integral spin and also of the rule 
that applies to the absorption and subsequent 
emission of a single particle of any spin. The root 
of the paradox lies in assigning a definite parity 
to wave functions of spin $4. The form of the 
matrix that represents inversion of the coor- 
dinate system® shows that the phase of these 
waves is changed by $2 rather than 0 or z, the 
phase of the large components increasing if that 
of the small decreases and vice versa. Thus two 
successive inversions lead to the original wave 
function but with opposite sign, the same result 
as that obtained by rotating the coordinates by 
360° about some axis. It is well known that there 
is no paradox in the case of the rotation because 
the space of all rotations, although continuous, 
is two-sided.’ 

A short discussion of the above results has 
been presented in a letter to the editor of The 
Physical Review.* It is pointed out in the letter 
that addition of scalar and axial vector operators 
increases the lifetimes of the 4n+2 nuclei relative 


7 Cf. E. Wigner, Gruppentheorie und ihre anwendung auf 
die quantenmechanik der atomspektren (Braunschweig, F. 
Vieweg & sohn, 1931), p. 99. 

(198i Critchfield and E. Wigner, Phys. Rev. 60, 412 
1941). 


to the 4n+1 nuclei in allowed transitions.® This 
improves the agreement of theory with experi- 
ment for C!® and F'* but makes it worse for He‘. 
Furthermore, the ratio of the transitions to 
excited and to ground states of the nuclei con- 
sidered by Grénblom® is predicted more closely 
by the antisymmetric theory than by the axial 
vector (or tensor) alone. The reason for this is 
again the relatively stronger transitions without 
change of spin which the antisymmetric theory 
predicts. 


Ill. MATRIX ELEMENTS 


The shapes of the forbidden spectra under the 
theory presented in this paper can be derived 
from the results of reference 3 on the scalar, 
axial vector, and pseudoscalar interactions. A 
slight extension is necessary to include some 
cross product terms in the first forbidden spectra. 
I have repeated the work of Konopinski and 
Uhlenbeck insofar as it applies to the antisym- 
metrical theory of beta-decay by the method 
which uses the eigenfunctions of both electron 
and neutrino in polar coordinates. Except for the 
addition of cross product terms between different 
kinds of interaction, the results are, of course, 
the same as in reference 3. Nevertheless an out- 
line of the method is presented, and the results 
are given in detail for several reasons. The chief 
reason is to have the contributions of various 
electron and neutrino angular momenta pre- 
sented separately because the possible influence 
of the nucleus on these contributions is very 
likely different. Beta-spectra may furnish a 
useful experimental means of determining the 
influence of nuclei on electron and neutrino 
waves. It appears certain from the work of Rabi 
et al.° that there is no spin dependent force 
between electron and nucleus except that due to 
magnetic moments; but it may well be that there 
are spin-independent forces in addition to the 
Coulomb field. Another reason for presenting the 
calculations in polar coordinates is to throw light 
on the detailed effect of the Coulomb field on the 
shape of the beta-spectrum and also on the life- 


* E. Wigner, Phys. Rev. 56, 526 (1939), and in the same 
issue, B. Grénblom, Phys. Rev. 56, 508 (1939). 

1 Rabi, Millman, Kusch, and Zacharias, Phys. Rev. 55, 
526 (1939); Kusch, Millman, and Rabi, Phys. Rev. 57, 
765 (1940). 
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times of the beta-activities. The latter questions 
are taken up in Section IV. 

In this section it will be assumed that there is 
no interaction between the electron and the elec- 
tromagnetic field of the nucleus. We have, 
therefore, to use the eigenfunctions of kinetic 
energy for both electron and neutrino in Eq. (11) 
and substitute H,’ in turn in Eq. (6). 

The “matrix element” of beta-decay theory, 
M, is related to our H’ by 


M=—i(2xh?R/gp.p,)H’, (12) 


where R is the radius of the sphere in which 
quantization is effected, p, is the momentum of 
the electron, and p, is the momentum of the 
neutrino. Using M in place of H,’ we find the 
probability of emission of a positron and neutrino 
(absorption from states y and ¢, respectively) 
with total angular momentum jh and 2-com- 
ponent mh to be 


P(Y, ) j,m= (g?/20%c*h’) Ep.(w— E)? 
X | Mv, b)im|*. (13) 


An electron and neutrino emitted into a state 
specified by 7 and m may cause a maximum 
change in the nuclear spin of jh. The nuclear spin 
may also change by integral multiples of h less 
than jh provided that the total angular mo- 
mentum is conserved. It is necessary to calculate 
the matrix element only for one particular m, 
and it is convenient to choose the highest value, 
m=j. The total probability of emission into 
(Wy, @); is the sum over all m, and this sum will 
be considered automatically accomplished in the 
square of matrix elements. 

Consider the particular case in which an elec- 
tron disappears from an S, state of energy —E, 
and the neutrino from a P, state of energy E—w. 
If the total angular momentum of these two 
particles is zero, the state will be denoted by the 
symbol 

(SyP4)o. 


The electron wave will always be given first. A 
more expanded form of the wave function is 


2-175, P,P — SOP PH), 


where superscripts denote the values of m for 
individual waves. Now we may calculate 
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M(S;, Pi)o by substituting this form into Eqs. 
(11) and (12) and evaluating the expressions at 
the nuclear radius. Only the leading terms in the 
series expansions in powers of r/X, or r/A, will 
be retained. Under these conditions 


M(S,P\)o=}a f a (14) 


Here fp; stands for /V*p;6dx. In general 


faq fvraaar. 


The factor a and the related factor 6 which 
appears later are defined by 


a=[(E+mc*)/2EF }}, 


15 
b=[(E- mc)/2E}!. 9) 


In some expressions a’, a’’, b’, and b” will be 
used for those a’s and b’s that are changed by 
the effect of the Coulomb field. Their values in 
the absence of a field are the same as the un- 
primed letters. 

The calculation for the state (P45;)o yields 


M(P,S\)o= —}b f = (16) 


In order to simplify the presentation of the 
formulas, we shall combine related matrix ele- 
ments into one equation when feasible. Thus 
Eqs. (14) and (16) may be condensed into 


M(S,P;)o/a= —M(PiS\)o/b=4 fo» (17) 


The allowed transitions may also cause emis- 
sions into states of unit angular momentum for 


which: 


M(S,P;):/a= -M(P\S)/b=4 fe (18) 


Matrix elements (17) and (18) comprise all 
the allowed transitions, and it is evident from 
Eq. (15) that the sum of squares of these matrix 
elements will be independent of light particle 
wave-lengths. We now present the matrix ele- 
ments that lead to first forbidden transitions. 
The latter are of two types: dipole matrix ele- 
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ments which are proportional to p3 and @ but also 
contain the first power of the momentum of one 
of the light particles, and pseudomonopole matrix 
elements, proportional to velocity terms of the 
heavy particle waves, i.e., to p; and ps2. In both 
cases there is a change in parity, so it is necessary 
to consider all possible changes in angular mo- 
mentum: 


M(S\S\)o= }ia f —s 


+ (ap, +b'p,)/6h f (e, r), 


M(P,P,)o=}ib f _" 
(19) 
~(bpat+a’p.)/6h f (c, r), 


M(S,5,):=(ap,—b'p.)/Oh f [osr—i(eXr)], 


Tei 
tj 


M(P,P,):= —(bp.—a'p.)/6h f Cosr-+i(eXr)], 


M(S,D1;)1/apn= —M(Dy5))1/bp. 


=M(P;P3):*/bp» 


—M(PuP;):*/ap. 
= (02/128) f [2psr+i(oxt)] 
(20) 
M(S,D4)2/ap, = —M(D1yS;)2/bp. 
= —M(P;P)2/bp,, 


=M(Pi,P;)2/ap. 


= (31/6h) f [or], (21) 
where 


[or Jic=4(oretonr,) —i(o, r) dix. 


Equations (19), (20), and (21) give all the 
matrix elements that lead to first forbidden 
transitions. Following are the results for second 
forbidden transitions. In these p;[rr] and 
[r(@Xr) ] are defined analogously to [er] above, 
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and the symbols [err] and Q will be used: 
Corr )ie=}(orpritosaretorurir,) 
1 


——(er?) ——r(a, r), 
15 5 


Q= J { psLer]—iLr(o xx) ]} /(12) in", 


(22) 
M(SjP1;)2= — (b'pb.pn/3 —ap,?/5)Q, 
M(P1,5;)2= (ap.p,/3 —b" p2/5)Q, 
M(P;Dy)2 = (a’'p.p,/3 — bp,?/5)Q*, 
M(DyP3)2= — (bp.pn/3 — a" p2/5)Q*, 
M(SjF25)2/apn?= —M(F5S;)2/bp2 
= —M(P,D2)2*/bp,? 
=M(D2,P;)2*/ap2 
= f (3oCrr] 
+2i[r(@ Xr) ]}/30v2h?, (23) 
M(P,,Dy;):/a= —M(DyP)2/b 
= (boba/6-3%) f psltr], (24) 


M (Sj F24)3/apn? = —M(Fo,5;)3/bp2 
—M(P,D2)3/bp,? 
=M(D2yP;)3/ap2 


=W-*(120)-4 f [Corr] (25) 
M(PyDy)3= —M(DyP)s 
= (psp./6h4) { Corr}. (26) 


IV. EFFECT OF THE COULOMB FIELD 


The Coulomb field of the final nucleus has an 
important effect on the beta-spectrum of most 
forbidden transitions. In fact, it is only for the 
very light nuclei and (or) rather large electron 
energies that the effect of the nuclear charge can 
be neglected. This has already been pointed out 
by Konopinski and Uhlenbeck.* The effect may 
be described as follows: In the absence of a 
field, electron waves may be classified as S 
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waves, P waves, etc., according to the trans- 
formation properties of the “large components,” 
the radial dependence of these waves being one, 
r, etc., (r” in general) near r=0. In the presence 
of the field, the angular dependence remains 
unchanged, but the radial dependence acquires 
two new properties. The first is that some of the 
solutions (total angular momentum 34h) are not 
bounded at r=0. The integral of the square of 
the wave function is bounded, however, and in 
calculating transition probabilities for beta- 
activities the value of the electron wave at the 
surface of the nucleus is used. This singularity 
is not important to the shapes of the beta-spectra. 

The second, and more important, effect of the 
electric field is to lower by unity the power of r 
with which a wave approaches r=0. Neglecting 
the small effect mentioned above, a P; wave ap- 
proaches r=0 for the most part as r itself, but 
there is a small admixture of a wave (with the 
same angular dependence) which approaches 
r=Oasr°, i.e., asa finite constant. The coefficient 
of this S part of the P, wave is proportional to 
Z/137 and is usually very small, but the more 
favorable radial dependence of this part near the 
origin makes it of importance compared with the 
P part itself. The physical interpretation of the 
electron state is that there is a finite probability 
of finding an electron extremely close to the 
center of attraction in spite of the fact that the 
electron has one unit of angular momentum. In 
the relativistic theory of the electron the increase 
in kinetic energy that arises upon concentrating 
the electron within a radius r of the center is 
proportional to ch/r. On the other hand, the 
decrease in potential energy in the Coulomb field 
is also proportional to 1/r so that it is possible 
to have the electron very close to the nucleus 
although it has angular momentum. 

So far as we carry the analysis of forbidden 
spectra in this paper it is sufficient to consider 
the amount of S part in the P waves and of the 
P part in D waves. There are P waves in both 
S, and P, eigenfunctions, however, and the coef- 
ficient of the S part will be different in the two. 
In order to determine the coefficients with which 
these anomalous parts of the wave appear the 
work of Dirac! on the radial functions is used. 


uP. A. M. Dirac, The Principles of Quantum Mechanics 
(Oxford, 1935), p. 265. 





Neglecting quantities of the order (2/137)? com- 
pared with unity we find the results are ade- 
quately expressed as modifications of the factors 
a and 6 in Eqs. (19) and (22) as already indicated 
by the use of primes. All other matrix elements 
are unaffected. The values of the primed letters 


are: 
a’ /a=1—3Ze?/2(E+mc) p, 


b’/b=1—3Ze/2(E— mc) p, 
a” /a=1—5Ze/4(E+mc)p, 
b’’/b=1—5Ze?/4(E— mc’) p, 


(27) 


where p is the radius of the nucleus at which the 
wave functions are evaluated. 


V. SHAPES OF SPECTRA 


The contributions of the individual matrix 
elements found above may be changed in dif- 
ferent ways by further, non-Coulombic interac- 
tion between electron and nucleus. In this section 
the electric field only will be taken into account. 
It will be further assumed that the nuclear wave 
functions are eigenfunctions of the time reversal 
operator” so that the operators in Eq. (11) that 
change sign under time reversal, @ and pe, lead 
to different states of the final nucleus than the 
operators that do not change sign, p; and ps3. 
This assumption eliminates most of the cross 
products of different operators in the squares of 
matrix elements. 

The probability of allowed transitions is cal- 
culated by adding the absolute squares of all 
matrix elements leading to such transitions. The 
positions of the nuclear particles do not matter 
to those matrix elements, Eqs. (17) and (18), and 
no change in parity takes place in the transition. 
Such transitions may thus be called monopole 
radiations. From the values of a and 0 as given 
in Eq. (15), the sum of squares of matrix ele- 
ments, Eq. (17), that apply if initial and final 
nuclei have spin zero is 


2X |M|*=4|Sps|* 0-0. (28) 


Similarly, if the nuclear spin changes by one unit 
the appropriate sum of squares is 


> |M|*=3| fel? At= +1. (29) 


12 E. Wigner, Géttingen Nachrichten (1932), p. 546. 
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If the transition involves no change in nuclear 
spin and the spin is not zero, the sum of (28) and 
(29) is to be used: 


D|M|?=2| Ss|?+2| So|? Ai=0 (60). (30) 


Since Eq. (30) represents a larger factor in the 
probability of decay than Eq. (29), transitions 
without change of spin will be more prominent 
than under the axial vector or tensor theory 
alone and this is desirable in Grénblom’s analysis 
mentioned above.’ Equation (13) represents the 
normal Fermi spectrum of electrons if M? is 
constant so that allowed transitions in the anti- 
symmetric theory give normal spectra. 

First forbidden transitions are conveniently 
divided into three classes: monopole, dipole, and 
mixed. The monopole radiations differ from the 
allowed emissions in that a change of parity is 
entailed, and for this reason the first forbidden 
monopole transitions are called pseudomonopole. 
The classification of beta-radiations as various 
multipole radiations is somewhat more con- 
venient than the broader: allowed, first, and 
second forbidden classification. This is especially 
true for light nuclei because the matrix elements 
characteristic of first forbidden transitions, say, 
may vary in magnitude over a much wider range 
than the difference between allowed and the 
strongest first forbidden.” 

Pseudomonopole matrix elements are propor- 
tional to p; and pz, i.e., to the small components 
of the nuclear waves. The two operators p; and 
p2 appear in Eq. (19) and, according to the 
assumption stated above, do not mix. The sum 
of squares of M applicable to pseudomonopole 
emissions is then 


| M|?=4| Sos|?+2| So 


The spectrum is normal. 

There are four distinct matrix elements that 
lead to dipole transitions. The sum of squares for 
each of these will be presented along with the 
allowed changes in nuclear spin. An indicated 
change of spin that is primed, as the 0’ in 
Ai=0', 1, means that the sum of spins of initial 


2 Ai=0. (31) 





18 Discussions of and formulas for the half-lives of beta- 
activities have been presented in three independent papers 
appearing in The Physical Review 61 (1942): C. L. 
Critchfield, p. 249 (light nuclei only); R. E. Marshak, p. 
431, and E. Greuling, p. 568. 
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and final nucleus must be as large as the unprimed 
number in the same Az. 


LX | M(a, r)|*=(1/36h?) 
X[(batcp?/E—3Ze*/2cp)? 


+m*c'p?/E*]| [(s, r)|*, 
Ai=0 (32) 
DX | M(psr) |?= (1/36h") 
X [(pa—cp2/E+3Ze?/2cp)* 
+2(p.2+p.2)+mictp?/E*) 
X|S psr|?, At=0',1 (33) 
> | M(e Xr) |?= (1/36h?) 
XO (pa—cp2/E+3Ze*/2cp)* 
+3(Pn? +p?) +mictp?/E* ) 
X|S(oXr)|2, Ai=0',1 (34) 


| M(or) |?=(1/12h)(p.2+p2) | f[or]|? 
Ai=0’,1',2. (35) 


Ordinarily the term containing Z? will far out- 
weigh other terms in a given expression and in 
every one except Eq. (35) the spectrum will be 
quite the same as the normal one. Deviations 
from the normal spectrum should be apparent 
in the high energy transitions in light elements 
if Ai<2 or in any case if Ai=2. 

The one mixed transition, between ps2 and 
(c, r) leads to 


2 | M(p2; (¢, 7)) |?=(1/6h) (Patcp?/E 
—3Ze*/2cp)| fp2|X|S(e,r)| At=0. (36) 


Again the Z term should predominate and the 
spectrum be normal. All first forbidden transi- 
tions cause a change in the parity of the nucleus. 

So-called second forbidden transitions can be 
of two types, pseudodipole and quadrupole. 
There are no pseudodipole transitions in the 
antisymmetric theory, however, because the only 
possible matrix elements of this type, pir and 
pet, do not permit At=2; and Ai<2, without 
change of parity, is allowed in monopole transi- 
tions. Thus all transitions are quadrupole, Ai=2 
or 3, in the second forbidden category. There are 
three distinct matrix elements and no cross 
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products under the above assumptions: 
DX | M(osrr) |? = (1/360h')(3p.4—4cpnp.*/E 
+10p.*p,? —4cpn*p?/E+3pn' 
+ (Ze*/cp)(6p.°—10cp,2p2/E 
+5Sp2pn—3p./E) + (15/8) (Ze?/cp)? 
X(4p.°+p2)]| Spsrr|? Ai=2, (37) 
DX | M(orr) |? = (1/360h*)(3p.'+ 10p.2p,? 
+3p,*)| forr|? Ai=2’,3, (38) 
LX |M(ro Xr) |?= {0 | M(osrr) |?/| S psrr|? 
—(1/3)20 | M(err) |*/| forr|*} 
| x|S(reXr)|? Ai=2. (39) 


If Z is fairly large or the energy particularly 
small, the normalization factor for waves in a 
Coulomb field should be included in the electron 
waves. An approximate form of the squares of 
the normalization factor is 


SO=F/(e—-1), €=2xZeE/hcp., (40) 


in which ¢ is positive for positron emission and 
negative for electron emission. The approxima- 
tion involved in Eq. (40) as well as throughout 
the calculations on the Coulomb effect is that in 
which quantities of the order (Z/137)? are 
neglected in comparison with unity. Our results 
thus apply only to light nuclei. The final form of 
P(E) may be written 


P(E) = (g?/2n°c*h") f(()Ep.(w-—E)?L|M|*, (41) 


in which the sum is taken over all elements that 
contribute to the transition. The integral of 
P(E)dE over all electron energies then gives the 
decay constant. These integrals have been pre- 
sented in the papers quoted in reference 13. 

The foregoing calculations have been carried 
through for positron emission. In order to apply 
them to electron emission, it is simply necessary 
to change the sign of Z throughout. 


VI. CONCLUSION 


Full application of the results of this paper 
would require many experimentally determined 
beta-spectra of forbidden transitions. Further- 


more, the spectra must be free of distortion such 
as introduced by thick sources, etc. The spectra 
of P® and Na*‘ as determined by Lawson" (cf. 
reference 3) are probably in the desired category. 
Both transitions are slower than allowed transi- 
tions, Na** by a factor 40 and P® by a factor 
2X10‘. Assignment of a definite multipole to 
these transitions is complicated by the fact that 
we do not know the magnitudes of the matrix 
elements and also by the intermultiplet character 
of the transitions (cf. Wigner, reference 9). Ac- 
cording to the estimates made of matrix elements 
in the first reference 13, Na®* could be a pseudo- 
monopole transition that depends upon p2. The 
Kurie plot predicted for such a transition would 
be a straight line, and in fact, the high energy 
end of the experimental curve does appear to be 
quite straight. The low energy end is probably 
complicated by another transition. 

The high energy end of the Kurie plot for P® 
is also straight but its lifetime is too long to 
ascribe the decay to the matrix element M(pze). 
In addition, the spectrum is probably simple so 
that the deviation at the low energy end cannot 
be explained as the influence of a second spec- 
trum. The lifetime of P® suggests the matrix 
element M(p;) according to the estimates in the 
first reference 13, but in order to obtain the ob- 
served spectrum, it would be necessary for the 
matrix element to become appreciably dependent 
upon the energy of the electron when this is low. 
This or any other elaborate explanation is to be 
avoided at this stage of the theory of forbidden 
spectra. 

Most of the work presented in this paper was 
done as a National Research Fellow at Princeton 
University and The Institute for Advanced 
Study. I am glad to acknowledge my debt to 
Professor Eugene Paul Wigner of Princeton Uni- 
versity for inspiration and assistance in develop- 
ing this theory and to Professor Furry of Harvard 
University for helpful criticism. The whole 
investigation is the direct result of the Seventh 
Annual Conference on Theoretical Physics 
sponsored by The George Washington University 
and The Carnegie Institution of Washington. 


4 J. L. Lawson, Phys. Rev. 56, 131 (1939). 
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On the Silicon Oxide Bands i 
i 
L. H. Woops s 
Ryerson Physical Laboratory, University of Chicago, Chicago, Illinois | 
(Received March 9, 1943) 
Pankhurst has described a band spectrum believed by him to belong to an oxide of silicon, ‘ 
possibly SiOz. A better excitation of the same bands has been obtained from a similar source (a a 
high voltage uncondensed discharge through a constriction in a quartz tube), but in helium gas t 
instead of hydrogen. These bands have been photographed on the 30-foot, 30,000 lines per inch c 
icago grating spectrograph. and near A as nm resolved in the first and secon 
Chicag ing spectrograph. A band 3840 has bee Ived in the fi d d 
orders and found to be a (0,0) transition, overlapped by a weak (1,1) transition, of the type 
2=-—»*2, having the constants 
»(0,0) = 26,015.05 cm-!, »(1,1) = 25,991.44 cm—, . 
By’ =0.7180 cm™, B,’ =0.704 cm“, V 
Bo” =0.7253 cm, B,” =0.712 cm™. V 
The coefficients of the spin doubling for the two states are yo’ = +0.012 cm™, yo’’ = +0.002 or C 
vo’ = +0.022 cm™, the value of yo’ not certainly fixed as between these two alternatives. The k 
doublet structure and the B values prove that the emitter is SiO*. Other bands at \4270 have t 
been resolved with weak intensity and tentatively ascribed to SiOz. P 
t 
c 
INTRODUCTION bands at 4270A have been resolved with weak r 
N 1940, R. C. Pankhurst! described a system wegen d and pronepin Teng = SOs; and ; 
of bands lying in the violet region of the still others a ys way oe ashe none ' 
spectrum, obtained from a quartz capillary in a nation. It wae uaaih to oe ‘ — yn 7 
hydrogen discharge tube through which passed a ments = as atmosp on sal ‘ Papo Sh Wane I 
heavy current. Part of the system, the region the a pooh — “on staan Wee ,; 
4400A to 4600A, had been previously described than is t eg : alae — to = thee mt 
by Cameron? and attributed by him to an oxide other _s ods peace = a v4 vor ' 
of silicon. Pankhurst’s published photographs Petter Photographs of the. weaker Sands, ut 
show the spectrum to be a series of more or less ‘®S€ a ave mecemmArey Deen Mmcennitely 
diffuse regions, the most intense of which had the Postponed. 
wave-lengths 4340A, 4270A, and 3840A. He con- EXPERIMENTAL i 
cluded from the general complexity of the spec- ; , 
om Get te sho aaa tee the wbeaet The general design of the apparatus was like c 
molecule SiO ; that used by Pankhurst, but with modified c 
a : ' electrodes. The bands were obtained by passing c 
The general interest in polyatomic molecules an alternating current of 2 amperes at 4000 volts c 
stimulated the present investigation of the bands, roe : eee ite Wetenentind atud : 
” Legend because — ron _ ap electrodes were fitted into the glass discharge s 
par ae ong a “go were par h tube against glass flanges and sealed with Picein. c 
os 7 on oman * A S iieen = 1g A quartz tube about 1 cm diameter at each end t 
a —e he e 30-foot Chicago grating but constricted in the middle of its length to a s 
ave made it possible to resolve the rotational 2-mm capillary was placed between the elec- j 
structure of these bands and to find their true trodes. Concentric with the axis of the quartz t 
emitters. Unfortunately, however, this investiga- capillary, there was a hole through each electrode P 
tion is incomplete. A band at 3840A was excited , permit the light from the capillary to pass to ‘ 
with great intensity in helium, and has been the window of the tube. The quartz tube, itself, f 
analyzed and positively assigned to SiO*+; other was held in position along the axis by two flat - 
1 R. C. Pankhurst, Proc. Phys. Soc. 52, 707 (1940). doughnut shaped steel disks, on the outer edge of f 
2 W. H. B. Cameron, Phil. Mag. 3, 110 (1927). which was a sleeve fitting snugly against the 1 
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inside of the glass discharge tube. Likewise the 
inside circumference of each ring extended in a 
sleeve slightly smaller in diameter than the 
quartz tube, which reached into the tube and 
held it firmly. These snugly fitting steel sleeves, 
devised to prevent the discharge from leaking 
around the quartz tube, forced the current to pass 
through the capillary and heat it. The whole dis- 
charge tube was cooled by immersion in a water 
tank. 

Although Pankhurst excited the spectra in an 
atmosphere of hydrogen, the present experiments 
were made in helium so that if the same spectrum 
were found it would be more certain that its 
carrier was in some way related to quartz. Also 
helium gives fewer atomic lines than hydrogen in 
the visible region, and interferes less with the 
silicon oxide bands. The tube was run for a long 
time with frequent charges of fresh helium to 
clean it, for these bands do not appear as long as 
many impurities are present. Even after many 
hours of operation, the tube had to be flushed 
with helium about every twenty minutes during 
an exposure. Hydrogen and water were the most 
persistent impurities. 

When the tube was clean, the capillary glowed 
brilliantly yellow in the helium atmosphere and 
remained red hot for several seconds after the 
current was shut off. The capillary of a quartz 
tube which had been used for a long time became 
coated with sputtered electrode metal, easily 
removable with acid, and with a glassy brown 
material which resisted HCL, H2SO,, and me- 
chanical scraping. Even aqueous HF did not act 
on it rapidly. It disappeared, however, if the 
capillary was heated in an oxygen flame, and the 
capillary became clear quartz again. This brown- 
ish material is thought to be SiO, an hypothesis 
supported by the disappearance of the brown 
color inside the capillary when the outside was 
heated to high temperatures. The identity of the 
substance was suggested by a review in Chemical 
Abstracts of an article? which could not be ob- 
tained in this country. The authors, Zintl and 
others, report that a bright brown powder, whose 
composition had the formula SiO, was condensed 
from vapor sublimed from a mechanical mixture 

* Zintl, Brauning, Grube, Krings, and Morawietz, Zeits. 


f. anorg. allgem. Chemie 245, 1-7 (1940); Chem. Abs. 35, 
1719 (1941). 


of SiO, and Si under conditions such that the 
components themselves were not volatile. Slow 
condensation of the vapor on a hot surface, con- 
ditions which obtained in the discharge tube used 
here, produced a “‘brown-black brittle glass or 
shellac like membrane”’ also of the composition 
represented by the formula SiO. The authors 
further report that both the powder and the 
glassy material acted like a compound in that 
each had physical and chemical properties differ- 
ent from those of the parent substances or of the 
mixture SiO.+Si. 

The same review mentions the reaction of SiO 
at room temperature with atmospheric oxygen, 
to form SiOz, a pyrophoric reaction if the SiO is 
powdered. This information indicates that the 
free energy of the reaction SiO0+40,—SiO, 
is much greater than that of the reaction 
CO+40.—COz;; that is SiOz, is much more stable 
against dissociation into SiO than is CO, against 
dissociation into CO. That SiO, is more stable 
than CO, against dissociation into the elements 
is indicated by the heats of formation. Since the 
spectrum of CO, is readily excited in discharges 
at low currents, it is quite logical that the spec- 
trum of SiOz could be excited under the same 
conditions. The great difficulty in exciting spectra 
of polyatomic molecules in discharges, namely 
that the high energy electrons break the mole- 
cules less stable than CO, into diatomic frag- 
ments, might be no difficulty in this case. 

After exposures varying from 1 to 30 seconds, 
on a Hilger E-1 spectrograph, showed the spec- 
trum to be clean and the band at 3840A to be 
very intense, an exposure of five hours was made 
on Process plates in the first and second orders of 
the 30-foot, 30,000 lines per inch, grating spectro- 
graph. The 3840A band was of sufficient intensity 
to analyze, but the rotational structure of the 
several bands at 4270A was barely discernible. In 
hydrogen gas Pankhurst observed the 4270A 
region to be more intense than the 3840A region 
by a ratio 5:3. Since the excitation potential of 
helium is higher than that of hydrogen, electrons 
in a helium discharge have greater energy than in 
a hydrogen discharge. That the 4270A bands be- 
come less intense in helium than the 3840A band 
of SiO* then indicates that these bands are not 
emitted by an ionized molecule. The spacing of 
the rotational structure in these bands is com- 
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Fic. 1. Photograph from Hilger E-1 spectrograph 
showing spectrum of helium and some impurities, chiefly 
water. The small inner spectrum is that of the excited 
molecules in the quartz capillary. The *2-»*2 band of SiO 
at 43840 in helium is much more intense than the group 
of bands at 44270. 


parable to or less than the spacing of the R 
branch in the 3840A band of SiO* where it is 
approaching a head. While it is possible that the 
4270A bands belong to neutral SiO, it is the 
belief of the present author, based on the fore- 
going evidence and supported by the pyrophoric 
formation of SiO, from SiO, that they are 
emitted by the triatomic molecule SiOz. 

The plates taken on the Hilger E-1 show a 
possible band at about 6530A of about the same 
intensity as the 4270A bands, but no attempt has 
been made yet to photograph it on the big 
grating. These plates also show, in the region 
3650 to 4270A, indications of the several bands 
reported by Pankhurst, but much less intense 
than the 4270A bands. It is easy to distinguish 
between the bands of silicon oxides and the re- 
maining spectrum on the plates from the small 
spectrograph, because the spectrum from the 


capillary is distinct from the spectrum of the 
wider part of the quartz tube. When light from 
the discharge was focused in a spot on the slit of 
the Hilger E-1, light from the capillary formed a 
sharp bright pinpoint of light in the center of the 
spot. On the photographic plate there were then 
two spectra: (1) a wide spectrum of helium lines 
and OH bands through the center of which ran 
(2) a narrow spectrum of bands of silicon oxides 
and lines of silicon. See Fig. 1. The group of 
silicon lines at about 2510A and the many 
ultraviolet bands of neutral SiO appeared 
strongly in the spectrum from the capillary. 


ANALYSIS 


The band at 3840A is a (0,0) vibrational transi- 
tion of the type 22-2, plus a corresponding (1,1) 
transition. The spin doublet structure is unre- 
solved at low K values, then increasingly widens 
in the P branch. See Fig. 2. In the R branch, the 
doublet structure widens but the rotational 
structure closes up as it approaches the band 
head, so the spin components for successive 
values of K approach each other, overlap, and 
finally form two separate heads. The origin of the 
1—1 band is displaced toward the red from the 
0-0 origin. The Condon parabola should be 
practically a straight line, since as will be seen 
later the B’ and B” values are nearly equal. No 
other bands were observed. The wave numbers of 
lines for the 0-0 and 1—1 vibrational transitions 
are given in the Tables I and II. 

Because the first P and R lines were diffuse or 
too faint, the location of vp was uncertain. To find 
its position, the number of missing lines was 
determined by linear interpolation, and an arbi- 
trary numbering was made. The actual spin 
doublet splitting, which in terms of the true K” 
values (which are to be determined) must be 





Fic. 2. The *2->*2(0-—-0) 3840A transition of SiO*. The spin doubling and the doublet heads are clearly visible. At the 
band origin the spin doublets of the (1->1) transition can be observed but are of much weaker intensity. 
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TaBLe I. The system *2-*2 of SiO* (cm~). 








(0-0) vibrational transition 


(1-1) vibrational transition 





kK” Pa Po Ra Ro ) al P. Ps Ra Rs 
3 26010.75 11 26007.15 
4 26009.20 26022.25 12 26008.32 
5 26007.79 26023.59 13 26009.58 
6 26006.29 26024.98 14 26010.86 
7 26004.77 26026.32 15 26012.00 
8 26003.20 26027 .64 16 26012.99 26013.30 
9 26001.68 26028.90 17 26014.14 26014.50 
10 26000.03 26030.19 18 26015.19 26015.63 
19 26016.34 26016.72 
11 25998.29 26031.43 20 26017.46 26017.90 25960.06 25960.51 
12 25996.83 26032.69 
13 25995.12 26031.95 26033.82 | 21 26018.54 26018.95 — 25958.40 
14 25993.51 35.18 35.11 | 22 26019.52 26019.95 25956.28 25956.72 
15 25991.83 36.41 36.28 | 23 26020.55 = 25954.55 25954.99 
16 25990.29 25990.16 — 37.56 | 24 26021.69 26022.16 25952.80 — 
17 88.57 88.40 38.78 38.68 | 25 26022.73 26023.14 25951.02 — 
18 86.90 86.74 39.98 39.88 | 26 _ _— 
19 85.22 85.04 41.14 41.00 | 27 — — 
20 83.51 83.34 42.23 42.13 | 28 — 25946.09 
29 25943.77 25944.35 
21 81.75 81.58 43.39 43.19 | 39 = 25942.54 
22 79.82 79.80 44.50 44.34 
23 78.24 78.05 45.57 45.41 | 34 25940.06 25940.69 
24 76.46 76.24 46.71 46.47 
25 74.59 74.39 47.72 47.49 
26 72.84 72.58 48.70 48.48 
27 71.10 70.75 49.78 49.52 
28 69.27 68.90 50.82 50.47 
29 67.37 67.08 51.85 51.54 
30 65.46 65.15 52.78 52.51 
31 63.60 63.14 53.76 53.45 
32 61.66 61.37 54.71 54.41 
33 59.77 59.43 55.65 55.26 
34 57.88 57.48 56.50 56.21 
35 55.94 55.56 57.49 57.07 
36 53.98 53.53 58.30 57.90 
37 51.96 51.56 59.17 58.79 
38 49.94 49.60 59.98 59.66 
39 48.02 47.58 60.82 60.43 
40 46.02 45.61 61.64 61.25 
41 44.02 43.56 62.45 62.01 
42 41.98 41.52 63.26 62.86 
43 39.90 39.43 64.01 63.59 
44 37.77 36.78 64.78 64.34 
45 35.72 35.25 65.55 65.03 
46 33.66 33.14 66.24 65.80 
47 31.49 31.08 66.96 66.48 
48 29.38 28.83 67.79 67.09 
49 27.12 26.78 68.34 67.78 
50 25.04 24.64 68.99 68.44 
51 22.88 22.44 69.66 69.08 
52 70.26 69.66 
53° 70.85 70.26 
54 71.41 70.85 











represented by the equations‘ 


Avi2(R) = (7 — 7") K+3(37'— vy”) 


Avio(P) = (y' —y")K—-4(y' +7”) 


K=1, 2,3--- 


K=0,1,2--- in R, 


in P, 


4G. Herzberg, Molecular Spectra and Molecular Structure 
(1939), p. 274. 


was then plotted. The two lines obtained by 
plotting Avi.(P) against m=K and —Ap,.(R) 
against m= —(K-+1), form a continuous curve 
whose value at m=0 is —4(y'+~7’’) and whose 
value at m= —1 is —}(3y'’—y’’). The intercept 
of the curve with the m axis determined the true 
K numbering of the lines. This gave »» with a 
small uncertainty. The slope of the curve gave 
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TABLE II. Constants of the 3840A (22—22) band 
of SiO* (cm~). 











, B B” D a v(v’, v’’) 
0 O 0.7180 0.7253 2.25X10-& 2.11 X10-* 26,015.05 =0,0) 
1 1 0.704 0.712 "= 25,991.44 = 1,1) 








~ the value | -yo’—o’’| = 0.011 cm. From the sum 


of the intercepts with the two ordinate axes, the 
value yo’ = +0.012 cm was obtained, and the 
other possible K numbering gave the value 
vo’ = —0.0005 cm. Consideration of the signs 
of the intercepts and their relation to the 
signs of the y’s leads to the conclusion that if 
vo’ = —0.0005, then yo’ = —0.011, a possibility 
which was discarded on the grounds that the spin 
splitting coefficients of the *2 levels of the 
analogous molecule CO* are positive.' If one 
assumes that SiO* and CO? are in this respect 
similar, that is, that yo’’>0, then one finds 
that the conditions |yo’—~7yo0’|=0.010 and 
|yo’| =0.012 require yo’>0. It follows that 
vo’ =+0.012 cm, and yo’=+0.002 cm or 
vo’ = +0.022 cm—. Further information is needed 
to fix the value. The determination of y from com- 
bination differences, namely y=4}{A.7,(K+4) 
—A:T2(K —})} was found inadequate because of 
the smallness of the y’s compared with experi- 
mental errors. 
The quantities 


[Ri(K+}4) —Pi(K+3) ]/(4K+2), 
[R2(K —43) —P2(K —}4) ]/(4K+2), 
[Ri(K —1+4) —P:(K+1+}) ]/(4K+2), 


5 L. H. Woods, Phys. Rev. 63, 431 (1943). 





and 
[R2(K—1—}3)—P2(K+1—})]/(4K+2) 


were plotted against K(K+1). From the slopes 
of the lines D’ and D” were determined, and from 
the intercepts at K=0, B’ and B” were de- 
termined.* The very approximate values wo’ = 811 
cm and wo’’ = 851 cm were calculated from the 
relation w*=4B*/D. Although this relation is an 
exact one, the values of D are not accurate 
enough to give good values of w. 

This same procedure was followed to determine 
the constants of the 1-1 transition (Table II). 
However, not enough lines could be measured to 
determine the D values because the transition is 
severely overlapped by the 0-0 band. The 
quantity vo,o—v1,1=23.61 cm agrees with the 
calculated difference of the values of w in order of 
magnitude. 

The B values and the doublet structure of the 
lines show that the band at 3840A is emitted by 
the ionized molecule SiO*. 

If the SiO—SiO* energy level diagrams are 
roughly analogous to those of CO—CO* the 
existence of a *II—* band of SiO* lying in the 
infra-red can be expected, probably lying at 
wave-lengths greater than 8000A. 

The writer is very grateful to Dr. S. Mrozowski 
who suggested this problem, for his kindness, 
pertinent suggestions, and frequent help, and to 
Professor R. S. Mulliken for his unfailing interest 
and valuable advice. 


6 See reference 4, p. 199. 
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Note on Spin Doubling in the Doublet Sigma States of CO* 


L. H. Woops 
Ryerson Physical Laboratory, University of Chicago, Chicago, Illinois 
(Received March 9, 1943) 


On the photographic plates of the CO2 spectrum made by S. Mrozowski from the 30-foot, 
30,000 lines per inch, Chicago grating, were found pictures of the first negative bands and the 
comet tail bands of CO* in which the *2 spin splittings were resolved. From the *2-+* transi- 
tion, the absolute differences of the spin doubling coefficients were measured. The absolute 
values of yo” and y;” for the lowest 22 state were measured from the *II-—>*= transition, and their 
algebraic signs determined by consideration of relative intensities of branches. Thus the signs 
and values of the spin doubling coefficients of *’ for v’ =0 to 4, and of #2” for v’’ =0 to 7 have 


been determined. 





INTRODUCTION 


HE first negative bands of CO+(?2—*2) are 
well known from the early days of spectro- 
scopy. Their extensive history, which is well 
related by Kayser,' includes no occasion when 
the spin doublets of these bands were photo- 
graphed in such a way that they showed meas- 
urable resolution. Recently in his comprehensive 
study of CO, on the 30-foot grating at the 
Ryerson Laboratory, S. Mrozowski has obtained 
photographs of the first negative bands in the 
second and third orders, and of the comet tail 
bands in the second order, in which the 2 spin 
splittings were resolved. These plates provided 
an opportunity to amplify the known experi- 
mental data on spin splitting with which to 
compare predictions of the theory developed by 
Van Vleck.? Such comparisons have already been 
made for several molecules by Mulliken and 
Christy* who were able to draw conclusions as 
to the nature of the electronic states involved. 
In the present investigation the absolute value 
of the difference of the spin splitting coefficients 
for the upper and lower electronic states | y’—’’ | 
was obtained for the 0-1, 0-2, 0-3, 1-2, 
1-3, 1-4, 1-5, 2-44, 2-45, 2-46, 3-55, 3-6, 
3—7 vibrational bands of the *2->*> transition. 
From an examination of the 5-0, 6-0, 7-0, 
and 12-—+1 vibrational transitions of the *II-—*2 
bands, the absolute values of yo’ and 7," were de- 
rived. The Rz and ®Q2, Q2 and &P: branches of the 
5—0 and 6-0 bands were identified by con- 


1H. Kayser, Handbuch der Spectroscopie, Vol. 8. 

2]. H. Van Vleck, Phys. Rev. 33, 467 (1929). 

*R. S. Mulliken and A. Christy, Phys. Rev. 38, 87 
(1931). 
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sideration of their relative intensities. The sign 
of yo’ then was determined by the sign of 
Av=R2—®Q2=Q.—°P:. With this knowledge, 
the signs and values of the spin splitting coeffi- 
cients were determined for the vibrational levels 
0 to 4 of the upper ?2 state and for the levels 0 
to 7 of the ground 22 state. y’ was found to 
increase very slightly with increase in v’ while 
yy’ showed a rather marked increase with v”’. 
Moreover y’>vy”">0. A similar evaluation of 
change of spin splitting coefficient with vibra- 
tional quantum number was made for the 
2>— "> transition in AlO by M. K. Sen,‘ who 
found y’>vy’’>0, vy’ slowly increasing and y” 
slowly decreasing with increase in v. His paper 
contains a nice discussion of the theory involved 
in this analysis, but the accuracy of all constants 
given therein, especially the B values, is over- 
estimated. 


ANALYSIS 


The spin splitting in the *2—*2 bands was 
originally measured on the third-order plates 
but this set of data was abandoned in favor of 
that from the second-order plates, which because 
of greater intensity showed spin splitting to 
much higher K values, giving more experimental 
points and therefore more accurate graphs. For 
each band the separation of the two components 
of the spin doublets 


Avp(K) = Pi(K+}4) —P2(K —}) 
=(7'-—y")K—-4(y' +7") 
Ave(K) = Ri(K+4) — Ro(K —}) 
= (y'—")K+4(37'-7") 


4M. K. Sen, Ind. J. Phys. 11, 251 (1937). 
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TABLE I. Data from graphs of doublet splitting. 











Theoretical 
values for 
ly’ —y”| v v’ exp. y values pure precession 
0.0180 0 1 vo’ = +0.0102 — 0.026 
0.0159 0 2 vi’ = +0.0110 —0.025 
0.0153 0 3 v2’ = +0.0122 —0.025 
v3’ = +0.0136 —0.024 
0.0172 1 2 ys’ = +0.0145 —0.024 
0.0158 1 3 vs’ = +0.0141 — 0.024 
0.0150 1 + ve’ = +0.0144 — 0.023 
0.0159 1 5 y7"’=+0.0151 — 0.023 
0.0149 2 4 vo’ = +0.0290 +0.022 
0.0153 2 5 vi’ = +0.0294 +0.024 
0.0150 2 6 v2’ = +0.0294 +0.024 
0.0153 3 5 v3’ = +0.0297 +0.024 
0.0153 3 6 vs’ = +0.0297 +0.024 








vo’ = +0.0102 : y:"" = +0.0109. 


was plotted against m, where for the P branch 
Avp(K) was plotted against m=K, K having 
the values 1, 2, 3, - and for the R branch, 
—Avr(K) was plotted against m= —(K+1), K 
having the values 0, 1, 2, ---. The result was a 
single continuous line from whose slope the 
absolute value of (7’—-’’) was determined. 

The differences measured in the *II—>? bands, 


Av= R2(K) —*®Q2(K) = Q2(K) — °P2(K) 
= Q0(K) — °Ri(K) =P:(K)—"Q:(K), 


all of which should be given by the equation 
Av=T,"(K+4) —T2" (K—43) =v" (K+), were 
plotted against m=K-+4. The slopes of these 
lines determined the absolute values of yo’’ and 
yi’. The photographs of these bands were made 
under conditions more favorable to production 
of CO, spectra than CO* spectra, and unfortu- 
nately the images of the *II-—*2 bands were 
neither intense nor free of other spectra. Conse- 
quently the values of y’’ for higher values of v 
could not be obtained. 

The sign of yo’’ was found by consideration of 
the sign of Av in the 5-0 and 6-0 transitions 
of the *II-»*> bands. Since A/B)=80 for COt, 
the *II state approximates Hund’s case a. As 
the rotational quantum number increases, how- 
ever, the spin vector uncouples from the inter- 
nuclear axis, and there is a transition towards 
case 6. The main branches then fulfill the 


selection rule AK=+1,0, because for them 
AJ=AK. They remain strong, but the satellite 
branches suffer decreasing intensity, since for 
them AJ#AK.5 

If the lower *Z level is regular, i.e., y’’>0, 
F(K—4) lies below F(K+4), Qz lies farther to 
the ultraviolet than does °R:, and R;» lies on the 
ultraviolet side of "Q2 in the *II,;2—*2 band. 
In the band *II3;2—*, @R; and ?Q; must lie on 
the violet sides of Q; and P,, respectively, if 
y’>0. That is, in the *IIj;2—*2 band, the 
stronger line should lie toward the violet for 
the Qe, °Re and Re, *Qz branches, and in 
*II3/2—"2 the stronger line should lie toward 
the red for the @R:, Q: and ?Q;, P: branches. 
This is true both in the 5-0 and 6-0 transitions 
of the *II—*> band; therefore yo’’ is positive. 
It is assumed that y;"’ has the same sign as yo”. 
The 12-1 transition had very low intensity 
and was so overlapped by other bands that the 
sign of y,’’ could not be verified directly from 
the plates. The data from the graphs of the 
doublet splitting and the calculated values of 
the spin splitting coefficients are given in Table I. 

The experimental values of the y’s have been 
obtained by simple algebraic operations with 
the data. Theoretical values of y were calculated 
from Van Vleck’s formula 


y =2AB,l(1+1)/v(I, 2) 


by assuming a relation of pure precession between 
the *2 level considered and the “II level of the 
comet tail bands. A is negative for the comet 
tail *II level. »(II, =) is obtained by projecting 
up or down from the *II(v’) level to ?Z(v’’) level, 
and is positive or negative according as the 7II 
level lies above or below the 22 level.? Hence y” 
is negative and y’ is positive in the theoretical 
case. 
The data indicate a fair degree of likeness to 
a state of pure precession for the upper *2 state 
of CO+, but very little for the lower *2 state. 
The writer is sincerely thankful to Dr. S. 
Mrozowski, and to Professor R. S. Mulliken for 
their guidance and advice. 


5 R. S. Mulliken, Rev. Mod. Phys. 3, 89 (1931). 
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The Retardation and Neutralization of Argon Ions in Helium 


W. J. Hamu 
Department of Physics, Washington University, St. Louis, Missouri 


(Received March 21, 1943) 


An apparatus is described for making direct measurements of the loss of forward velocity of 
argon ions in single encounters with helium atoms. Results are given in the form of collision 
cross sections for given percentage retardations, for ions of energies of 1000, 800, 600, 400, and 
200 ev. The cross sections are zero for energy losses greater than 33.1 percent and increase both 
with decreasing percentage loss and decreasing beam speed. A correction to the readings neces- 
sary to obtain the results proves to be of unusual interest because the correction is itself the 
cross section for charge exchange between the argon ions and the helium atoms. Quantitative 
values of the cross sections for this unusual charge exchange are thus obtained. The values run 
roughly from 3.5 to 4.5 10~'* cm* per molecule in the energy range studied, compared with 
28 X 107-"* cm? for A* ions in A and 63X10-"* cm? for N2* ions in Ne. All the results are in 
sufficiently precise form to invite detailed theoretical analysis. 





INTRODUCTION masses results in small angular scattering. One 


NVESTIGATIONS of the collision of ions 0 the investigators of the Dempster group, J.S. 
and atoms have led to a number of experi- Thompson,’ sent Cs* ions through helium and 


mental procedures for studying the processes of hydrogen in apparatus similar to that employed 
scattering, neutralization, and retardation, which by Mayer‘ in his classic studies of the absorption 
affect a beam of ions shot through a gas. Various Of electron beams in gases. He developed a 


experimenters! using the Dempster mass spectro- method ; for deriving the Crees sections for 
graph technique, and in the majority of cases retardation from the resulting data on multiple 


employing the alkali ions emitted by Kunsman collisions. The present paper describes ented 
filaments, have published data on the relative ™ent devised to measure retardation cross 
importance of these processes for different ion- Sections directly for At ions in helium gas under 
gas combinations. The experimental procedures conditions such that only single collisions occur 
were frequently such that quantitative specifi- 1 the test region to any appreciable extent. 
cation of cross sections for the various processes 
was not possible. This was particularly true 
when retarding collisions were present, wherein The circuit diagram of Fig. 1 illustrates the 
the projectile ions suffered loss of forward speed essential parts of the apparatus. “S’’ represents 
without undergoing appreciable deviations. the ion source, which is a modification of the 
Theoretical considerations by A. V. Hershey? “High Efficiency Ion Source’ described by 
indicated that results of the retardation of ions 
in passing through a gas would involve fewer 


APPARATUS AND EXPERIMENTAL PROCEDURE 











experimental difficulties and be easier to analyze og 
than angular scattering results. Retardation 4 
measurements are most readily made by using 





ions very much heavier than the gas atoms of 
the retarding gas as such a combination of 








1A. J. Dempster, Phil. Mag. 3, 115 (1927); R. B. y. 
Kennard, Phys. Rev. 31, 423 (1928); I. W. Cox, Phys. 
Rev. 34, 1426 (1929); K. H. Bracewell, Phys. Rev. 54, 
639 (1938). 

2 Private communication. The writer is pleased to 
acknowledge his indebtedness to Dr. Hershey for furnishing 
his memoranda on the subject of retardation measure- 
ments. Many of his suggestions are incorporated in the 3]. S. Thompson, Phys. Rev. 35, 1196 (1930). 
experimental procedure. *H. F. Mayer, Ann. d. Physik 64, 451 (1921). 

















Fic. 1. Schematic diagram of the apparatus. 
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Fic. 2. Details of retardation chamber electrode assembly. 


Finkelstein.’ The source consists essentially of a 
combined magnetic and electric electron trap 
which maintains a high electron density in the 
source and serves the double purpose of using 
the electrons in a highly efficient way for 
ionization and of neutralizing the ordinarily 
inhibitory space charge of the positive ions. The 
source as used in this experiment is capable of 
currents of 2 to 10 ma of A* ions and has the 
highly desirable characteristic of operating at 
extremely low gas pressures. The beam energy 
used ranged between 200 and 1000 ev. Under 
these conditions the beam was not too sharply 
defined but was of such great intensity that it 
could be collimated with slits and homogenized 
with an electrostatic velocity selector to give a 
sharp beam of 10-7 amp., highly suitable for the 
present experiment. 
“H’’ represents the velocity selector.* ‘‘R’’ is 
; the retardation chamber, details of which appear 
h in Fig. 2. The source and retardation electrodes 
were housed in flanged Pyrex pipe sections 
(diam. 3 in.) which were attached to a heavy 
} steel box of diamond-shaped cross section con- 
taining the 127° selector. This steel box shielded 
the selector from the axial magnetic field used 
| with the source and served to connect the entiré 
i system directly onto the mouth of a three-stage 
fractionating diffusion pump’ of 80 liters per 
sec. speed. Rubber gasket technique was used 
throughout in making vacuum seals. Argon of 





5 A. T. Finkelstein, Rev. Sci. Inst. 11, 94 (1940). 

*A. L. Hughes and V. Rojansky, Phys. Rev. 34, 291 
| (1929). 
7™L. Malter and N. Marcuvitz, Rev. Sci. Inst. 9, 92 
(1938). 
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99.6 percent purity was fed into “S’’ through a 
variable leak,’ and helium of 98.2 percent purity 
into ‘‘R.”’ The gas pressures were measured by a 
D-79510 ionization manometer connected to a 
control and measuring unit similar to that 
described by Bowie.’ 

Power (a.c.) was supplied to the apparatus 
through a 1 kva Sola constant voltage trans- 
former, and direct current for the source and 
source magnet coils was supplied by a generator 
maintained at constant voltage by a General 
Electric type CD stabilizer. General Radio 
“Variacs’’ were used for adjusting accelerating 
voltages and source-cathode heating current. 
A shelf of ‘“B” batteries suspended by long 
insulating strips of ‘Amphenol A,” a polystyrene 
derivative, supplied the retardation potentials. 
These batteries, the switches, and the associated 
wiring which carried the positive ion currents 
were well shielded and were supported where 
necessary with ‘‘Amphenol A.” 

The retardation electrode assembly shown in 
Fig. 2 consisted of a can, a grid, and a plate. 
The grid and plate were moved as a unit by the 
sliding motion of the supporting rod through a 
double Wilson seal."° The region between the 
two seals was kept filled with helium at atmos- 
pheric pressure. For most measurements the 
can and grid were grounded so that the space 
between S; and G was practically field-free. The 
circuit connections were so arranged that At ions 
of a chosen energy, Vi, would be shot into He 
atoms in the field-free space where they would 
suffer their collisions along a path whose length 
could be varied at will, after which they would 
enter the analyzing region between the grid and 
the plate. The plate was held at a positive 
potential of V2 volts with respect to ground so 
that all ions which had energy less than V2 
electron volts would be unable to reach the 
plate. Throwing switch 5S, (Fig. 1) upward 
connected the entire retardation electrode as- 
sembly as a Faraday cage by means of which the 
total ion current J» in the beam was measured. 

Ideally, if the interaction of an A+ ions with 
a helium atom caused the former to be so 
retarded that its residual energy was less than 


8 R. D. Fowler, Rev. Sci. Inst. 6, 26 (1935). 


*R. M. Bowie, Rev. Sci. Inst. 11, 265 (1940). 
10R. R. Wilson, Rev. Sci. Inst. 12, 91 (1941). 





lil 


el 


be 


th 


Or 





oro Pp PS B 


O VS 


2 Gl bm oe 








ARGON IONS 


V2 electron volts, the retarding field would 
prevent it from reaching the plate. When J,, 
the plate current, is plotted as a function of V2, 
any ordinate of the curve is proportional to the 
number of ions retaining energies equal to, or 
greater than, the corresponding voltage V2. 
Energy and momentum considerations show 
that for the case of an argon ion colliding with 
a helium atom the minimum velocity retained by 
the ion is the fraction (M—m)/(M+m) =0.8179 
of its initial velocity, M and m being the masses 
of the ion and atom, respectively. Hence the 
greatest percentage loss of energy that can occur 
at such an impact is 33.1 percent, this maximum 
loss corresponding to a perfectly elastic, head-on 
collision. Any possible inelastic phenomena such 
as excitation or ionization would have the effect 
of decreasing this maximum possible percentage 
loss. The maximum deflection that the ion can 
experience is given by: @max=sin~! (m/M), or 
5° 46’. The geometry of the electrodes was such 
that for the longest path lengths used (13 cm) 
all ions so deflected could reach the plate if they 
were not rejected by the analyzing field. 

The graph of J, against V2 as it would appear 
under ideal experimental conditions is shown in 
Fig. 3 as a dashed curve. A typical actual 
experimental result is shown as a solid curve. 
The difference between the two curves may be 
analyzed as follows: The negative value of J, 
at values of V22 V; could only result from the 
emission of secondary electrons from the grid by 
impact of the positive ions and by fast neutral 
atoms which might be present as the result of 
the Kallman-Rosen ‘‘Umladung”’ or charge ex- 
change process. A choice of 40X40 mesh grid 
of 0.007-in. diameter nickel wire was made fol- 
lowing preliminary studies on field penetration. 
At the large values of V2, all ions of the beam 
eventually strike the grid because of its rela- 
tively negative potential, the ions passing be- 
tween the meshes on the first passage being 
deflected back to it. All secondary electrons 
liberated are collected by the plate. Hence the 
plate current is solely due to the secondary 
electrons and is just equal to the correction to 
be made. At the other extreme of voltage when 
V2 approaches zero, the ions which pass through 
the grid are not deflected back to it, so that 
only a fraction, 48/73, as many secondary 
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electrons are liberated, the fraction being derived 
from the effective grid open-area to metal-area 
ratio. For intermediate V2 values, the correction 
was taken as a linear function of V2, varying 
between the secondary electron current observed 
at high V2 and approximately two-thirds of this 
value. An assumption made in developing this 
correction procedure was that neutrals and ions 
of the same velocity were equally efficient in 
ejecting secondaries. Other studies" support this 
assumption. 

The generous slit dimensions adopted to 
permit use of a 10~-’-ampere ribbon of ions and 
galvanometer measuring technique allowed some 
leakage of source gas into the retardation 
chamber despite the great speed of the pump. 
The effects of this argon impurity in the target 
gas (He) were eliminated by taking double sets 
of readings with and without helium present, 
which yielded corresponding ‘‘helium data’’ and 
“foreign gas data.’ The effect of a slight amount 
of helium leakage into the source chamber with 
the consequent generation of a He* ion beam 
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Fic. 3. Graph of the current reaching the plate P, as a 
function of the retarding voltage V2. The dashed curve is 
an ideal one; the solid curve is experimental (helium 
pressure, 1.07 10-* mm Hg; 430 v At ions; for J», one 
division = 2.4 10-* amp.). The difference is attributed to 
the emission of secondary electrons from the grid G. 


was considered negligible because of the rela- 
tively low ionization efficiency of He by electrons 
of optimum energy for ionizing argon.” For a 
similar reason the ion beam must have con- 
tained a negligible percentage of A** ions." No 
spectrographic investigation of beam quality 
was attempted. 
A. Rostagni, Zeits. f. Physik 88, 69 (1934). 


2 P. T. Smith, Phys. Rev. 36, 1293 (1930). 
1 W. Bleakney, Phys. Rev. 36, 1303 (1930). 
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The following procedure was used in taking 
readings: Ions of energy V; were shot into the 
helium-free retardation chamber, and the re- 
tarding voltage between grid and plate was made 
to have values such that energy losses of p 
percent or greater were studied, where p was 
successively 0, 5, 10, 15, 20, 33, and also several 
values between 33 and 100. For each V2 setting, 
values of the plate current were recorded as the 
grid-plate assembly was moved in the retardation 


TABLE I. Experimental values of o, the cross section 
for specified percentage loss of energy, in units of 107!* cm? 
per molecule. . 











Energy of ion 1000 800 600 400 200 ev 
Energy loss r 
(p% or more) 
33 3.41 3.46 4.02 4.32 4.93 
25 3.599 3.747 4.411 4.900 5.635 
20 3.996 4.204 4.902 5.770 7.830 
15 4.42 4.72 5.49 6.86 _- 
10 5.08 5.57 6.53 8.67 — 
5 — 6.94 —_— _- — 








tube so as to alter the path length from 4 to 
13 cm. The current entering the can was main- 
tained at a constant value for a given ion 
velocity by frequent measurement of its value 
as a Faraday cage current and by proper 
adjustment of the source-magnet current so as 
to increase or decrease its focusing action on the 
beam. Similar experiments were next performed 
with helium in the retardation chamber, the 
other variables being set at the exact values 
obtained when no helium was present. 

The data obtained were analyzed in the 
following manner: Let Jo’ represent the value of 
the Faraday cage current, the total entering 
beam strength where path length is zero, and 
let J,’ represent the plate current in the ‘foreign 
gas run” when no helium has been admitted to 
the retardation chamber. Let J) and J, represent 
similar quantities in the experiments with 
helium present. Then, if sufficiently low pressure 
exists so that the use of single collision equations 
is justified, the foreign gas effect is given by: 
I,' =I)’ exp (— 8x), where 8 is the absorption 
coefficient representing chiefly the loss of At 
ions in argon by exchange-neutralization and 
also by scattering. When helium is admitted to 
the tube, the law of the beam-intensity decrease 
becomes: J,= J») exp —(8+ap)x, where a is the 
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cross section (Nic) of the atoms in 1 cm? of the 
helium at 1 mm Hg pressure for the ions of 
velocity determined by Vi, and for a percentage 
energy loss determined by V2. Taking the 
logarithms of these equations and subtracting 
we find: 


(log I,—log J,’)+constant = — apx= — Nipxe. 


Hence the curves obtained by plotting (log J, 
—log J,’) as ordinate against x as abscissa should 
be straight lines and possess a slope of value 
—ap, or —Nipo. 


RESULTS 


Argon ions of energies 1000, 800, 600, 400, 
and 200 ev were shot through helium gas 
maintained at a pressure of 1.07X10-* mm Hg. 
Observed values of Jo and J, for different path 
lengths and different V2 were tabulated, and 
the values of J, corrected for secondary electrons 
were calculated. A preliminary run was always 
made to determine the corresponding J,’ values. 
Curves of (log J,’—log J,) versus the path length 
x were then plotted. From such plots, the slopes 
which are proportional to the cross sections were 
derived. Table I gives the total cross section 
values derived from such graphs for the energy 
and loss criterion combinations of chief interest 
in this study. Cross sections remained constant 
when V,— V2 exceeded 33.1 percent. 


INTERPRETATION OF DATA 


At a single impact, an argon ion colliding with 
a helium atom cannot lose more than 33.1 


‘TABLE IT. Values of o,, the neutralization cross section in 
units of 10-'* cm* per molecule, for At ions in He. 








200 ev 
4.93 


600 
4.02 


400 
4.32 


800 
3.46 


Energy of ion 1000 


Cu 3.41 











percent of its energy. However, in no case did 
the slope of the (log J,—log J,’) vs. x graph, and 
hence o, assume a zero value when V2 was set 
for energy losses greater than 33.1 percent. The 
procedure for taking the readings eliminated 
most of the normal causes for such an apparent 
discrepancy in the results. Thus a_ possible 
divergence of the beam was doubly eliminated, 
first, because with a divergent beam a graph of 
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(log I,—log J,’) against x would not have been 
a straight line, which it was, and second, because 
the procedure of making the correction for 
foreign gas also eliminated geometric errors. 
Multiple collisions were eliminated as a pos- 
sible cause of the apparent error because the 
(log I,—log I,’) versus x curves would not have 
been straight lines if multiple collisions had 
occurred to a measurable extent. Foreign gases 
were adequately handled by the correction 
procedure used. Neither foreign gases nor foreign 
ions would be likely to cause an effect like the 
one observed, namely, a constant correction 
independent of the retarding voltage. Finally, 
it is repeated that no inelastic collisions in the 
gas can account for the observed effect. 

Only one possible explanation for the observed 
weakening of the beam as a function of path 
length but not of retarding voltage seems to 
remain. The effect will appear if A+ ions are 
removed from the beam by charge exchange 
with the helium atoms, the so-called Kallman- 
Rosen ‘‘Umladung” process. Accepting this 
interpretation, the results of the present experi- 
ment provide an excellent quantitative measure 
of the cross section for such charge exchange. 
Thus the o@ calculated above is composite: It 
should be written ¢=0o,+¢,, where o,, the cross 
section for exchange neutralization, is the value 
of o for the 33.1 percent loss criterion, and a,, 
the retardation cross section, is the difference 
remaining when go, is subtracted from the 
remaining table entries. 

The o, values are listed in Table II. These 
values may be compared with those of 63 x 10-"* 
cm? for N.* in nitrogen and 28X10~'* cm? for 
A* in argon. 

The occurrence of charge exchange between 
an ion of argon and an atom with a very different 

TABLE III. Values of o¢, for A* ions in He, the cross 


section for specific percentage loss of energy, in units of 
10-'6 cm? per molecule. 








Energy loss 
in percent 
33 
25 
20 
15 
10 
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Fic. 4. Collision cross section for energy loss by re- 
tardation plotted against percentage energy loss for A* 
ions of various speeds passing through He. 


ionization potential, like helium, has been the 
subject of some debate. Such measurements as 
existed indicated that the cross section must 
certainly be very small. In collisions between 
ions and atoms of different elements, two 
possibilities arise: The ionization potential of 
the projectile may exceed that of the target 
atom, or vice versa. Quantum mechanical con- 
siderations by Morse and Stueckelberg™ lead to 
the prediction that there is a certain likelihood 
of electron transfer in both cases which reaches a 
maximum when the ionization potentials of the 
two atoms are the same. Their results also 
indicate that the probability of transfer of an 
electron is an even function of the difference of 
ionization potentials provided the kinetic energy 
of the ion relative to the atom is large, the same 
therefore for an At ion meeting a He atom as for 
a Het ion meeting an A atom. The present 
experiment provides a most promising means for 
investigating this question further. 

The o, values, the cross sections for energy 
losses at single encounters, are listed in Table III 
and plotted versus percentage energy loss in 
Fig. 4. The ordinate of any curve gives the 





4 P. Morse and E. Stueckelberg, Ann. d. Physik 9, 589 
(1931). 
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cross section within which the ion must pass if 
it is to incur at least the percentage energy loss 
indicated as abscissa. It is believed that the 
curves will lend themselves in a very satisfactory 
manner to analysis of the forces existing between 
the particles. 

The author wishes to acknowledge gratefully 


the guidance and valuable suggestions of Doctor 
Robert N. Varney who directed this problem by 
correspondence from the Naval Proving Ground 
at Dahlgren, Virginia; also the courtesies ex- 
tended by Professor A. L. Hughes, Professor G. 
E. M. Jauncey, and Doctor Joseph Keller of 
this department. 
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Diffusion Rates of Carbon in Iron-Molybdenum and Iron-Tungsten Alloys 


R. SMOLUCHOWSKI 
Research Laboratory, General Electric Company, Schenectady, New York 
(Received March 24, 1943) 


Measurements are reported of the influence of tungsten and molybdenum on the diffusion 
rate of carbon in face-centered iron at 1000°C. Both of these elements slow down the diffusion of 
carbon, the influence of tungsten being more than twice that of molybdenum. These results are 
compared with previous measurements on iron-cobalt alloys and other known data. No con- 
nection seems to exist between the variation in lattice parameter of face-centered iron and the 


variation in diffusion rate of carbon. 


N a previous paper! a strong accelerating 

influence of cobalt on the diffusion rate of 
carbon in face-centered iron was reported.* In 
continuation of this study, experiments were 
made with iron-molybdenum and iron-tungsten 
alloys. 

There are relatively few binary systems in 
which measurements of this type can be made 
conveniently. The main condition is the existence 
of a sufficiently wide range of homogeneous face- 
centered solid solution. It should be possible to 
vary the concentration of carbon and of the 
alloying element, within relatively wide limits, 
to suit the accuracy of the measurements. Both 
in tungsten-iron and in molybdenum-iron alloys 
no single phase region with face-centered lattice 
exists beyond a few percent of the added element. 
Carbon, however, increases this range somewhat. 
Accordingly, experiments were made with the 
following alloys; the tungsten pair: 1.1 atomic 
percent (about 3.5 wt. percent) tungsten with 


1 R. Smoluchowski, Phys. Rev. 62, 539 (1942). 

* Note added in proof: Mr. Malcolm F. Hawkes and 
Dr. Robert F. Mehl of the Carnegie Institute of Technology 
kindly informed me that their experiments do not show 
this strong accelerating influence of cobalt. This point 
needs thus further experimental clarification in particular 
in view of reference 6 quoted in reference 1 of this paper. 


2.6 atomic percent (0.58 wt. percent) carbon and 
with 0.22 atomic percent (0.05 wt. percent) 
carbon; the molybdenum pair: 1.1 atomic per- 
cent (about 1.9 wt. percent) molybdenum with 
3.1 atomic percent (0.7 wt. percent) carbon and 
without carbon, and finally a pair of alloys with 
the same molybdenum content but with 4.2 
atomic percent (about 0.9 wt. percent) carbon, 
and without carbon. The first two groups of 
alloys, the tungsten alloys and the molybdenum 
alloys with lower carbon content, fall within the 
region of homogeneous face-centered lattice at 
1000°C.2 The third group, the molybdenum 
alloys with higher carbon, fall into the region 
where at 1000°C some of the carbides do not 
dissolve. These alloys were included in the 
experiments in order to see to what extent the 
presence of carbides affects the observed diffusion. 

The experimental method was essentially the 
same as that used in the previous investigation. 
The welded samples were heated for three days 
at 1000°C and the concentration of carbon in 
layers parallel to the interface was determined. 


2S. Takeda, Tech. Rep. Sendai 9, 483 (1930); 9, 627 
(1930); and 10, 42 (1931). E. C. Bain, The Alloying 
Elements in Steel (Am. Soc. Metals, Cleveland, Ohio, 
1939). 
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An interesting phenomenon was observed in the 
case of tungsten alloys. In order to prevent 
decarburization, all samples were enclosed in 
quartz tubes filled with argon. In each tube were 
two samples; one with iron containing molybde- 
num or tungsten and another with pure iron, in 
order to be able to compare the diffusion rates 
in iron and in iron alloy under identical condi- 
tions. Determination of the carbon concentration 
in the first heat-treated tungsten-iron samples 
showed that the carbon-rich side was somewhat 
decarburized and the carbon-poor side was 
carburized on the surface. There was a radial 
gradient of carbon concentration. No such effect 
was observed on the pure iron sample which 
was enclosed in the same quartz tube and thus 
underwent the heat treatment simultaneously 
with the iron-tungsten sample. An explanation 
of this effect, on the assumption that the argon 
used was impure, is difficult to accept because 
of the absence of this decarburization and 
carburization in the pure iron sample. Such an 
effect is, of course, very undesirable because it 
affects the concentration gradient due to volume 
diffusion and thus falsifies the results. In order 
to prevent it, the samples were covered electro- 
lytically with a 0.3-mm thick layer of copper 
before heat treatment. This protective measure 
was applied to all samples, and it eliminated the 
occurrence of these surface phenomena. The 
computation of the diffusion coefficients was 
made by the usual Grube method. Because of 
the presence of the alloying elements on both 
sides of the interface, the computations were 
simple, and the use of the corrective formulae 
derived in the previous paper was not necessary. 

The numerical results at various carbon 
concentrations are given in Table I. 

The decelerating influence of both tungsten 
and molybdenum on the diffusion rate of carbon 


TABLE I. Diffusion coefficients of carbon in iron- 
molybdenum and in iron-tungsten alloys at 1000°C. 
DX 10? cm? sec.~!. 











Solvent 
Atomic 1.1 at. 1.1 at. 
percent Pure percent percent 
carbon Fe Mo Ww 
0.6 2.9 2.5 1.1 
1.1 3.2 2.6 1.2 
2.2 3.7 2.8 1.25 
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Fic. 1. Diffusion coefficients of carbon in various iron 
alloys at 1000°C. (Data for Mn and Ni were obtained 
from reference 3. Crosses indicate values for Mo based on 
reference 4.) 


is clearly seen. For the same atomic percentage 
of added element, tungsten is much more 
effective than molybdenum. In Fig. 1 are plotted 
the diffusion coefficients of carbon in various 
iron alloys at 1000°C. Former measurements on 
iron-cobalt alloys! are included, also Wells and 
Mehl’s data on iron-nickel alloys and iron- 
manganese alloys’ and the recent measurements 
of Ham, Parke, and Herzig on iron-molybdenum 
alloys.‘ The latter were made with 0.47 atomic 
percent (0.80 wt. percent) molybdenum at 
various temperatures. The corresponding points 
in Fig. 1 were obtained by interpolation from 
Table III in Ham, Parke, and Herzig’s paper at 
1.1 and 2.2 atomic percent carbon. We see that 
the agreement is satisfactory, the influence of 
molybdenum being somewhat stronger than that 
reported here. The difference easily may be 
due to experimental errors and interpolation 
inaccuracies. 

The diffusion coefficients given in Table I 


3C. Wells and R. F. Mehl, Trans. A. I. M. E. 140, 279 
(1940). 

‘J. L. Ham, R. M. Parke, and A. J. Herzig, Am. Soc. 
Metals, Preprint, 1942. 





Fic. 2. Activation ener- 
gies of the diffusion of 
carbon in various iron 
alloys. 
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and Fig. 1 were obtained with alloys which form 
solid solutions at 1000°C. As mentioned before, 
experiments were made also with another mo- 
Ivbdenum alloy which had a little more carbon 
than the solubility limit at that temperature. 
The diffusion coefficients obtained from this 
alloy were within 15 percent of those given in 
Table I. This shows only that the few carbides 
present were quickly dissolved when the concen- 
tration of the surrounding solid solution de- 
creased and thus they did not affect appreciably 
the total diffusion rate. In general, one would 
expect deviations from the simple concentration 
distribution law and other complicating effects. 

In the previous paper, it was shown that the 
Dushman-Langmuir equation gives values of the 
activation energy Q in good agreement with 
experimental data in the case of iron-cobalt 
alloys. The same formula was applied here, and 
we assumed as before that the ‘‘free path” 6 
was 1.8A. We shall not repeat the details and 
discussion of the limitations of this procedure. 
The results for 1.1 atomic percent carbon are 
shown in Fig. 2 where the values for various 
iron alloys, including those investigated earlier, 
are plotted against the composition. A strong 
increase of the activation energy due to the 
presence of tungsten and molybdenum is ap- 
parent. It should be stressed that it is perfectly 
justifiable not to take into account the variation 
in 6 with composition. A change in the lattice 
constant of one percent changes the value of Q 
by about 50 calories, which is, of course, entirely 
negligible as compared with the effects here 


Mi observed. 
| It seemed interesting to see if there is any 


R. SMOLUCHOWSKI 


correlation between the changes in the lattice of 
iron brought about by the dissolved elements 
and the change of the rate of diffusion of carbon. 
Practically no data are available on the lattice 
constants in the face-centered lattice range of 
the iron-tungsten and iron-molybdenum alloys. 
In order to obtain a rough estimate of the 
influence of these two elements on the lattice of 
face-centered iron, measurements were made of 
the density and thermal expansion in the 
proximity of 1000°C. It turned out that both 
tungsten and molybdenum expand the iron 
lattice at a rate which is close to that computed 
from atomic radii, that is, about 2 percent per 
10 percent of added element. Let us compare 
the change in the diffusion rate of carbon in 
iron due to the presence of other metallic 
elements as pictured in Fig. 1 on the one hand, 
and the change in lattice constants on the other. 
It turns out that there is no simple relation: 
cobalt, which leaves the iron lattice constant 
practically unaltered,* accelerates the diffusion 
markedly. On the other hand, carbon, which 
accelerates its own diffusion almost as effectively 
as cobalt, increases the lattice constant of iron 
rapidly, because of its interstitial location. 
Molybdenum and tungsten, which both have 
practically identical atomic radii, decrease this 
rate to a very different extent. The influence of 
tungsten is more than twice as strong as that 
of molybdenum. Nickel has practically no in- 
fluence on the rate of diffusion, but it decreases 
the lattice constant of the face-centered iron.® 
Manganese, like nickel, leaves the diffusion rate 
practically unchanged, but increases the lattice 
constant of face-centered iron.’ In the previous 
paper it was stated that the simple model of 
carbon atoms diffusing through the crystal 
lattice of iron, in which only the geometrical 
factors play a role, is untenable and that only 
a study of the influence of the alloying elements 
on the electronic interactions between the iron 
lattice and the carbon atoms can correlate the 
observed facts. The results here reported add 
support to this conclusion. 


5 A. Osawa, Sci. Rep. Tohoku 19, 109 (1930). 

*E. A. Owen, E. L. Yates, and A. H. Sully, Proc. 
Phys. Soc. 49, 315 (1937): H. Esser and G. Miiller, Arch. 
f. d. Eisenh. 7, 265 (1933-1934). 

7E. Ohman, Zeits. f. physik. Chemie B8, 81 (1930). 
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Investigations of Ferromagnetic Impurities. II 


F. W. Constant, R. E. Farres, AND H. E. LENANDER 
Duke University, Durham, North Carolina 


(Received February 23, 1943) 


Continuing previous studies on minute ferromagnetic impurities in ordinary metals, an 
attempt was made to render such impurities non-magnetic through heat treatment. In the case 
of copper and brass it was found that heating to 900°C and 750°C, respectively, in hydrogen for 
two hours or more, followed by quenching, would eliminate the ferromagnetism. This method 
was not successful with silver. These results are interpreted in terms of the crystal structures 
of the metals and of iron, which is believed to be the chief magnetic impurity. In the case of 
copper, reheating to 900°C in an atmosphere containing oxygen caused the ferromagnetism to 
return strongly. This is explained as due to the diffusion of the oxygen into the copper at high 
temperatures. The effect of heat treatment at other intermediate temperatures on the magnetic 
hardness and hysteresis of the impurity was also investigated. 





INTRODUCTION 


i an earlier paper! (Part I), results are given of 
measurements on the ferromagnetic impuri- 
ties found in various commercial materials such 
as copper, silver, brass, aluminum, etc. The 
method involved suspending a specimen by a 
quartz fiber in the magnetic field of a pair of 
Helmholtz coils. In this way it was possible to 
measure either the remanent magnetization of 
the ferromagnetic impurities left after a specimen 
had previously been exposed to the much stronger 
field of an electromagnet, or, by means of a 
second set of magnetizing coils at right angles to 
the Helmholtz coils, to measure a hysteresis loop 
for a given specimen. As a result of this work it 
was found that most commercial metals con- 
tained ferromagnetic impurities which could 
easily be detected magnetically although com- 
pared with the remanent magnetization of pure 
iron, the magnetization of these impurities could 
represent an amount of impurity far too small to 
be detected by chemical analysis. The hysteresis 
measurements showed that the impurities were 
magnetically ‘‘hard.”’ 

Since in some cases, such as galvanometer coils, 
standard inductances, etc., magnetic impurities 
could be objectionable, it seemed of interest to 
investigate whether such impurities could be 
rendered non-magnetic by heat treatment. It was 
also heped that in this manner more would be 


1F, W. Constant and J. M. Formwalt, Phys. Rev. 56, 
373 (1939). 


learned about the state of the magnetic impurity 
in the parent metal. For this purpose three 
common metals were chosen, namely copper, 
brass, and silver. The effect of heat treatment in 
various atmospheres on both the remanent 
magnetization and hysteresis was investigated. 


EXPERIMENTAL 


A given specimen was surface cleaned, mag- 
netized, and its remanent magnetization meas- 
ured in the usual way. The coercive force was 
also found by using increasing negative fields 
until the remanent magnetization was reduced to 
zero. The specimen was then placed in a pro- 
tecting box of molybdenum (which showed no 
ferromagnetism), and this in turn was introduced 
into a Leeds and Northrup ‘“Fyrestan”’ tube 
situated in a resistance furnace. The tube could 
be evacuated and various gases run through at 
one or two atmospheres pressure. The furnace 
was then heated up to the required temperature 
and held there for eight hours after which the 
tube was opened, and the box containing the 
specimen was removed quickly and quenched in 
oil. In a few cases the specimen was allowed to 
cool overnight in the furnace. Although heating 
to high temperatures would demagnetize a speci- 
men, it did not necessarily follow that the 
specimen was non-magnetic. Hence after heat 
treatment each specimen was again placed in a 
magnetizing field and the remanence, coercive 
force, or hysteresis measured. 
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Fic. 1. Ratio of remanence after heat treatment of 
copper in hydrogen to the original remanence expressed 
in percent. 


RESULTS 
1. Copper 


Samples were taken from No. 12 D.C.C. 
magnet wire. Asa test for uniformity two or three 
samples were run through together; the resulting 
agreement was always good, differences not ex- 
ceeding a few percent. The remanence was meas- 
ured both before and after heat treatment, the 
same magnetizing field being used each time. The 
new remanence was computed as a percent of the 
original, and this was plotted against the temper- 
ature of the 8-hour run. Figure 1 shows the 
results for heating in a hydrogen atmosphere, 
which was naturally tried first in order to keep 
the specimen “‘clean.’’ After a slight initial de- 
crease the curve rises to a peak at about 525°C 
for which the new remanence is some five and a 
half times as great as the original. Above 525°C 
the curve falls rapidly to zero at 875°C. Speci- 
mens heated to 900°C or higher for two hours or 
more were completely non-ferromagnetic when 
cooled. In the case of copper slow cooling gave 
practically the same results as quenching. 

In order to understand better the curve in 
Fig. 1, hysteresis loops were run on samples be- 
fore and after heating to 300°, 500°, and 700°C, 
respectively. These loops were obtained by meas- 
uring the period of the specimen when suspended 
parallel to the magnetizing field. This field ex- 
erted a restoring torque proportional to the 
magnetization which aided the restoring torque 
of the fiber and hence shortened the period of 
oscillation. Heating to 300°C produced only a 


slight change in the hysteresis loop. Figure 2 
shows loops before and after heating to 500°C 
and the heat treatment has evidently made the 
magnetism ‘‘harder,”’ increasing both the rema- 
nence and coercive force. This would largely ex- 
plain the maximum in the curve of Fig. 1 at about 
525°C. Heating to 700°C, on the other hand, 
decreased the remanence somewhat, the coercive 
force remaining about the same (55 oersteds). 
As a further check on the above results a speci- 
men was magnetized in a field of about 4500 
oersteds and then placed in a weak field in the 
opposite direction which could gradually be in- 
creased until the remanence disappeared. This 
gave approximately the coercive force for a 
larger loop between plus and minus 4500 oersteds. 
Here again heating to 500°C increased the 
coercive force from 215 to 320 oersteds whereas 
heating to 700°C only increased it slightly. 
Copper samples were next heat treated in a 
nitrogen atmosphere. The nitrogen used was 
ordinary tank nitrogen. The percent of the 
original remanence after various heatings is 
shown in Fig. 3. The results differ surprisingly 
from those shown in Fig. 1 for a hydrogen 
atmosphere. First, the strong peak at 525°C has 
disappeared, and secondly, around 900°C an in- 
crease in remanence has actually occurred. 
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Fic. 2. Hysteresis of copper before and after 
heating to 500°C. 
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Finally, temperatures of at least 1000°C were 
necessary to eliminate the ferromagnetism. Since 
air gave similar results to nitrogen, it was be- 
lieved likely that the tank nitrogen contained 
sufficient oxygen impurity to cause the effect. 
Attempts were made to purify the nitrogen by 
passing it over red hot copper and drying it and 
repeatedly evacuating and refilling the system. 
Asa result the 900°C peak was almost eliminated. 

Another interesting point was that a sample 
which had been rendered non-magnetic in hydro- 
gen at 900°C could be made more magnetic than 
originally by reheating in air or tank nitrogen at 
900°C. This was not a surface effect. This cycle 
could be repeated indefinitely unless tempera- 
tures exceeding 1020°C were used in which case 
the magnetism returned with more difficulty. 

Heating in hydrogen at 900°C or higher always 
eliminated the ferromagnetism. Although re- 
heating in air at 900°C would cause it to return, 
heating or annealing at lower temperatures such 
as 500°C, whether in air or hydrogen, caused no 
such return. The conclusion was that oxygen 
diffusing into the copper above about 850°C was 
responsible for the return of, or increase in, 
ferromagnetism between 850° and 950°C. 

Finally, following the work of Smith?* and 
Bitter,**® it was found that cold rolling non- 
magnetic specimens caused a return of the 
ferromagnetism. 
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Fic. 3. Ratio of remanence after heat treatment of 
copper in nitrogen to the original remanence expressed 
in percent. 


2. Brass 


The results of heat treating samples of a brass 
rod in hydrogen at various temperatures are 
shown in Fig. 4 in which, as in Fig. 1 for copper, 
the change in remanence is plotted. The brass 


2C. S. Smith, Phys. Rev. 57, 337 (1940). 

?C. S. Smith, Am. Inst. Min. Metal. Eng., Tech. Pub. 
No. 1394 (1942). 

4F. Bitter and A. Kaufmann, Phys. Rev. 56, 1044 
(1939). 

5F. Bitter, A. Kaufmann, C. Starr, and S. T. Pan, 
Phys. Rev. 60, 134 (1941). 
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specimens were about ten times as magnetic as 
the copper, but the curve is similar except that 
the peak occurs earlier and magnetism disappears 
when the heating is carried to about 750°C. 
Specimens from a different source of brass gave a 
similar curve except for a much higher peak (650 
percent) at 400°C. Above 750°C quenching was 
necessary if the last 5 percent of remanence were 
to be removed. 

Hysteresis loops were also taken for brass. 
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Fic. 4. Ratio of remanence after heat treatment of 
brass in hydrogen to the original remanence expressed in 
percent. 


However, it was impossible to obtain a loop of 
any size for untreated specimens because the 
magnetization was sufficiently ‘‘soft’’ to rotate 
with the rotation of the specimen or of the 
resultant field. But after heating to 400°C or 
500°C loops could be obtained, indicating, as for 
copper, an increase in magnetic hardness. 
Reheating non-magnetic brass specimens at 
intermediate temperatures did cause a return of 
the magnetism. Probably due to the greater 
amount of impurity in the brass, the solubility 
limit for the impurity became a factor and the 
non-magnetic state was then less stable at 
intermediate or room temperatures. 


3. Silver 


Similar heat treatments in hydrogen carried 
out on specimens of silver wire gave much less 
striking effects, as shown in Fig. 5, where the 
remanence is again plotted. Heating up to the 
melting point failed to reduce the ferromagnetic 
remanence by more than half, and this decrease 
was probably due to the magnetic impurity be- 
coming magnetically softer. 


DISCUSSION 


The ferromagnetic impurity was taken to be 
iron, perhaps picked up by the specimen when 
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4 
mt work and that of Smith? and Bitter* § that copper ‘du 
Pe containing up to more than one percent iron in all 
<<: ecueghalilltiaa solid solution at 900°C may be cooled to room co 

temperature retaining the iron in solution in a 
metastable, non-magnetic form. Annealing at 








100° 200° 300° 400" 800° 600° 700° 600° 900% 7 : 
Fic. 5. Ratio of remanence after heat treatment of intermediate aa tee com _ restare the 
silver in hydrogen to the original remanence expressed in Magnetism although it may precipitate the irom -_ 
percent. out in a non-magnetic state. Cold work, however, he 
: ; restores the magnetism, and Smith has showr i 
drawn through dies, or during cold rolling. The that there is at the same time a reappearance of; 
fact that the original remanence was partly a body-centered lattice, presumably through ~ 
retained if copper and brass specimens were modification of the iron. Hence the mechanism o! 
heated to temperatures up to 550°C and not rendering a specimen of copper non-magnetic is 
remagnetized after cooling indicated that the imply to heat it at a sufficient temperature for a 
Curie temperatures for the impurities in these  cuficient time to dissolve the iron impurity in the 
metals did not lie below 550°C. This made it very copper lattice. Subsequent cooling or annealing 
unlikely that nickel could be the chief ferro- — }¢Jow 850°C will not alone restore the magnetism. 
magnetic impurity. Further, contamination with The case of brass would seem similar to that of 
cobalt was much less likely than with iron. copper. The lattice constant of brass (face- , 
Finally, brass gave a positive chemical test for centered cubic) is almost exactly that of face-cen- 
iron. tered iron. The solubility of iron in brass takes 
Although the work of Tammann and Oelsen® place even more readily than in copper, a lower 
would indicate around one percent impurity when temperature for brass than iron being sufficient to 
900°C is required to render a copper-iron alloy dissolve a greater amount of impurity in a give» 
non-magnetic, it is very unlikely that in our case time. : 





= much was present. The 900°C is given as a safe The lack of complete disappearance of thé ; 

upper limit for rendering copper wire non- magnetism upon heating silver indicates lack of 
| magnetic. For oe thing, quenching at 900°C €X- complete solubility of the impurity in the silver. pe 
| posed the specimen to air (oxygen), which could Although silver is face-centered cubic in structure el 
| act on, it as mentioned above. Further, longer its lattice constant, 4.60A, differs considerably ar 
annealing at 800°C was found equally effective in fom that of y-iron. As even the liquid phases of wi 
| eliminating the magnetism. silver and iron do not mix little solubility can be pe 
The phase diagram for copper-iron alloys as expected for the solid phases. So 
; given by Hansen’ shows that at 900°C up to two Of final interest is the comparison of the effect in 
percent of iron will dissolve in copper as a solid of an oxygen and a hydrogen atmosphere during pe 
solution with the face-centered cubic structure pest treatment of copper at 850°C or higher. pi 
characteristic of copper. This is reasonable since goth gases have definite effects. However, con- os 
| the lattice constant of copper is 3.62A* and at siderably higher temperatures are required before t fe 
| high temperatures iron is also face-centered cubic oxygen can diffuse through copper than in the th 
in structure with a lattice constant of 3.60A. At case of hydrogen, hence the effect of oxygen is w 
} room temperature the solubility of iron in copper noticed only at 850°C or higher. The simplest he 
| is practically nil since iron is then normally body- explanation of the results given above would th 
:| centered cubic in structure. This latter structure seem to be that (1) the iron impurity exists sc 
| is associated with the ferromagnetic state of iron. partly as an oxide in the copper and this oxide is ki 
i However, it now seems evident from the above reduced by hydrogen giving an increase in ne 
i magnetism at around 550°C, (2) the iron dis- C 
Hf Chesie 196 287 (1930) Oelsen, Zeits. f. anorg. allgem. solves in the copper lattice at higher tempera- tk 
¥ 7M. Hansen, Der Aufbau der Zweistofleigierungen tures, (3) introduction of oxygen at higher SC 
F| OW oe Weebod, The Structure of Crystals (Little temperatures causes the dissolved iron to be . Sl 


and Ives, New York, 1924). reprecipitated out as an oxide, and (4) reintro- la 
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‘duction of hydrogen again reduces the oxide and __ perature to allow the magnetic impurity to dis- 
allows the iron once again to dissolve in the solve in the copper lattice. A temperature of 
copper. 900°C is safe for a two-hour treatment. This 
method was unsuccessful in the case of silver, but 
might work for certain other metals with which 
Commercial specimens of copper and brass iron can form solid solutions. The magnetism can 
may be rendered completely non-magnetic by be made to return in copper by reheating in air 
heating in hydrogen to a sufficiently high tem- above 850°C, or by annealing plus cold work. 


CONCLUSIONS 
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Magnetic Properties of Solid Solutions 
III. The Paramagnetic Alloys of Copper and Nickel 


. A. R. KAUFMANN AND C. STARR 


Massachusetts Institute of Technology, Cambridge, Massachusetts ‘ 
(Received December 29, 1942) 


The magnetization of the paramagnetic alloys of copper with nickel has been studied as a 
function of temperature and applied field. It is found that a Curie-Weiss equation cannot 
represent the results and an equation of the form x =a7+6+(c/T) is found to be satisfactory. 
Alloys in the range from 30 to 38 percent nickel showed a considerable field dependence of the 





susceptibility at low temperatures. 


HE magnetic properties of copper-nickel 
alloys have been investigated by many 
people but most of the work has been done at 
elevated temperatures and on those alloys which 
are definitely ferromagnetic. The present work 
was intended primarily to cover the low tem- 
perature properties of those alloys which are not 
ferromagnetic at room temperature, although the 
investigation was extended to elevated tem- 
peratures in order to get a more comprehensive 
picture of the magnetic behavior. In this alloy 
series there are compositions which are neither 
ferromagnetic nor truly paramagnetic even in 
the Curie-Weiss sense and it is this phenomenon 
which was looked for in the investigation. It was 
hoped that such a study would reveal not only 
the magnetic moment of nickel atoms in dilute 
solution in copper but also contribute to the 
knowledge of conditions which lead to ferromag- 
netism as the concentration of nickel is increased. 
Copper-nickel alloys are particularly suited for 
this purpose since they form a complete series of 
solid solutions in which the added nickel atoms 
simply replace the copper atoms on the crystal 
lattice of the latter. 


PREPARATION OF ALLOYS 


The materials used for the alloys were Mal- 
linkrodt’s reagent copper' and reagent nickel?* 
which were melted together in quantities of about 
40 grams. Melting was carried out in an alundum 
crucible and in a hydrogen atmosphere, the tem- 
perature of the furnace being about 1500°C. The 
alloys were chill cast in a vacuum directly from 
this melting temperature. These castings were 
given a homogenizing anneal at 900°C in hy- 
drogen for 36 hours and then were cold swagged 
with intermediate anneals into the form of 
7s-inch diameter Tods. Short pieces weighing 


TABLE I. Composition of alloys, weight percent. 














Alloy Ni Fe Alloy Ni Fe 

J27 —_ 0.00 J25 20.34 0.016 
J7 0.94 0.0023 j24 29.98 0.024 
J8 4.99 0.0046 J23 35.59 0.030 
j9 10.10 0.0094 j22 38.00 0.029 
326 10.99 0.012 y21 50.16 0.034 








1Sb —0.01 percent, As —0.001 percent, Fe —0.004 per- 
cent, Pb —0.01 percent, Sn —0.01 percent. 

2Co —0.01 percent, Fe —0.10 percent, Pb —0.01 per- 
cent. 
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TABLE II. Susceptibility per gram of alloy. 








Weight percent nickel 





T-°K o* 0.94 4.99 10.1 11.0 20.3 30.0 35.6 38.0 50.2 
14.0 — 0.034 — 0.023 0.125 0.581 1.27 4.377 14.8f 44.77 
20.4 — 0.032 0.104 0.362 0.483 1.10 3.62T 11.0 31.5t 
77.3 — 0.073 —0.052 0.052 0.216 0.274 0.67 1.81 4.01 7.70¢ Ferromag, 
295 — 0.0796 — 0.062 0.042 0.185 0.221 0.553 1.277 1.993 2.84 11.8 
523 — 0.0786 — 0.060 0.047 0.232 0.592 1.222 1.675 2.11 
773 — 0.0748 — 0.055 0.063 0.262 0.635 1.200 1.561 1.85 
1023 — 0.0693 — 0.045 0.079 0.300 0.687 1.213 1.515 1.74 
1273 — 0.0645 — 0.043 0.103 0.339 0.735 1.230 1.497 1.69 











* Copper with 0.0025 percent iron described in preceding paper, Phys. 


Rev. 60, 134 (1941). 


t+ Measurements showing field dependence. Initial susceptibility given. 


about 4 grams were cut from the rods and were 
used as the specimens for measurement. Im- 
mediately adjacent pieces were cut off for 
chemical analysis for iron and nickel, the results 
of which are given, in Table I. The specimens 
were given a final anneal in hydrogen and cooled 
in the furnace in order to remove the effects of 
cold work. 


‘ METHOD OF MEASUREMENT 


The same technique and the same magnet were 
used for these measurements as has been de- 
scribed in the preceding papers* of this series. 
For measurements below room temperature the 
specimens were suspended from the balance with 
a silk thread while a sling made of molybdenum 
foil was used for this purpose at elevated tem- 
peratures. No attempt was made to make a cor- 
rection for the force on the silk thread but a 
correction was made whenever the molybdenum 
sling was used. 


EXPERIMENTAL RESULTS 


The susceptibility of the alloys was determined 
by plotting the magnetization against the field 
and taking the slope of the straight line which 
is found between 10 and 30 kilogauss. This 
procedure automatically corrects for any traces 
of ferromagnetic particles since these always are 
practically saturated in fields of this magnitude. 
The ferromagnetic inclusions in all cases were 
almost negligible. 

Certain of the alloys, which were normally 
paramagnetic at room temperature, began to 
show progressively larger saturation effects in 


3F. Bitter and A. R. Kaufmann, Phys. Rev. 56, 1044 
(1939); Bitter, Kaufmann, Starr, and Pan, Phys. Rev. 60, 
134 (1941). 


high fields as the temperature was lowered. As 
with the copper-iron alloys, this effect was con- 
sidered to be a true property of the alloy and not 
just simply the result of ferromagnetic inclusions. 
In such cases the initial slope of the magnetiza- 
tion curves was taken as the initial susceptibility 
since the room temperature measurements 
showed an insufficient amount of ferrromagnetic 
particles to affect appreciably this quantity. An 
example of this type of behavior is shown in Fig. 
1 for the case of an alloy with about 40 percent 
nickel. 

The values of mass susceptibility as determined 
by the above methods are given in Table II. The 
measurements showing field dependence are 
indicated and for these the initial susceptibility 
is given. 

It is to be noted, as discussed in an earlier 
paper,’ that the relative precision of the values 
given in Table II for different specimens is not 
greater than about 1 percent because of the dif- 
ficulty in locating the specimens at exactly the 
same point in the magnet. On the other hand, the 
relative accuracy of the values obtained for any 
given specimen at different temperatures above 
room temperature is much greater than this, 
since the position of the specimen was not 
altered throughout the run. This position could 
not change appreciably as a result of thermal 
expansion since the heated zone in the furnace 
is only a few inches long. This conclusion is of 
importance in evaluating the small but con- 
sistent changes in susceptibility which are ob- 
served above room temperature. Fortunately, 
this conclusion is substantiated by the measure- 
ments themselves since it is found that both dia- 
and paramagnetic susceptibilities change in the 
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proper direction as the temperature is varied; 
that is, the paramagnetic values become greater 
when the diamagnetic values decrease. If the 
observed changes with temperature were due 
simply to a shift in position of the specimen, then 
the dia- and paramagnetic susceptibilities would 
both either increase or decrease depending on the 
direction of the shift. 


DISCUSSION OF RESULTS 


The susceptibility at room temperature is 
plotted as a function of nickel content in Fig. 2, 
and the results of previous workers are shown for 
comparison. Our results lie along a smooth curve 
which is in qualitative agreement with previous 
measurements* ® except for the work of Williams® 
who found a bulge at about 30 percent nickel. 
The quantitative agreement of the various 
workers is not too good. It seems likely that the 
reason for this is to be found in the purity and 
treatment of the alloys rather than in the 
methods of measurement. For example, at room 
temperature 0.1 percent iron will change the 
susceptibility of copper to the same extent as 
4 percent of nickel and hence the iron impurities 
of the alloys up to about 20 percent nickel may 
influence considerably the value obtained for the 
susceptibility. Above 30 percent nickel the sus- 
ceptibility increases very rapidly with concen- 
tration and, therefore, it becomes increasingly 
important to be sure that the alloys are well 
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Fic. 1. Magnetization curves [specific intensity of mag- 
netization (¢) ] at low temperatures. 


‘G. Gustafsson, Ann. d. Physik 28, 121 (1937). 
5M. Wheeler, Phys. Rev. 56, 1137 (1939). 
°E. H. Williams, Phys. Rev. 38, 828 (1931). 
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Fic. 2. Mass susceptibility as a function of nickel 
content at room temperature. 


homogenized before measurement. The fluctua- 
tion in nickel concentration throughout a speci- 
men as solidified from the liquid state can easily 
be great enough to give an erroneous value for 
the susceptibility corresponding to the average 
composition of the specimen. 

It is not beyond the realm of possibility that 
at least a partially ordered structure might occur 
in the composition range of Cu;Ni and that such 
an occurrence would affect the susceptibility. 
This might account for the irregularity in the 
susceptibility versus concentration curve as 
reported by Williams.* In order to check this, the 
authors slowly cooled the 30 percent nickel alloy 
from 820°C to 220°C over a period of ten days, 
since this should allow an ordered structure to 
form. The susceptibility of this specimen agreed 
with the value obtained previously to within the 
experimental error. As a further test the sus- 
ceptibility of the slowly cooled specimen was 
measured at 50°C intervals up to 450°C, since 
the breaking up of an ordered structure should 
give a discontinuous change in the observed 
values. No such effect was observed and hence 
it was concluded that the specimen did not have 
an ordered structure, or at least that if there 
were an ordered structure it caused less than a 
1 percent change in the susceptibility. 

From Fig. 2 it is evident that the susceptibility 
does not vary linearly with the concentration of 
nickel and that the deviation from linearity 
becomes more and more marked as the ferro- 
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Fic. 3. Non-linear variation of susceptibility with 
concentration in dilute solutions. 


magnetic compositions are approached. This 
behavior is true at all temperatures down to 
14°K and it is true even for small concentrations 
of nickel as shown in Fig. 3. 

The observed low temperature susceptibility 
values were corrected for the susceptibility of 
copper and the dissolved iron impurities, and the 
reciprocals of these quantities were plotted 
against the absolute temperature as shown in 
Fig. 4. It is apparent that a Curie-Weiss law does 
not describe the results for the low nickel alloys 
and that it tends to become more satisfactory 
as the nickel content approaches the ferromag- 
netic compositions. In view of this behavior no 
attempt was made to calculate Bohr magneton 
values for the nickel atoms by means of the 
Curie-Weiss equation although this was done by 
a different method as will be explained shortly. 

Because of the failure of the Curie-Weiss equa- 
tion to describe the results below room tem- 
perature, the measurements were extended to 
1300°K. The susceptibility values over the entire 
range of temperature are plotted in Fig. 5. It is 
immediately apparent that the susceptibility of 
all the alloys increases at high temperature, 
although this effect is concealed by the tendency 
toward a Curie-Weiss behavior in the higher 
nickel alloys. This high temperature behavior is 
undoubtedly the same phenomenon as_ has 


already been reported for pure copper’ although 


7 Bitter, Kaufmann, Starr, and Pan, Phys. Rev. 60, 134 
(1941). 


the effect seems to become much greater as 
nickel is added. 

For purposes of comparison the measurements 
of Gustafsson‘ are also shown in Fig. 5. It is 
apparent that the temperature range covered by 
his work did not allow him to observe the phe- 
nomenon now reported and that his work is in 
agreement with ours over the corresponding 
ranges of temperature. 

It is obvious from Fig. 5 that the Curie-Weiss 
equation cannot represent the results at high 
temperatures and hence an expression of the 
form 


x=aT+b+(c/T) (1) 


was fitted to the curves by substituting experi- 
mental values at three temperatures and solving 
for the constants. This procedure is quite arbi- 
trary, of course, and can only be justified if the 
curve calculated from the constants represents 
the experimental curve accurately. In the present 
case there was satisfactory agreement between 
the computed and the experimental curves for 
all the alloys from 10 to 35 percent nickel. That 
is, the difference between the two curves was not 
greater than a few percent except at 14°K. At 
this temperature and for the other alloys at other 
temperatures as well, the difference was as large 
at 10 percent. In spite of this it seems worth 
while to present briefly the results obtained and 
accordingly the constants are tabulated in 
Table IIT. 

The values in Table III were found when the 
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Fic. 4. Reciprocal of mass susceptibility corrected for 
copper as a function of temperature. 
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susceptibility data obtained at 20.4°K, 290°K, 
and 1270°K are used for the simultaneous solu- 
tions. A choice of different temperatures would 
lead to somewhat different values for the con- 
stants. Also, if one of the terms in Eq. (1) is 
small compared with the other terms, then that 
constant is relatively more subject to error 
because of the limited experimental accuracy. 
This reason presumably accounts for the erratic 
variation of the constants a and b when c becomes 
large. 

The first term in Eq. (1), which indicates an 
increase in susceptibility with increasing tem- 
perature, is not accounted for by any of the 
present theories of magnetism. Such behavior 
has been reported, however, for the alkali 
metals* and has been attributed by Stoner to a 
decrease in the number of electrons per unit 
volume caused by the thermal expansion of the 
crystal. This explanation is fairly consistent with 
the present observations on pure copper since 
it is possible to calculate a temperature coef- 
ficient for the mass susceptibility which is only 
about six times smaller than the value for a 
which is given in Table III. This calculation was 
made with the equation for the susceptibility of 
free electrons and on the assumption that the 
change in the number of electrons per unit volume 
is proportional to the temperature coefficient of 
thermal expansion for copper. Such a calculation 
is not nearly as satisfactory for the alloys since 
the coefficient of thermal expansion changes only 
slightly with addition of nickel to copper, but the 
constant a changes by at least a factor of 10 
when only 20 percent of nickel is added. How- 


TABLE III. Constants in Eq. (1) and Bohr magneton 
numbers (pz). 











Wt. 
on 

Alloy nickel a X10” bx<106 cX106 PB 
J27A 0.00 0.237* —0.097* 

JIA 0.94 0.217 —0.071 0.787 0.199 
J8A 4.99 0.67 0.017 1.756 0.129 
J26A 10.99 1.38 0.160 6.34 0.165 
J25A 20.3 2.22 0.444 13.29 0.176 
J24A 30.0 0.91 1.073 51.95 0.286 
J23A 35.6 0.31 1.302 197.7 0.512 
J22A 38.0 5.42 0.510 632.0 0.886 








* Values at 290°K, 770°K, and 1270°K used in simultaneous solution. 


SE. C. Stoner, Magnetism and Matter (Methuen & 
Company, London, 1934), p. 509. 
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Fic. 5. Mass susceptibility as a function of temperature 
for different nickel contents. 


ever, the other factors which may change with 
alloying such as the number of free electrons per 
atom and the shape of the energy band could 
easily mask any effect due to differences in 
expansion coefficient. 

The constant b of Eq. (1) may be considered 
as the sum of the diamagnetism and temperature 
independent electron paramagnetism. It is inter- 
esting to note, first of all, that this treatment of 
the experimental data leads to a value of —0.097 
x<10-* for the temperature independent mass 
susceptibility of copper instead of the ordinary 
room temperature value of —0.085X10-* and 
that this new value is in somewhat better agree- 
ment with the calculated value for copper.*® 

The values of } increase continuously with 
nickel content until c becomes so large that } 
can no longer be satisfactorily determined. At 


*E. C. Stoner, Magnetism and Matter (Methuen & 
Company, London, 1934), p. 511. 
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30 percent nickel 6 has become about 15 times 
as large as the theoretical free electron suscep- 
tibility of copper. In terms of the band theory of 
susceptibility this would mean that the number 
of states at the top of the energy band has in- 
creased by this factor since the electronic sus- 
ceptibility is proportional to the number of states. 
However, an increase of this magnitude hardly 
seems compatible with the published energy 
band shape for copper'® in spite of the fact that 
the alloys are approaching an unfilled d band. 
If the nickel and copper atoms of the alloys form 
a common set of bands, then the d band should 
not begin to be unfilled until 40 percent nickel is 
reached. The present determination of ) means 
either that this perfect sharing of electrons does 
not occur in the alloy or else one must conclude 
that 5 does not have the significance attributed 
to it here. 

The last term in Eq. (1) offers further evidence 
that a common 3d band for the copper and 
nickel atoms does not describe the results ac- 
curately. The reason is simply that this term 
indicates a normal paramagnetic susceptibility 
whereas the band picture would alloy only a tem- 
perature-independent electronic susceptibility up 
to 40 percent nickel. There can be no doubt that 
the low nickel alloys exhibit a large temperature 
dependence, but whether or not this is correctly 
represented by a term in (1/7) is open to dispute. 
At any rate the constant c as here determined 
shows a continuous variation with nickel content 
becoming rapidly larger as a concentration of 40 
percent is approached. If one regards c as the 
ordinary Curie constant, the Bohr magneton 
values per nickel atom may be obtained from the 
formula pg=2.84(cM/K)', where M is the 
molecular weight of nickel and K is the fraction 
of nickel by weight in the alloy. From Table III 
it can be seen that the Bohr magneton values 
remain fairly constant up to about 20 percent 


J. C. Slater, Phys. Rev. 49, 537 (1936). 


nickel and then increase rapidly. The significance 
of this is not clear since the values of pg are 
about ten times smaller than that found for pure 
nickel above its Curie temperature. 

Finally, we must discuss the field dependence 
of the susceptibility which was observed at low 
temperatures in the alloys with 30 percent and 
more of nickel. The careful measurements of 
Weiss and Forrer™ have shown that the mag- 
netization curve of pure nickel shows no discon- 
tinuity at the Curie temperature and that the 
curve becomes linear within the experimental 
error only after 40 degrees above the Curie tem- 
perature has been reached. It seems reasonable 
to suppose that this same phenomenon exists 
also in the ferromagnetic copper-nickel alloys. 
From this it follows that in varying the com- 
position from ferrromagnetic to paramagnetic 
alloys at any given temperature the paramag- 
netic compositions close to the boundary will 
show field dependence. Our observations must be 
part of this phenomenon in spite of the fact that 
the Curie temperature has been supposed to 
reach absolute zero at 40 percent nickel. 

The ordinary Weiss theory will lead to a large 
field dependence of the susceptibility above the 
Curie temperature as may readily be verified by 
numerical substitution in the equation. For 
nickel this calculated field dependence becomes 
smaller than 1 percent at 30°C above the Curie 
temperature and this is in agreement with the 
measurements. (The calculated magnetization, 
however, is only about one-half as large as the 
measured values in fields up to 20,000 oersted.) 
This result substantiates the idea that the ob- 
servations on low nickel alloys are to be corre- 
lated with the tendency toward ferromagnetism. 
A more detailed discussion of this problem will 
be undertaken at a later time after more obser- 
vations of the same sort have been made on 
similar alloy systems. 


'P. Weiss and R. Forrer, Ann. de physique 5, 153 
(1926). 





a= _~_ ao -_ tr a ac 


a alUhhFrlC<iC rr |)FelCceklC helCUCllClCO 


=a, —. ae -_ _~ — 








we \¥ 


PHYSICAL REVIEW VOLUME 63, NUMBERS tit AND 12 JUNE 1 AND 15, 1943 


Electrical Conduction and Recrystallization in Thin Pb Films 
Deposited at Low Temperatures* 


Epcar L. Armi** 
Cryogenic Laboratory, Norman Bridge Laboratory of Physics, California Institute of Technology, Pasadena, California 


(Received March 26, 1943) 


By evaporation in high vacuum, thin Pb films were deposited on glass plates cooled by liquid 
hydrogen. Beginning of electrical conduction was observed in films of a nominal thickness of 7A. 
Aging of these films at various temperatures revealed the existence of two regions of different 
behavior; decrease of resistance below and increase of resistance above an equilibrium tem- 





perature. Possible explanations are offered. 





S numerous workers have shown, thin metal 
films condensed out of their vapor are 
unstable and undergo structural changes while 
aging. Various factors have an influence on these 
changes by affecting the mobility of the metal 
molecules involved, in particular the nature of 
the backing, thickness of the film, temperature, 
presence of gases, etc. Thus either an increase or 
decrease of the electrical resistance of the film 
may be observed. 

Starting with investigations by Anderson! and 
Foster® on the electrical conductivity of thin Ag 
and Pb films, the range of the present research 
was extended to temperatures of liquid and solid 
hydrogen. It was the purpose of this investigation 
to study the behavior of thin films of heavy 
metal atoms deposited at temperatures suffi- 
ciently low to eliminate disturbing thermal 
effects. These films were found to possess a 
“‘critical’’ temperature (in the neighborhood of 
100°K or higher, depending on the individual 
data of the film investigated), below which the 
electrical resistance decreases with time until it 
reaches a stable value, and above which it 
increases at a rate rapidly growing with higher 
temperatures. 


APPARATUS AND METHODS 


Figure 1 shows a lengthwise cut through the 
evaporation chamber. Constructed of Pyrex 
glass, it has wide bore ground joints at both ends. 


*E. L. Armi, Thesis, June 1941, California Institute of 
Technology. 

** Now at the Department of Physics, Alabama Poly- 
technic Institute, Auburn, Alabama. 

1A, B. C. Anderson, Thesis C.I.T. (1933); A. B. C. 
Anderson and A. Goetz, Phys. Rev. 45, 293 (1934). 

2M. G. Foster, Phys. Rev. 57, 42 (1940) and Thesis 
C.1.T. (1939). 


The Pb to be evaporated is held in an iron 
crucible and placed in an electrical furnace. Stray 
evaporation is cut only by a cylindrical shield 
around the axis of the beam. Two shutters are 
provided and operated by an electromagnet from 
the outside, the first closing the opening of the 
cylindrical shield, the second immediately in 
front of the deposit surface. 
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Fic. 1. Evaporation chamber. Lengthwise cut and detail 
of shutter and electrodes. 
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Fic. 2. Continuous deposits. Curve I room temperature 
deposit. Curve II low temperature deposits. 


Deposits were made on 0.2 mm microscope 
cover glasses provided with platinized contact 
strips, the film area being 1.3 cm by 1.4 cm. The 
backs of the cover glasses were platinized and 
copper plated before being inserted; in order to 
assure good heat transfer with the copper cooling 
reservoir, they were soldered to the latter, in 
vacuum, by induction. 

Film thicknesses were calculated from the ex- 
perimental evaporation data, with the density of 
Pb taken as that of the metal in bulk (11.3 
g/cm*). Evaporation rates were established by 
running the furnace at fixed temperatures for 
extended periods of time and weighing on a 
microbalance the deposits obtained; the values 
thus found (between 1 X 10-8 and 1X 10-* g cm~ 
min.—' for temperatures in the neighborhood of 
500°C) were in good agreement with those cal- 
culated from kinetic theory.* 

Vacua were measured either with an ionization 

’ The formula for effusion from a small opening was used. 
p= (dN/dt)(2MRT)*(760) /dSdw cos 0(1.0132 X 10°). In this 
p is taken from the vapor pressure formula logis p= 
—52.23A/T+B. (A =188.5, B=7.82) (p vapor pressure in 
mm Hg, M molecular weight, R gas constant per mole, 
T absolute temperature, dN /dt number of moles per second 


effusing from an opening of area dS cm? into a solid angle 
dw, at an angle @ with the normal to dS.) 


gauge connected directly to the evaporation 
chamber or with a McLeod manometer. A high 
frequency induction furnace was used to outgas 
metal parts before evaporation was started. All 
glass parts were carefully cleaned and a fresh 
cover glass inserted for each run. Vacua obtained 
were of the order of 10-* mm Hg. 

The electrical circuit, described previously by 
Anderson! and Foster,’? permits resistance deter- 
minations up to 10" ohms by use of a leakage 
compensation method. Five condensers of dif- 
ferent sizes are provided to cover a wide range of 
resistance values. A string electrometer serves as 
null instrument. 

The experiments fall into two groups, depend- 
ing on the temperature of the film. In the first 
group, deposits were made on a backing at room 
temperature. In the second group, the tank 
behind the backing was filled with a cooling 
agent. Liquid air was siphoned directly into the 
reservoir and, by reducing the pressure over the 
liquid air, it was solidified; similarly, after pre- 
cooling with liquid air, liquid hydrogen was 
siphoned into the revervoir and solid hydrogen 
was obtained by reducing the vapor pressure. 
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Fic. 3, Interrupted deposit (room temperature). Solid 
lines show depositing periods, broken lines interruption 
periods. At A, a small quantity of air was admitted. 


























ELECTRICAL CONDUCTION 


RESULTS AND DISCUSSION 
(a) Room Temperature Deposits 


Curve I on Fig. 2 may be regarded as typical 
for room temperature deposits. At the beginning 
of the depositing period, the resistance measured 
is that of the glass plate at room temperature; 
after heat influx from the furnace has acted for 
several minutes, a small drop in the resistance of 
the glass plate is observed. No effect of the 
deposited Pb is seen until the film has reached 
a considerable thickness; in this case only after 
80 minutes depositing time (corresponding to a 
thickness of the film of 160A) a point is reached 
where the resistance decreases rapidly with each 
new layer of metal atoms added. 

Figure 3 shows how the resistance of a room 
temperature film increases when the deposit is 
interrupted (dotted lines); the increase is the 
more rapid the higher the resistance at the time 
of interrupting. At A, a small amount of air is 
admitted during the depositing process (to a 
pressure of 10-* mm Hg), with the immediate 
effect of stopping the decrease of resistance; only 
after a considerable length of time does the 
resistance resume its downward move. 

On Fig. 4, the behavior of films under the 
influence of air is studied in more detail. While 
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Fic. 4. Air admission to room-temperature films. Film I, 
initial resistance 210° ohms; film II, initial resistance 
6X10° ohms. At A, B, and C, air was admitted to the 
pressures indicated. 


in film I (initial resistance 2X10’ ohms) the 
resistance of the background is reached in a few 
minutes, a considerable conduction is maintained 
in film II (initial resistance 6105 ohms) after 
exposure to air even at full atmospheric pressure, 
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It can be assumed that, apart from the phe- 
nomena discussed further below, gas adsorption 
and oxidation play an important role. Thus, 
while very thin films are easily attacked by gases 
even under reduced pressure,‘ heavier films may 
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Fic. 5. Aging of several low-temperature films at low 
temperatures. 


conserve their inner layers intact for electrical 
conduction. Also, while a film is being deposited, 
bound gas atoms will impede the formation of 
continuous metallic layers. The maintenance of a 
high vacuum during the depositing process is 
therefore imperative, as the thickness at which 
conduction sets in will depend largely upon the 
vapor pressure in the evaporation chamber at 
the time of deposit. 


(b) Low Temperature Deposits 


Curve II on Fig. 2 represents a group of low 
temperature deposits, with the backing cooled 
by liquid hydrogen. Here conduction sets in at 
thicknesses as small as 7A. Warming up the films 
to room temperature destroys all traces of con- 
duction; it is thus possible to deposit several 
films on top of one another, having essentially 
the same curve of resistance against time, with 
only a lengthening of the exposure time required 
to obtain conduction in successive deposits. 

Figure 5 displays the behavior of some of these 
films, after the molecular beam is shut off. In all 
cases a decrease of resistance is found, in contrast 


4M. C. Johnson and T. V. Starkey, Proc. Roy. Soc. 40, 
126 (1933). 
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Fic. 6. Aging of several low temperature films at 
raised temperatures. 


to the observations on room-temperature de- 
posits. After about one hour, stable values of 
resistance are finally reached, at a level which 
depends upon the initial resistance and the tem- 
perature of the film. Raising the temperature 
will fasten the ‘‘reaction’’ and reduce the value 
of the final resistance reached. 

Plotting the observed decrease of resistance 
for any given temperature, it is found that a law 
of the form 


log (R—Rr)=log (Ro— Rr) — Krt* 


is obeyed, in which Rr is the final resistance and 
Kr the “reaction constant” for a given tempera- 
ture, and the exponent s has a value of approxi- 
mately 3. 

Turning to Fig. 6, it is seen that by raising the 
temperature to higher values than in Fig. 5, an 
increase instead of a decrease of the electrical 
resistance is observed, the rate of this increase 
going up rapidly with the temperature. A “‘cri- 
tical’ temperature is found, separating the two 
regions of opposite behavior; this “‘critical’’ 
temperature depends upon the initial resistance 
of the film (and hence upon its thickness and 
structure). 

To interpret this behavior, it may be assumed 
that the reduced mobility of Pb molecules at 
low temperatures permits only small-scale re- 
grouping from the random positions into which 
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these molecules had been deposited; thus, by a 
slow formation of crystallites and bridges, a 
gradual increase of the electrical conduction is 
achieved. With raised temperature and increased 
mobility, that process is speeded up, resulting in 
a maximum conductivity at the ‘“‘critical’’ tem- 
perature. Above this temperature, the agglomera- 
tions, growing larger at the expense of their 
immediate neighborhood, become surrounded by 
zones of relative depletion, which in turn diminish 
the possibilities for electrical conduction. Finally, 
at high enough temperatures, no current will 
pass, as the film then is broken up entirely. 

Accordingly, in films deposited at higher tem- 
peratures, much heavier deposits will be needed 
before conduction can be established, since the 
agglomeration and recrystallizing process oc- 
curring will inhibit any film formation until 
enough material is deposited. The thickness at 
which conduction sets in, in thin films, may thus 
be linked to the mobility of the constituent 
molecules. Complicating the phenomenon is the 
presence of bound gas atoms which act as a 
hindrance to the formation of a continuous metal 
film. The results obtained here seem to confirm 
very satisfactorily the views expressed in recent 
papers by Suhrmann and Berndt® and Weber and 
O’Brien. ® 

Mention may finally be made of the influence 
of electric currents upon the conductivity of the 
films. Curve III in Fig. 6 shows how, by passing 
a current of 5 mA, the resistance is made to drop 
by several powers of 10 and the ‘“‘critical’’ tem- 
perature increased by about 40°, a phenomenon 
which may be explained by the high local current 
densities occurring. 

The author wishes to express his appreciation 
to Professor A. Goetz for his stimulating direc- 
tion of the work, to Dr. A. B. Dember for un- 
tiring help and advice and to Drs. A. B. C. 
Anderson and M. G. Foster, whose apparatus 
and methods were used as a starting point in this 
investigation. 


5 R. Suhrmann and W. Berndt, Zeits. f. Physik 115, 


17 (1940). 
6A. H. Weber and D. F. O’Brien, Phys. Rev. 60, 574 


(1941). 
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MINUTES OF THE MEETING OF THE NEW YORK STATE SECTION HELD 
AT ROCHESTER, NEW YorRK, APRIL 2-3, 1943 


HE Spring Meeting of the New York State Section of the American 
Physical Society was held at the Strong Memorial Hospital and at the 
Memorial Art Gallery, Women's Campus, University of Rochester, on Friday 
evening, April 2, 1943 and Saturday, April 3, 1943. An attendance of over 100 
members and visitors was recorded at each of the sessions. 

On Friday evening the Section was guest of the University of Rochester 
School of Medicine and Dentistry. Dr. Stafford L. Warren discussed the work 
of the Biophysical Laboratories of the School of Medicine, and after the lecture 
trips were conducted through the various laboratories. 

At the Saturday meeting the following program of invited papers was 
presented : 


Address of Welcome. PRESIDENT ALAN VALENTINE, University of Rochester. 
Microwaves. GEORGE K. Town, Stromberg- Carlson Manufacturing Company. 
Vacuum Pumps in the Development of Industry. K. Hickman, Distillation Products. 
Resistance of Glass to Thermal and Physical Shock. C. D. OuGHToN, Bausch and Lomb Optical 
Company. 
Animal Electricity. W. O. FENN, School of Medicine and Dentistry, University of Rochester. 
Discussion: Education of Physicists for War and Post-War Conditions. 

Colleges and Universities. R. C. Gisss, Cornell University. 

Secondary Schools. E. L. MANNING, State Education Department. 

Industry. J. T. LirrLeton, Corning Glass Works. 

W. R. FREDRICKSON 


Secretary 
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MINUTES OF THE MEETING OF THE METROPOLITAN SECTION HELD AT 
COLUMBIA UNIVERSITY, NEW YorRK, APRIL 2, 1943 


MEETING of the Metropolitan Section of the American Physical Society, 

being the only one for the season 1942-1943, was held on Friday, April 2, 

1943 at Columbia University. The meeting began with dinner in the Men’s 

Faculty Club. The session for the presentation of papers was held in the Pupin 

Physics Laboratories beginning at eight o’clock. The following papers were 
presented : 


Contributed Paper: Formation of Centers of Condensation in Supercooled Phases. PeTer G. 
BERGMANN (Abstract below). 

Invited Paper: The Elasticity of High Polymers as Bearing on the Synthetic Rubber Problem. 
HERMANN Mark (Introduced by Karl K. Darrow). 


The following officers were elected for the season 1943-1944: Chairman: 
HELEN A. MESSINGER; Vice Chairman: L. H. GERMER; Secretary-Treasurer: 
W. S. Gorton; Elected Members of the Executive Committee: W. E. LAMB, JR., 
C. M. SLACK. 

W. S. GorTON 
Secretary-Treasurer 
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a “critical’’ size below which the droplet is more likely to 


Heaviside calculus. 


The Formation of Centers of Condensation in Super- dissolve than to grow. To overcome this difficulty, the . 
cooled Phases. PETER G. BERGMANN, Lehigh University— author has considered the change to which an arbitrary size ' 
In the investigation of the spontaneous formation of amore _ distribution of droplets or crystals is subject in the course : 
stable phase in the interior of a supersaturated or super- of time. If the number of droplets containing between yu ' 
cooled phase, it is usually assumed that a count of the and u+dy moles at the time ¢ is Z(u, t)du, then Z satisfies 
centers of condensation which become visible per unit time the equation 
furnishes us directly with the rate at which new centers are aZ 1 2 OZ , Z df(u) ) 
being formed. This assumption is objectionable, because = MSP) (= a o)|. ; 
new centers will not necessarily continue to grow, but may _ where o(u) is the rate at which droplets of size » annex or 
redissolve before they become visible. In fact, it is well drop one molecule, and f(u) is the free energy per droplet. ‘ 
known that for every degree of supercooling, there exists It is proposed to discuss this equation with the help of t 

: 
i 
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MINUTES OF THE MEETING OF THE SOUTHEASTERN SECTION HELD AT 
AUBURN, ALABAMA, APRIL 2-3, 1943 


HE ninth annual meeting of the Southeastern Section of the American 
Physical Society was held at the Alabama Polytechnic Institute, Auburn, 
Alabama, on Friday and Saturday, April 2-3, 1943. Approximately sixty-five 
members and guests attended the meeting. Local arrangements were made by 


a committee headed by Fred Allison. 


The regular program consisted of eight papers, abstracts for five of which 


~ As = oe es UrrlUlUC ClO 


-_ =~ w> — AS 


are appended hereto. Abstracts of other papers will be found in the June issue 
of the American Journal of Physics. Additional features of the program were 
a talk by M. H. Trytten on ‘“‘The Manpower Situation in the Profession of 
Physics’”’ and a panel discussion on ‘‘Physics and the War Emergency,” led 
by A. E. Ruark. 

At the business meeting the election of the following officers for the year 
1943-44 was announced: Chairman, C. W. Edwards; Vice Chairman, R. C. 
L. Mooney; Secretary, E. Scott Barr; Treasurer, C. B. Crawley; and Member 





of the Executive Committee, K. L. Hertel. 
The invitation extended by the University of Tennessee to hold the 1944 
meeting of the Section in Knoxville, Tennessee, was accepted. 


E. Scott BARR, Secretary 


ABSTRACTS OF CONTRIBUTED PAPERS 


1. The Single Scattering of Fast Electrons in Light 
Gases. T. N. GAUTIER,* University of North Carolina.— 
The single scattering of RaE electrons has been studied in 
air, hydrogen, and argon. The scattering angles were 
mainly between 10° and 30° and the primary energies were 
between 200 kev and the upper limit of the RaE spectrum. 
Compared to the scattering in argon, the scattering in 
hydrogen and air was found to agree with the predictions 
of the Mott formula for nuclear scattering and the Mdller 
formula for electronic scattering, within the limits of 
instrumental and statistical error. The probable error was 
about 5 percent for the hydrogen-argon comparison, and 
about 2 percent for air-argon comparison. 


* Now at the National Bureau of Standards. 


2. Measurement and Analysis of Noises by an Improved 
Method. Danret S. ELLiott anp JoHn L. McLucas, 
Tulane University—In a previous paper an improved 
method for studying noises was outlined. This method 
consists primarily of two parts: (1) a device for tabulating 
rapidly varying intensities so as to determine the manner 
in which the intensities are dispersed; and (2) facilities 
for recording noise samples on a phonograph disk “in the 
field” and subsequently reproducing them in the laboratory 
to be analyzed by means of a commercial sound analyzer. 
In the present study a number of noises of various kinds 
have been analyzed by this method. Furthermore a 
modification has been made in the apparatus which now 
makes possible studies where the ordinary commercial 
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electric lines are unavailable. Comparisons have been 
made of the same situations using the energy from the 
ordinary commercial electric lines and from the new 
modification which uses a portable energy source. 


3. Generalization of Bateman’s Differential Equations. 
ARTHUR E. Ruark, University of North Carolina.—Bate- 
man worked out the system of differential equations 
governing the functions W,(¢t), which represent the chance 
that events will occur in time ¢; the chance that an event 
will occur in time dt is the only datum required to make 
the solution definite. The generalization of Bateman’s 
ideas to the case where W, depends on several variables 
is briefly presented. 


4. X-Ray Diffraction Measurements with Geiger- 
Mueller Tube. R. Pepinsky, Alabama Polytechnic Insti- 
tute, AND P. Wetsz, Bartol Research Foundation.—The 
time-rate of development of crystalline conditions due to. 
physical or chemical action in materials is often of interest. 
Examples are the dependence of hysteresis or drift in 
rubber or other polymers upon the rate of molecular 
alignment due to stretching, the chronology of polymorphic 
changes in crystalline forms, or structural changes in 
alloys with various treatments. Because of exposure times 
required for photographic recording of x-ray diffraction 
patterns, reactions occurring in minutes or seconds have 
in general been unobservable. The possibility of use of 
Geiger-Mueller tubes as rapid indicators of diffraction 
maxima and their history during reactions was explored. 
A preliminary survey was made with a compact GM-tube 
circuit for audible counting of x-ray quanta, the apparatus 
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being of general interest because of its simplicity and wide 
applicability. A GM tube with argon-ether filling and 
0.0005’ Pyrex window was used, the x-ray beam pene- 
trating the length of the cylinder. The counter supplied 
positive pulses to the grid of a tube normally biased to 
cut-off; with practically no current demanded by the 
circuit, its B voltage could be obtained from the GM-tube 
high voltage supply. Observations with this apparatus 
indicated that the GM tube can be successfully utilized 
for the diffraction studies. 


5. The Near Infra-Red Absorption Spectra of Some 
Vegetable Oils. E. S. BARR AnD W. R. Harp, Tulane 
University—The near infra-red absorption spectra of 
castor, cottonseed, olive, palm, and peanut oils have been 
measured in the region between ly and 5yz. A pronounced 
absorption at 3.434 was observed in each case and at- 
tributed to the C-H valency vibration. The other observed 
absorption positions did not in general appear at the 
same positions, which permits differentiation among these 
materials. Since cell thicknesses of 0.06 mm proved best, 
only small amounts of the oils were needed. The changing 
composition of such oils with aging seemed likely to 
produce corresponding changes in spectra. This was 
investigated, particularly for peanut oil, and marked 
changes were found in the samples subjected to various 
amounts of aging. Differences resulting from variations in 
the preparation of the oils were considered. For example, 
natural and “‘tasteless’’ castor oil spectra differed in the 
region of 2.54. Work on these materials is being continued, 
and it is planned to include other oils of vegetable origin 
in the general project. 
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MINUTES OF THE MEETING OF THE NEW ENGLAND SECTION HELD AT 
WELLESLEY, MASSACHUSETTS, APRIL 3, 1943 


HE twenty-first regular meeting of the New 
England Section of the American Physical 
Society was held at Wellesley College, Wellesley, 
Massachusetts, on Saturday, April 3, 1943. Forty- 
six members were in attendance. The invited 
papers were as follows: 


The M.K.S. System of Units. WeNnpeLL H. Furry, 
Harvard University. 

Physics at Wellesley College. Louise S. McDoweE Lt, 
Wellesley College. 

“Automatic Volume Control’ of the Ear. Haro.p P. 
Knauss, Ohio State University. 

Entropy and Probability. K. K. Darrow, Bell Telephone 
Laboratories. 


Four short papers were presented, one of these 
being by title. Those pertaining to research are 
as follows: 


Radiations from the 100-Day Arsenic Activity. L. G. 
E.utiott AND M. Deutscu, Massachusetts Institute of 
Technology, Cambridge.—Sagane et al.' have isolated an 
arsenic activity of about 100-day half-life from germanium 
bombarded with deuterons which they assigned to As” 
because it was found to emit negative electrons. By study- 
ing the radiations in the magnetic lens beta-ray spectrom- 
eter we found these electrons to be due to the internal 
conversion of a 0.052+0.003-Mev gamma-ray. The con- 
version in the K shell was found to be about 5 times as 
strong as that in the L shell. K x-rays were also observed. 











In order to determine whether the gamma-rays follow 
orbital electron capture or are due to an isomeric transition, 
coincidences between conversion electrons and K x-rays 
were observed. The electrons due to conversion in the L 
shell were found to be accompanied, on the average by half 
as many x-rays as those from the K shell. We conclude 
from this that the 100-day activity decays by orbital elec- 
tron capture followed, within less than 10~* second, by a 
highly converted 0.052-Mev gamma-ray. The previous as- 
signment must therefore be in error. Positrons are emitted, 
if at all, in less than 2 percent of all the disintegrations. 


1 Sagane ef al. Proc. Phys.-Math. Soc. Japan, 21, 660 (1939). 


Absorption Bands in 240-185 my Region in Spectra of 
Aromatic Amino Acids. GLApys A. ANSLOW AND Sawa C. 
NassAR, Smith College—An analysis of the end absorption 
in the spectra of the aromatic amino acids reveals three 
independent bands, which in phenylalanine set in near 225, 
205, and 198 mu. In tyrosine, tryptophane, and also in 
histodine these are shifted to longer wave-lengths, more in 
tyrosine than in the other two. The latter two bands have 
intensities and locations similar to the 205 and 185 mu 
bands typical of the cyclic hydrocarbons, and may originate 
from energy absorption in the respective rings. These bands 
shift more than the weaker 225 my band, which tends to 
merge with the 205 band. A detailed investigation of the 
longer wave-length band, which occurs in all fatty acids 
between 225 and 240 my, reveals that the tyrosine band is 
unique in having a long wave-length component. This must 
result from energy absorption in the ring near the OH 
substitution, as is also indicated by its existence in the 
spectrum of hydroquinone. Thus, the seat of the absorption 
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is in the C—C—OH configuration, and is independent of 


xX 
the nature of X. 


Definition and Measurement of Absolute Electrode Po- 
tentials. Davip C. GRAHAME, Amherst College (Read by 
title) —The absolute electrostatic potential of a point in 
space is defined unambiguously in terms of the properties 
of a unit plus charge. By a slight extension of this definition, 
the electrostatic potential of a conducting surface can also 
be defined unambiguously. But the potential of a point 
within a conductor varies enormously from point to point 
if one takes into account the discrete structure of matter, 
as one must in many problems. To avoid this difficulty 
one may define the potential of an isotropic conductor as 
the potential of a point in a cavity within the conductor. 
This definition is in accord with most usage and does not 

-conflict with the basic definition of potential mentioned 
above. Its particular advantage is that it makes possible an 
experimental determination of the electrostatic potential differ- 
ence between two phases. By a very slight modification of the 
original Kelvin method it becomes possible to carry out 
measurements of liquid junction potentials and of absolute 
electrode potentials without theoretical difficulties. From 
such measurements one can calculate individual ionic ac- 
tivities, energies of hydration of gaseous ions, and related 
quantities. 


Another paper, concerned with instruction, has 
been sent to the American Journal of Physics. 
MILDRED ALLEN, Secretary-Treasurer. 
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And of its Ohio Section 


CoLuMBUS MEETING, APRIL 30 AND May 1, 1943 


HE spring meeting of 1943, being the 254th K. Lark-Horovitz, Purdue University: Biological Ap- 
meeting of the American Physical Society, Plications of Radioactive Tracers. 


ace : Henry Crew, Northwestern University: Galileo, the 
was held in the buildings of the Ohio State Pioneer Physicist. 


University in Columbus on Friday, April thir- Louis T. More, University of Cincinnati: Newton’s 
tieth, and Saturday, May first. As all but our Philosophy of Nature. 
very newest members will recall, the spring R. BowLInG BARNES, American Cyanamid Company: 


Applied Infra-Red Spectroscopy. 

: . . . J. R. DowninG, DuPont Experimental Station: Appli- 
held in Washington, and only in 1942 was it cation of Eafra-Red 8 ay to Chantest Resussdh. 
transposed to Baltimore, the latter city becoming H. H. Niesen anv Ety E. BELL, Ohio State University: 


itself overcrowded before the advent of 1943. Automatic Recording Vacuum Infra-Red Grating Spec- 
Though not distrustful of the hospitality of the trometer. 

Ohio State University, the officers of the Society 
were hardly overconfident of the prospects of a 


meeting was for many years without exception 


Sixteen ten-minute papers were contributed, 
three of these (2, 10, 11) being read by title: 


meeting held so far from its traditional place, Tho shutuuctn ene eemaniiel ot ten cael tae 
and their fears seemed at first to be justified by ‘ PPS 
proceedings. 


the were ot contributed —- The am The dinner of the Society was held on Friday 
dispelled their premonitions: The meeting was ea” ; ape 
- evening in the Deshler-Wallick Hotel, President 
large by any standards excepting only those of a , 
sien Winsiiines ti f th t with 166 Hull presiding. The after-dinner programme 
en Lee See women re consisted of talks by Charles F. Kettering of the 


registrants and two hundred in attendance at General Motors Research Corporation and by 
the dinner. The Local Committee, headed by R. L. Edwards of Miami University 


Dean Alpheus W. Smith, made excellent prepa- 
rations for the meetings. The Ohio Section of 
the Society, and Section F of the Ohio Academy 
of Sciences, merged their meetings with ours. 
Part of the success of this meeting was due to 
the invited papers, the givers of which are to be 
thanked for their speeches and for the journeying 
they did. Those of R. B. Barnes, J. R. Downing, 
H. H. Nielsen,-and E. E. Bell formed a sym- 
posium on applied infra-red spectroscopy, ad- 
mirably introduced by the first-names; Norman 


Wright took the chair. Two—those of Henry ~~ -T.. Society has lost through death one Fellow 
Crew and L. T. More—would have been given : : ms a 

‘ (A. F. Kimball, General Electric Company) and 
at the New York meeting of December, 1942 as . 

, : one Member (H. A. Pidgeon, Bell Telephone 
part of the Galileo-Newton symposium of the 
A.A.A.S., had that organization not cancelled 
its meeting. It is a satisfaction to have provided Elected to Fellowship: P. Pringsheim. 

a hearing for these papers; I am informed that Elected to Membership: Kenneth L. Andrews, Ernest Carl 


they are to : cienti NM _  Barkofsky, Gordon Edward Becker, Winston H. Bostick, 
y to appear in the Scientific Monthly Wayne M. Bowers, Joseph M. Brogan, W. Byron Brown, 


The titles of these five, and that of K. Lark- Royal C. Bryant, Richard C. Burns, Robert Van Duyne 
Horovitz which opened the invited programme, Campbell, H. W. Cobb, Frederick Coensgen, Philip 
are subjoined: Cooperman, Charles D. Coryell, Donald W. Crouch, Elsie 
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The Council met at intervals throughout 
Friday and Saturday. It elected one candidate 
to Fellowship and seventy-five to Membership. 
Their names follow. On receipt and after dis- 
cussion of a petition signed by forty Fellows 
and Members, the Council voted to establish a 
Division of Electron and Ion Optics in the 
American Physical Society under Article LX of 
the Constitution thereof. This Division is now 
in process of organization, and members wishing - 
to join it will shortly have an opportunity. 


Laboratories). 
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Margaret Dollmann, Richard Ehrlich, Milton Walter 
Franke, Urie Funaroff, Kenneth William Gardiner, 
Thomas N. Gautier, Jr., H. E. Greene, Leonard Gross, 
Dwight T. Hamilton, E. J. Hellund, J. Ross Heverly, 
David Lawrence Hill, Rosalie C. Hoyt, Albert A. Hudimac, 
Margaret Catherine Jenkel, W. E. Johnson, Albert E. 
Keheler, Jr., John M. Kopper, Ralph Werner Krone, 
Jerome Kurshan, Theodore James La Chapelle, Jr., Leo 
R. Landrey, Christian LaPointe, Howard Carlton Latham, 
Lester E. Lighton, Gustave A. Linenberger, Jr., Bernard 
Lippmann, Luther Leo Lowry, Arthur L. Lutz, Lloyd B. 
Macurdy, I. Paul Maizlish, N. Marcuvitz, Roland E. 


Meyerott, Charles Frank Miller, James F. Miller, Jr., 
Arnold R. Moore, Roger M. Morrow, Vernon W. Myers, 
I. Opatowski, Doris H. Paul, T. F. Protzman, Robert S. 
Rasmussen, Charles Ravitsky, W. Joseph Riegger, Emery 
H. Rogers, Sylvan D. Rolle, Ben Z. Rubin, James J. 
Shipman, William C. Skinner, O. J. Snow, H. Cecil Spicer, 
Edwin Garnet Stanford, Lorenzo Sturkey, George Clifford 
Sziklai, J. A. Van Horn, Frank Verbrugge, C. V. Weaver, 
Mitchel Weissbluth, E. M. Wise, Adrian Ziolkowski. 


KarRL K. Darrow, Secretary 


ABSTRACTS OF CONTRIBUTED PAPERS 


1. On the Linear Stress-Strain Curves of the Hexagonal 
Metals. J. S. KOEHLER, Carnegie Institute of Technology.— 
Taylor! supposed that in the work-hardened state a crystal 
contains many line dislocations, all of which are parallel to 
one another. He obtained the parabolic stress-strain curve 
of the face-centered cubic metals by assuming that the 
dislocations are present in a two-dimensional square lat- 
tice. We have been able to obtain linear stress-strain curves 
by using a rectangular dislocation lattice. The resulting 
formulas have been used to fit the measured stress-strain 
curves of single crystals of Mg, Zn, Cd.? The dislocation 
spacing in highly work-hardened hexagonal metals is about 
10-* cm in agreement with previous estimates. In the case 
of Cd at low temperatures the stress-strain curves show a 
sudden change of slope. Our expressions fit the lower 
straight-line section only. The change in slope may possibly 
be due toa change in the mechanism of plastic deformation. 
It may be that twinning plays a role in the deformation 
occurring at large stresses. 


1G. I. Taylor, Proc. Roy. Soc. 145, 362 (1934). 
2 E. Schmid and W. Boas, Kristallplastizitat (Springer, 1935), p. 158. 


2. Cathode Drop and Current Density in Welding Arcs. 
G. H. Fett, Department of Electrical Engineering, Univer- 
sity of Illinois——Experiments to determine the cathode 
drop of welding electrodes in air at atmospheric pressure 
and to compute the electron density in the arc were con- 
tinued using the method previously described.! The cathode 
drop for certain d.c. welding electrodes used for hard- 
surfacing was greatly influenced by the coating, the values 
changing from 9.4 to 14 volts after the coating was re- 
moved in the extreme case. Mild steel electrodes showed 
little or no change with coating removal. The electron 
density computed values was not proportionately influ- 
enced by the coating. Experiments involving currents less 
than one ampere indicate that with these electrodes the 
current density increases with a decrease in current. At 
values of current above 5 amperes the current density is 
constant for any given electrode, and is influenced by the 
type of electrode more than the coating. The energy lost 


by the arc due to convection varies directly with the arc 
perimeter. Hence the current density must increase as the 
current decreases in order that the arc temperature within 
the core be held constant. This proportional increase in loss 
with reduction in current may account for the fact that arc 
currents below a certain minimum, usually about 0.3 am- 
pere, cannot be maintained despite the use of fairly high 
potentials. 


1G. H. Fett, Welding J. Res. Supp. 21, 27 (1942). 


3. Magnetic Potential Energy. F. W. WarBurTON, Uni- 
versity of Kentucky.—Treatment of magnetic energy as ex- 
ternal (potential) energy in the sense of a modification of 
electrostatic energy due to electron motion, makes use of 
the classical Lagrangian function for the energy of the elec- 
tron in a magnetic field. However, a Hamiltonian function 
defined in terms of the mechanical momentum rather than 
the classical canonical momentum differs from the classical 
Hamiltonian by including explicitly the magnetic poten- 
tial energy, and equals the total energy. The change in 
magnetic potential energy upon adding a magnetic field is 
equal in a first approximation to the Larmor precession 
increase in rotational kinetic énergy. The two give a change 
in energy proportional to e/mc and provide a simple inter- 
pretation of the gyromagnetic experiments and the so- 
called electron spin. The omission of the canonical mo- 
mentum and the proof that the Hamiltonian remains 
constant in time, is not violation of conservation of energy, 
and it permits an electron-proton pair to release part of 
their energy in radiation in accord rather than at variance 
with the added requirements of quantum theory. 


4. A Physical Mechanism for the Nebular Red Shift. 
L. H. THomas, Ohio State University.—It seems possible to 
account for the nebular red shift by a second-order effect 
of ordinary coherent refraction. When a wave train, re- 
garded as plane and infinite in lateral extent, but finite in 
length, advances through a refracting medium, the parts 
of which are not elastically anchored in position, the me- 
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dium receives as the front of the wave train advances over 
it an average velocity in the direction of propagation of 
the wave train. This results from the second-order effect 
of the magnetic field of the wave train on the first-order 
transverse velocities, lagging in phase, of the material par- 
ticles. If there is no “‘frictional’’ resistance, this velocity is 
lost as the back of the wave train passes, the medium is left 
undisturbed, and the wave train is propagated undisturbed. 
If “frictional’’ resistance (which may be the statistical 
effect of waves passing in other directions) damps out the 
velocity in the direction of propagation, the wave train as 
a whole is gradually reddened. The whole process is co- 
herent and occurs even in the approximation of geometrical 
optics. It would not be expected that it should be possible 
to observe it under laboratory conditions, because a ve- 
locity of the medium in the direction of propagation of the 
wave train is then opposed by elastic restoring forces. 


5. A Recording Relay For An Infra-Red Spectrograph. 
Ety E. BELL, Ohio State University—A recording photo- 
electric bridge has been constructed for use with an infra- 
red, grating spectrograph which is able to reproduce the 
2.7-mu water vapor band with detail equivalent to that 
previously obtained with point-by-point measurements 
under similar conditions. A single base, mounting a light 
source, a Moll bifilar galvanometer, and a pair of photo- 
tubes is supported inside of a heavy metal tank containing 
the spectrograph. Here the relay is relatively free from 
temperature fluctuations and from vibrations. An unbal- 
ance of light on the photo-tubes is applied to commercial 
recording milliammeters through a bridge circuit similar 
to that of McAlister, Matheson, and Sweeney' except that 
an additional stage of amplification has been added. A 
provoked potential of 10-* volt in the thermocouple is 
recorded as readily distinguishable from the background. 


1E. D. McAlister, G. L. Matheson, and W. J. Sweeney, Rev. Sci. 
Inst. 12, 314 (1941). 


6. The Infra-Red Spectrum of Methyl Fluoride at 9.5u. 
K. P. YATES AND V. MILLER, Ohio State University.—The 
infra-red band, vs, at 9.54 in the spectrum of methyl fluoride 
was first investigated under high dispersion by Bennett 
and Myer in 1928 who showed it to be of the parallel type. 
The data seemed to indicate, however, that some of the 
rotational lines, particularly in the P branch, were complex 
in structure. This has been investigated by the writers 
under higher dispersion than that available to Bennett and 
Myer. The curve obtained agrees with that of Bennett and 
Myer in general details, but most of the rotational lines 
in both the P and R branches are split up into component 
structure. The effect appears to be related to the con- 
vergence of the lines in the band and it seems reasonable to 
ascribe it as due principally to a Coriolis resonance inter- 
action between the parallel type vibration vs at 1048 cm 
and the perpendicular type vibration vs at 1200 cm=. The 
contour of the band vs of CHF is at present being obtained 
in order to determine in what manner the resonance inter- 
action affects this vibration. 
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7. Raman Spectra of Diisobutylene, Cyclohexene, and 
Dipentene. Forrest F. CLEVELAND, Illinois Institute of 
Technology.—In continuing work on. olefinic hydrocarbons, 
Raman frequencies, relative intensities, and depolarization 
factors have been obtained for diisobutylene, cyclohexene, 
and dipentene. The relative intensities and depolarization 
factors were obtained by use of a Gaertner microdensi- 
tometer. Diisobutylene is a mixture of the two isomers 
2,4,4-trimethyl-1-pentene and 2,4,4-trimethyl-2-pentene, 
the present sample containing only about 15 percent of the 
latter isomer. Raman frequencies and estimated intensities 
for the two isomers have been obtained previously (Rank 
and Bordner), but the present investigation provides the 
first polarization data and more accurate measurements of 
the intensities. Raman frequencies and estimated intensi- 
ties have also been reported for cyclohexene and dipentene, 
but only qualitative polarization data for cyclohexene 
seems to exist in the literature. Dipentene has a ring struc- 
ture identical with that of cyclohexene and it is thus of 
interest to compare the spectra of these two compounds. 
The olefinic frequencies at 1653(447)0.14, 3022(508)0.19, 
and 3058(19) in the cyclohexene spectrum appear at 
1682(219)0.35, 3012(54)P, and 3050(26)P in the spectrum 
of dipentene. The olefinic frequencies previously observed 
at 1303(52)0.64, 1414(108)0.81, 1653(203)0.20, 1674(55)P, 
2978(253)P, and 3077(47)0.70 for 2-methyl-1-heptene, 
which contains the group X YC = CHg, appear at 1236(104)- 
0.63, 1410(155)0.62, 1644(161)0.33, 1658(51)0.60, 2977- 
(790)P, and 3072(32)0.71 for diisobutylene and at 1290- 
(46)0.72, 1376(127)0.70, 1648(213)0.38, 2983(115)P, and 
3078(16)D for dipentene. 


8. Magnetolysis and the Electric Field Around the Mag- 
netic Current. FELIX EHRENHAFT, New York City.—Be- 
tween the vertical cylindrical poles (magnetodes) of an 
electromagnet of soft Swedish iron, whose bases form a 
horizontal gap (pole diam. 8 mm, gap 1-2 mm) acidulated 
water (one percent H2SO, by volume) is decomposed into 
oxygen and hydrogen gas. Hereby the old problem of the 
alchemists about the decomposition of water through the 
magnet has been solved. As long as the two poles immersed 
in the solution are not magnetized we get pure hydrogen, 
but as soon as the two poles are magnetized we get a mix- 
ture of hydrogen and oxygen (about 2-12 percent oxygen). 
The amount of gases evolved magnetically is proportional 
to the magnetic field (first magnetolytic law). Microscopic 
observation shows that magnetically evolved gas bubbles 
carry either a North or South magnetic charge.* Mag- 
netolysis and electrolysis have been superposed. After a 
swarm of positive electrically charged hydrogen bubbles 
are produced, the magnetodes were short-circuited elec- 
trically. Each of these positive charged gas bubbles moves 
in a circle around the gap between the magnetodes, through 
which a constant magnetic current flows, reversing its di- 
rection on reversion of the magnetic field, exactly as a 
single magnetic pole would circulate around the constant 
electric current, reversing its direction with the reversal of 
the electric field (Oersted, Faraday). The magnetic current 
is surrounded by circular electric lines of force (second 
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method of measuring the magnetic current). Photographs 
demonstrating these phenomena will be shown. The experi- 
ments can be seen at C. Zeiss Inc., New York, whom the 
author thanks for their hospitality. 

* Papers by F. Ehrenhaft: Science 1941, 1912; Phys. Rev. 1941, 


1942, 1943; J. Frank. Inst. 1940, 1942; Phil. Sci. 1941—some of which 
with L. Banet. 


9. Combined Electron Diffraction and Microscope Ap- 
paratus. A. F. Presus, The Ohio State University. (Intro- 
duced by Alpheus W. Smith.)—A modified magnetic elec- 
tron microscope is described which permits the observation 
of both highly magnified images and transmission electron 
diffraction patterns of identical microscopic specimen areas. 


10. Investigation of Large Cosmic-Ray Bursts in Iron. 
R. E. Lapp, The University of Chicago.—During the past 
year experiments have been carried out with a Model C 
ionization chamber operated specifically to record cosmic- 
ray bursts. The steel chamber was completely shielded by 
35 cm of iron and bursts with 150 to 3000 particles were 
recorded. Mesotrons producing bursts of such magnitude 
must have energies up to 6X10" ev. Comparison of the 
size frequency distribution curve of large bursts in 35 cm 
of iron was made with the theoretical calculations of 
Christy and Kusaka.* The observed burst frequency was 
in agreement with a mesotron having spin 0 or 4, but was 
in definite disagreement with spin 1. In addition, the re- 
cording system of the ionization chamber was coupled so 
as to record simultaneously coincidences of G-M counter 
coincidence sets placed above the apparatus. The counters 
in each 2-fold set were separated by from 2 to 10 meters 
to detect air showers. Preliminary results indicate that 
bursts under 35 cm of iron are not coincident with air 
showers to an accuracy of less than 10 percent. With the 
chamber completely unshielded and with negligible ma- 
terial in the vicinity, coincidences were observed between 
bursts and air showers. The experiments are being con- 
tinued with 2-, 3-, and 4-fold coincidence sets and experi- 
ments are also in progress with a 12-cm iron shield around 
the ionization chamber. 


*R. F. Christy and S. Kusaka, Phys. Rev. 59, 414 (1941). 


11. Measurements of Cascade Showers Produced by 
Ionizing and Non-Ionizing Radiation in the Stratosphere. 
Juttus TABIN AND MARCEL SCHEIN, University of Chicago. 
—By using a cosmic-ray balloon apparatus, electron- 
produced and photon-produced cascade showers were 
measured below a 2-cm lead plate. The showers were de- 
tected by means of four Geiger-Miiller counters placed 
below and an additional anticoincidence counter placed 
above the lead. The apparatus reached an altitude corre- 
sponding to 2.5 cm of Hg pressure. At a pressure of 3 cm 
Hg (first radiation unit from the top of the atmosphere) 
it was found that 71 percent of the showers below the lead 
were produced by non-ionizing and 29 percent were pro- 
duced by ionizing radiation. 


12. Factors Influencing the Plateau Characteristics Of 
Self-Quenching G-M Counters. W. Spatz, New York Uni- 
versity.—Investigations of factors influencing the starting 
potential and plateau slope of self-quenching G-M counters 
show: I. Impurities such as air or oxygen increase the 
plateau slope and starting potential of argon-alcohol 
counters. II. Operation of the counter causes changes in 
plateau characteristics as follows: (a) Immediately after 
use, the plateau slope increases slightly, (b) With continued 
use the plateau slope becomes steeper. (c) For moderate 
use the counter recovers when inactive, but never to the 
original characteristics. (d) With further use, the counter 
loses its plateau and does not recover. (e) The pressure in 
an argon-alcohol counter was observed to increase as a 
function of the total number of counts recorded. III. For 
argon-alcohol counters the useful life in these experiments 
was found to be about 10°-10'° counts; for argon-methane, 
107-108 counts. The lifetime may be seriously shortened by 
operating at high overvoltages, by allowing breakdown to 
occur, or by using any circuit permitting large currents to 
flow. The observed changes of characteristics are presum- 
ably due to the decomposition of the organic vapor by the 
discharge. 


13. Threshold for the Nuclear Photo-Effect. G. C. 
BALDWIN AND H. W. Kocu, University of Illinois——With 
the method previously described! for controlling the peak 
energy of the x-rays from the 20-Mev betatron, (y, ) reac- 
tions have been observed and their thresholds determined 
in the following elements since the last report: Zn (39 min.) 
11.6+0.4 Mev; N (10 min.), 11.1+0.5 Mev; In (1.1 min.), 
9.5+0.6 Mev; Ag (2.3 min.), 9.34+0.5 Mev; Se (17 min.), 
9.740.5 Mev. These values are upper limits for the 
thresholds. Ag and Se also show strong 25- and 57-minute 
activities, respectively. Analysis of decay curves taken at 
varying energies indicates that these activities have higher 
thresholds than those listed above. An improved deter- 
mination of the O"*(y, 2)O" threshold has been made, giv- 
ing 16.3+0.4 Mev. This and the nitrogen threshold are in 
fair agreement with values predicted from mass data. Chief 
sources of error in the method are (1) gradual change in 
energy calibration with time and (2) low ratio of activity 
to counter background near threshold, making estimation 
of a lower limit difficult. 


1G. C. Baldwin and H. W. Koch, Phys. Rev. 63, 59 (1943). 


14. The Relative Intensities and Characteristic Radia- 
tions of Radioactive Scandium. Cari T. Hispon, M. L. 
Poot, AND J. D. Kursatov, Ohio State University.—A 
reinvestigation of the radioactive isotopes of scandium has 
yielded a new isotope of half-life 3.4 days. This isotope 
emits electrons of 0.46 Mev but no gamma-rays. Other data 
obtained by appropriate bombardments of K, Ca, Sc, Ti, 
and V with alpha-particles, deuterons, protons, slow and 
fast neutrons include the following: 


Isotope Half-life Energy of 8 in Mev Energy of y 
Sc# 3.92 hr. 1.13 positrons 1.65 Mev 
Sc# 3.92 hr. 1.33 positrons 1.33 Mev 
Sc# 2.44 day 1.33 positrons 0.28 and 1.33 
Sc# 1.83 day 0.57 electron 1.33 Mev 
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No evidence was found for the existence of the 13.2-day 
nor the 2.62-day periods previously assigned to Sc and 
Sc*’, respectively. The Ca(p, 2) reaction yielded Sc** and 
Sc but not Sc“. With deuteron bombardment of calcium 
the most abundantly produced isotope is Sc; the isomers 
of Sc“, however, are not produced by this reaction in equal 
proportions as they are by the Sc(m, 2”) and K(a, m) reac- 
tions. The scandium oxide used was of high purity and 
kindly furnished by Professor Quill of the University of 
Kentucky. 


15. The Artificial Radioactivities of Cerium. M. L. Poot 
AND J. D. Kursatov, Ohio State University—Four new 
radioactive isotopes have been produced and assigned posi- 
tions in cerium and praseodymium. Ce"™®, which has a half- 
life of 140 days, is obtained by the reactions Ba!*"(a, n) 
and La!**(d, ). The fact that gamma-rays only (of 0.21 
Mev) are emitted by this isotope is consistent with the 
isomeric assignment. A 30-day activity is assigned to 
Ce"! and has been obtained by the reactions Ba™*(a, m), 
Ce°(d, p), Cel°(m, y), Ce™*(m, 2m), and Pr™(n, p). The iso- 
tope emits beta-particles of 0.65 Mev and gamma-rays of 
0.2 Mev. Ce™! is of particular interest because of the un- 
usual arrangement of stable isotopes in the cerium region. 
A 36-hour activity; assigned to Ce, has been produced by 
the reactions Ce™(d, p) and Ce™*(n, 7). This substance 
forms a chain reaction by disintegrating into Pr'**, Pr in 
turn disintegrates with a half-life of 13.5 days by the emis- 
sion of 0.95 Mev beta-particles. No gamma-rays are 
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emitted. The study of the above radioactive isotopes and 
the clarification of certain meagerly detailed published data 
on radioactive cerium has been made possible and expedited 
by improved methods used in the separation of cerium from 
praseodymium and lanthanum. 


16. Progress Report on the Radioactivities in the 
Illinium Region. J. D. KuRBATOV AND M. L. Poo, Ohio 
State University —A continuation of comparative studies 
of the radioactivities of Pr, Nd, Il, Sm, and Eu has led to 
the following results. The nuclear reaction Nd(p, ”) pro- 
duces the periods 2.7 hours and 5.3 days. Both activities 
disintegrate with the emission of gamma-radiation and 
2-Mev electrons; this suggests an isomeric position for the 
two isotopes in. illinium. The nuclear reaction Nd(d, n) 
produces, in addition to the above, a period of 16 days. The 
disintegration takes place with the emission of gamma- 
radiation and 1.7-Mev electrons. The same activation 
yields a long period of about one year. Gamma-rays of 0.7 
Mev constitute the principal radiations emitted. The acti- 
vation of samarium with deuterons results in several 
periods. The Sm(d, m) reaction yields a 40-day activity and 
a one-year activity. Both are associated with europium and 
disintegrate mainly by the emission of gamma-rays. Ab- 
sorption measurements of the gamma-rays of this one-year 
period presents results comparable to those of the one year 
period obtained by deuteron bombardment of Nd. Thus a 
chain reaction of four elements is suggested. 
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